## Glimpse of the Box Experiment dataset:
## Rows: 2,795
## Columns: 20
## $ Date <dttm> 2022-09-27, 2022-09-27, 2022-09-27, 2022-09-27,…
## $ Time <dttm> 1899-12-31 09:47:50, 1899-12-31 09:50:07, 1899-…
## $ Data <chr> "Box Experiment", "Box Experiment", "Box Experim…
## $ Group <chr> "Baie Dankie", "Baie Dankie", "Baie Dankie", "Ba…
## $ GPSS <chr> "-28.010549999999999", "-28.010549999999999", "-…
## $ GPSE <chr> "31.191050000000001", "31.191050000000001", "31.…
## $ MaleID <chr> "Nge", "Nge", "Nge", "Nge", "Nge", "Nge", "Nge",…
## $ FemaleID <chr> "Oerw", "Oerw", "Oerw", "Oerw", "Oerw", "Oerw", …
## $ `Male placement corn` <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
## $ MaleCorn <dbl> 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, …
## $ FemaleCorn <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
## $ DyadDistance <chr> "2m", "2m", "1m", "1m", "0m", "0m", "0m", "0m", …
## $ DyadResponse <chr> "Tolerance", "Tolerance", "Tolerance", "Toleranc…
## $ OtherResponse <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
## $ Audience <chr> "Obse; Oup; Sirk", "Obse; Oup; Sirk", "Oup; Sirk…
## $ IDIndividual1 <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
## $ IntruderID <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, "Sey…
## $ Remarks <chr> NA, NA, "Nge box did not open because of the bat…
## $ Observers <chr> "Josefien; Michael; Ona; Zonke", "Josefien; Mich…
## $ DeviceId <chr> "{7A4E6639-7387-7648-88EC-7FD27A0F258A}", "{7A4E…
I am now using the View function to have a sight on the entire dataset and glimpse to display a summary of my dataset
I have 20 variables (here columns) and 2795 trials (here rows)
I will now make a brief summary of each variables and their use before creating a new dataframe (df) with my variables of interest that I will call Bex
The highlighted variables are the ones I will use for Bex. I will then clean the data before heading to the statistical analysis and the interpretation of the results
Date : “Date” is in a POSIXct format which is appropriate for the display of time
Time : “Time is coded” in a POSIXct format
Data : chr “Data” is coded as character
Group : chr The data is coded in r as a character
GPSS : num “GPSS” is coded as numerical
GPSE : num “GPSE” is coded in as numerical
MaleID : chr “MaleID” is coded as character
FemaleID : chr “FemaleID” is coded as character
Male placement corn: dbl “Male placement corn is coded in r as double
It gives the amount of corn given to the male of the dyad before the trials
Within a session it happened that we gave more placement corn to attract the monkeys again to the boxes. This lead to an update of the number in the same session. The number found at the end of the session is the total placement corn an individual has received
I will fuse this column with male corn as the data has been separated between these two variables. This is due to a mistake when creating the original box experiment form in cybertracker
This variable could be related to the level of motivation of a monkey but as it is not directly related to my hypothesis I may not use this column. I will re-consider the use of this column later on
In regards of this possibility I will change the format of the variable to numerical
MaleCorn : dbl “MaleCorn” is coded in r as double
FemaleCorn : dbl The data is coded in r as double
DyadDistance : chr The data is coded in r as character
DyadResponse : chr The data is coded in r as character
Create a table with each combination existing
Decide what is more important
Ex:
OtherResponse : chr “The data”OtherResponse” is coded as character
Audience : chr “Audience” is in r as character
IDIndividual1 : chr “IDIndividual1” is coded in r as character
IntruderID : chr “IndtruderID” is coded as character
Remarks : chr The data is coded in r as character
Observers :chr The data is coded in r as character
DeviceID :chr “The data”DeviceID” is coded in r as character
Since I do not want to work with the whole dataset, I’m gonna select the variables of interest using the function select
I will keep Time, Date, Group, MaleID, FemaleID, MaleCorn, Male placement corn, FemaleCorn, DyadDistance, DyadResponse, OtherResponse, Audience, IDIndividual1, IntruderID, Remarks
## Rows: 2,795
## Columns: 15
## $ Time <dttm> 1899-12-31 09:47:50, 1899-12-31 09:50:07, 1899-…
## $ Date <dttm> 2022-09-27, 2022-09-27, 2022-09-27, 2022-09-27,…
## $ Group <chr> "Baie Dankie", "Baie Dankie", "Baie Dankie", "Ba…
## $ MaleID <chr> "Nge", "Nge", "Nge", "Nge", "Nge", "Nge", "Nge",…
## $ FemaleID <chr> "Oerw", "Oerw", "Oerw", "Oerw", "Oerw", "Oerw", …
## $ MaleCorn <dbl> 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, …
## $ `Male placement corn` <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
## $ FemaleCorn <dbl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
## $ DyadDistance <chr> "2m", "2m", "1m", "1m", "0m", "0m", "0m", "0m", …
## $ DyadResponse <chr> "Tolerance", "Tolerance", "Tolerance", "Toleranc…
## $ OtherResponse <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
## $ Audience <chr> "Obse; Oup; Sirk", "Obse; Oup; Sirk", "Oup; Sirk…
## $ IDIndividual1 <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
## $ IntruderID <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, "Sey…
## $ Remarks <chr> NA, NA, "Nge box did not open because of the bat…
## Number of rows with common NAs in MaleCornOld and 'Male placement corn': 1499
## Number of occurrences of 0 in MaleCorn: 1499
## Number of remaining NA values in MaleCorn: 0
I have found 1499 NA in common between MaleCornOld and ‘male placement corn’, 1609 NA in Male placement corn and 2685 in MaleCorn old
For the merge of MaleCornOld and Male placement corn, I used different conditions: 1.In this code, a new variable MaleCorn is created. If there is a missing value in Male placement corn, it takes the corresponding value from MaleCornOld; otherwise, it takes the value from Male placementcorn. 2.If there are no value in both MaleCornOld and Male placement corn (NA,NA) for a given row, I would like the code to display 0 as it means that no placement was given
In this way, I should not loose any data, minimize the mistakes and already transform the NA’s of this variable into a number which will remove the remaining NA’s which are meant to be 0
After the merge I found that there were no NA’s remaining in the “New” Male Corn and that 1499 0’s where found in the column which corresponds to the amount of common NA’s found previously between the “Old” Male Corn and male placement corn
## Number of remaining NA values in FemaleCorn: 0
Now in order to see where are located the missing points in the data, I’m going to print the variables with and without NA’s
The function sapply is used to apply the function sum for NA’s to each column of the data frame, so each variable
## Variables with Missing Data:
| x | |
|---|---|
| MaleID | 19 |
| FemaleID | 60 |
| DyadDistance | 33 |
| DyadResponse | 47 |
| OtherResponse | 2758 |
| Audience | 924 |
| IDIndividual1 | 2143 |
| IntruderID | 2737 |
| Remarks | 2181 |
## Variables with No Missing Data:
| x | |
|---|---|
| Time | 0 |
| Date | 0 |
| Group | 0 |
| FemaleCorn | 0 |
| MaleCorn | 0 |
We can see that out of the 14 variables we have in Bex we have 9 variables with missing data which are Male ID, Female ID, DyadDistance, DyadResponse, OtherResponse, Audience, IDIndividual1, IntruderID, Remarks: I will proceed to clean these variables one by one
MaleID 19
FemaleID 60
DyadDistance 33
DyadResponse 47
OtherResponse 2758
Audience 924
ID Individual1 2143
IntruderID 2737
Remarks 2181
Before making treating the NA’s in the dataset I will make a backup of the data at this point:
Since most of the time we did not have any remarks it is understandable that this variable contains 2181 NA’s out of 2795 rows
I will first transform every missing data in the column Remark into No Remarks and then check that the amount of “No remarks” found
After the changes we can effectively see that we have 2181 “No Remarks” and we have no missing data left in that column, I will treat this column by hand once all the NA’s have been removed from the dataset
## Number of 'No Remarks' in the 'Remarks' column: 2181
##
## No Remarks Remarks
## 2181 614
## Number of 'No Intrusion' in the 'Intruder ID' column after replacement: 2737
## Number of NAs replaced in IDIndividual1: 2143
## Number of remaining NA values in IDIndividual1: 0
## Number of changes made in 'Audience': 924
## Remaining NA values in 'Audience': 0
## Number of changes made in 'OtherResponse': 2758
## Remaining NA values in 'OtherResponse': 0
## [1] "1899-12-31 09:47:50 UTC" "1899-12-31 09:50:07 UTC"
## [3] "1899-12-31 09:53:11 UTC" "1899-12-31 09:54:28 UTC"
## [5] "1899-12-31 09:55:19 UTC" "1899-12-31 09:56:56 UTC"
## [1] "09:47:50" "09:50:07" "09:53:11" "09:54:28" "09:55:19" "09:56:56"
## Warning: NAs introduits lors de la conversion automatique
## # A tibble: 69 × 16
## Time Date Group MaleID FemaleID FemaleCorn DyadDistance
## <chr> <dttm> <chr> <chr> <chr> <dbl> <dbl>
## 1 12:09:34 2022-09-27 00:00:00 Baie Da… Xia Piep 7 NA
## 2 12:13:28 2022-09-27 00:00:00 Baie Da… Xia Piep 7 NA
## 3 16:02:32 2022-09-15 00:00:00 Ankhase Sho Ginq 6 NA
## 4 10:46:33 2023-08-17 00:00:00 Baie Da… Xia Piep 0 NA
## 5 09:30:17 2023-07-29 00:00:00 Baie Da… Xin Ouli 0 NA
## 6 12:08:51 2023-07-11 00:00:00 Baie Da… Xia Piep 0 NA
## 7 13:30:07 2023-06-29 00:00:00 Baie Da… Sey Sirk 0 NA
## 8 09:54:24 2023-06-27 00:00:00 Ankhase Sho Ginq 0 NA
## 9 10:13:56 2023-06-23 00:00:00 Ankhase Sho Ginq 0 NA
## 10 09:39:04 2023-06-15 00:00:00 Ankhase Sho Ginq 2 NA
## # ℹ 59 more rows
## # ℹ 9 more variables: DyadResponse <chr>, OtherResponse <chr>, Audience <chr>,
## # IDIndividual1 <chr>, IntruderID <chr>, Remarks <chr>, MaleCorn <dbl>,
## # Intrusion <dbl>, AmountAudience <dbl>
## Number of NA values in DyadDistance column (using second approach): 69
## Rows with NA values in DyadDistance column: 24, 27, 95, 492, 744, 971, 1113, 1130, 1164, 1261, 1341, 1396, 1491, 1583, 1683, 1693, 1717, 1718, 1719, 1724, 1725, 1739, 1755, 1756, 1757, 1764, 1779, 1782, 1792, 1799, 1800, 1840, 1841, 1868, 1869, 1888, 1891, 1892, 1896, 1911, 1912, 1915, 1918, 1919, 1952, 1953, 1958, 1980, 1981, 1984, 1986, 1996, 2000, 2009, 2054, 2104, 2105, 2191, 2233, 2234, 2287, 2437, 2569, 2579, 2580, 2643, 2676, 2709, 2729
We have 69 missing values in DyadDistance. I will look at each row in it’s context as the actual distance of the box was always dependent of the previous trials. I will start with the bigger number as for now the oldest trial is at the last row while the closest one is in row 1.
Now that I have looked at each missing line and saw which ones to keep, I decided to create a new variable called Distance. I will also to create a new variable called No trial.
For the variable Distance I will replace each row where there was missing data with a value and I will delete the ones where no values could be assigned. This will allow me to have no missing data and find a number to each trial that has been done
Before making the changes i’m gonna make a backup called BackupbeforeDistanceNA
## Number of NA's in DyadDistance after replacements and deletions: 1
## Data size after deletions: 2748
## Row index with NA in DyadDistance: 1925
*It seems that there is still the row 1925 with an NA in DyadDistance
## Row index with NA in DyadDistance:
*In this modification, I added a check to see if the columns Dyadistance and Distance already exist in your dataframe (Bex). If they do, it prints a message saying that the modification has already been applied, and no changes are made. If they don’t exist, it proceeds with the modifications. This way, running the code multiple times won’t cause redundant changes.
Before cleaning Female and Male ID, here is a list of every dyad of the box experiment and their respective groups. This will help us find the missing names when only one individual is missing out of the duo (either male or female):
Sirk & Sey - BD
Ouli & Xin - BD
Piep & Xia - BD
Oerw & Nge - BD
Oort & Kom - BD
Ginq & Sho - AK
Ndaw & Buk - Ak
Xian & Pom - AK
Guat & Pom - Ak
Note that the 4 letter codes correspond to the femaleID, the 3 letter codes to the males ID and the 2 letter codes to the group name of the monkeys
I need to check where are the NA’s in both FemaleID and Male ID by looking at the rows where data is missing. Since every trial was made with a Dyad and never with an single individual, treating these two columns together makes more sense. If both individuals are missing I may have to delete the row.
## Row numbers with missing values in FemaleID: 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 1693 1694 1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710 1808 1809 1810 1811 1812 1813 1814 1815 1816 1817 1818 1819 1820 1884 1885 2619 2620 2621 2622 2623 2624 2625 2626 2627 2628 2629
## Number of missing values in FemaleID: 59
## Row numbers with missing values in MaleID: 1693 1694 1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710
## Number of missing values in MaleID: 18
## Number of rows with missing values in both FemaleID and MaleID: 18
## Row numbers with missing values in both FemaleID and MaleID: 1693, 1694, 1695, 1696, 1697, 1698, 1699, 1700, 1701, 1702, 1703, 1704, 1705, 1706, 1707, 1708, 1709, 1710
## Number of missing values in FemaleID not in MaleID: 41
## Row numbers with missing values in FemaleID not in MaleID: 865, 866, 867, 868, 869, 870, 871, 872, 873, 874, 875, 876, 877, 878, 879, 1808, 1809, 1810, 1811, 1812, 1813, 1814, 1815, 1816, 1817, 1818, 1819, 1820, 1884, 1885, 2619, 2620, 2621, 2622, 2623, 2624, 2625, 2626, 2627, 2628, 2629
FemaleID has 41 NA’s while they are 18 NA’s in Male ID
In these missing data, we have 18 NA’s that are in common between FemaleID and MaleID which represents the totality of the missing values in MaleID
All the missing data in MaleID are found in consecutive rows, from row 1693 to row 1710 and are from the group Noha (NH) on the 19th of april 2023. We can also see that trials had bee made in the same day, and looking at the time of the experiment, the previous trials made and the audience we can see that these NA’s in female and male ID we can asses that the individuals involved were Xian for the female ID and Pom for the MaleID. I will thus replace these values using a condtion. These NA’s in Noha (Trial 1693 to 1710) are the only NA’s that MaleID has and are the only NA’s of female ID in Noha. I will thus replace every NA of MaleID NA in Noha with Pom and every Female ID NA in Noha with Xian
## Number of remaining NA values in MaleID after replacement: 0
## Number of remaining NA values in FemaleID after replacement: 41
## Number of rows with missing values in both MaleID and FemaleID after replacement: 0
In order to clean FemaleID, I will use the data from the now complete MaleID. I will use conditions stating that depending which name is found in MaleID when there is an NA in FemaleID, a certain name will have to replace the NA in female ID
Before automating the process I will check manually the data to see if they are any exceptions or mistakes
## Rows with missing values in FemaleID:
## # A tibble: 41 × 16
## Time Date Group MaleID FemaleID FemaleCorn DyadDistance
## <chr> <dttm> <chr> <chr> <chr> <dbl> <dbl>
## 1 09:31:55 2023-07-22 00:00:00 Ankhase Buk <NA> 7 1
## 2 09:33:14 2023-07-22 00:00:00 Ankhase Buk <NA> 7 1
## 3 09:34:07 2023-07-22 00:00:00 Ankhase Buk <NA> 7 0
## 4 09:34:51 2023-07-22 00:00:00 Ankhase Buk <NA> 7 0
## 5 09:36:59 2023-07-22 00:00:00 Ankhase Buk <NA> 7 0
## 6 09:38:13 2023-07-22 00:00:00 Ankhase Buk <NA> 7 1
## 7 09:39:26 2023-07-22 00:00:00 Ankhase Buk <NA> 7 0
## 8 09:41:11 2023-07-22 00:00:00 Ankhase Buk <NA> 0 0
## 9 09:42:17 2023-07-22 00:00:00 Ankhase Buk <NA> 0 0
## 10 09:44:06 2023-07-22 00:00:00 Ankhase Buk <NA> 0 1
## # ℹ 31 more rows
## # ℹ 9 more variables: DyadResponse <chr>, OtherResponse <chr>, Audience <chr>,
## # IDIndividual1 <chr>, IntruderID <chr>, Remarks <chr>, MaleCorn <dbl>,
## # Intrusion <dbl>, AmountAudience <dbl>
If there is NA in femaleID, we will replace the value with - Sirk if MaleID is Sey - Ouli if MaleID is Xin - Piep if MaleID is Xia - Oerw if MaleID is Nge - Oort if MaleID is Kom - Ginq if MaleID is Sho - Ndaw if MaleID is Buk
## # A tibble: 20 × 3
## MaleID FemaleID Count
## <chr> <chr> <int>
## 1 Xia Piep 576
## 2 Sey Sirk 557
## 3 Kom Oort 338
## 4 Sho Ginq 278
## 5 Pom Xian 259
## 6 Buk Ndaw 245
## 7 Xin Ouli 159
## 8 Nge Oerw 153
## 9 Piep Xia 35
## 10 Oort Kom 29
## 11 Ouli Xin 27
## 12 Oerw Nge 19
## 13 Sirk Sey 17
## 14 Buk <NA> 15
## 15 Sey <NA> 13
## 16 Nge <NA> 11
## 17 Buk Ginq 6
## 18 Pom Guat 5
## 19 Xin Oort 4
## 20 Kom <NA> 2
## Number of NA values in MaleID: 0
## Number of NA values in FemaleID: 0
## Rows with missing values in DyadResponse: 871, 1163, 1219, 1339, 1579, 1888, 1962
## Lines with missing values in DyadResponse:
## # A tibble: 7 × 16
## Time Date Group MaleID FemaleID FemaleCorn DyadDistance
## <chr> <dttm> <chr> <chr> <chr> <dbl> <dbl>
## 1 09:39:26 2023-07-22 00:00:00 Ankhase Buk Ndaw 7 0
## 2 10:13:56 2023-06-23 00:00:00 Ankhase Sho Ginq 0 4
## 3 08:34:45 2023-06-17 00:00:00 Baie Dan… Kom Oort 3 2
## 4 08:54:12 2023-06-09 00:00:00 Baie Dan… Xia Piep 1 0
## 5 13:35:08 2023-05-03 00:00:00 Baie Dan… Kom Oort 5 3
## 6 13:27:30 2023-01-18 00:00:00 Ankhase Buk Ndaw 5 4
## 7 08:36:49 2022-12-13 00:00:00 Baie Dan… Kom Oort 8 4
## # ℹ 9 more variables: DyadResponse <chr>, OtherResponse <chr>, Audience <chr>,
## # IDIndividual1 <chr>, IntruderID <chr>, Remarks <chr>, MaleCorn <dbl>,
## # Intrusion <dbl>, AmountAudience <dbl>
Row 871: The previous row was tolerance at 1m and the next tolerance at 0 which means that the row 871 should be Tolerance for DyadResponse
Row 1163: The value can not be found from the other rows so I will delete row 1163
Row 1219: The previous row was not approaching at 2m and the next is tolerance at 2m and tolerance at 1m, which means that the row 1219 should be Tolerance for DyadResponse
Row 1339: The previous row was tolerance at 0m while the next one was tolerance at 0m, which means that the tow 1339 should be Tolerance for DyadResponse
Row 1579: The value can not be found from the other rows so I will delete row 1579
Row 1888: The value can not be found from the other rows so I will delete row 1888
Row 1962: The value can not be found from the other rows so I will delete row 1962
## Number of remaining NA values in DyadResponse: 0
## Final check of NA values in Bex:
## Time Date Group MaleID FemaleID
## 0 0 0 0 0
## FemaleCorn DyadDistance DyadResponse OtherResponse Audience
## 0 0 0 0 0
## IDIndividual1 IntruderID Remarks MaleCorn Intrusion
## 0 0 0 0 0
## AmountAudience
## 0
Since I have removed all the missing data from the different columns, I now have to correct potential mistakes that can be found and create new variables to be able to manipulate better my data.
Since the column remarks contains corrections and additional information, I will treat it now
Before that lets check how many remarks we have in our dataset, how many of the main keywords we can find and make a visual representation of it
## Number of 'No Remarks' entries: 2139
## Number of actual remarks entries: 603
## Total number of keyword occurrences in the Barplot: 822
## Glimpse of the Bex Before treating Remarks:
## Rows: 2,742
## Columns: 16
## $ Time <chr> "09:47:50", "09:50:07", "09:53:11", "09:54:28", "09:55:…
## $ Date <dttm> 2022-09-27, 2022-09-27, 2022-09-27, 2022-09-27, 2022-0…
## $ Group <chr> "Baie Dankie", "Baie Dankie", "Baie Dankie", "Baie Dank…
## $ MaleID <chr> "Nge", "Nge", "Nge", "Nge", "Nge", "Nge", "Nge", "Nge",…
## $ FemaleID <chr> "Oerw", "Oerw", "Oerw", "Oerw", "Oerw", "Oerw", "Oerw",…
## $ FemaleCorn <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 7, 7, 7, 7…
## $ DyadDistance <dbl> 2, 2, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 2, 2, 1, 1, 0, 0…
## $ DyadResponse <chr> "Tolerance", "Tolerance", "Tolerance", "Tolerance", "To…
## $ OtherResponse <chr> "No Response", "No Response", "No Response", "No Respon…
## $ Audience <chr> "Obse; Oup; Sirk", "Obse; Oup; Sirk", "Oup; Sirk", "Sir…
## $ IDIndividual1 <chr> "No individual", "No individual", "No individual", "No …
## $ IntruderID <chr> "No Intrusion", "No Intrusion", "No Intrusion", "No Int…
## $ Remarks <chr> "No Remarks", "No Remarks", "Nge box did not open becau…
## $ MaleCorn <dbl> 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3…
## $ Intrusion <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0…
## $ AmountAudience <dbl> 3, 3, 2, 1, 2, 2, 2, 1, 1, 2, 6, 6, 3, 2, 2, 2, 2, 2, 2…
## [1] "/Users/maki/Desktop/Master Thesis/BEX 2223 Master Thesis Maung Kyaw/IVPToleranceBex2223"
Before treating all the data in the Remarks I will create a few columns to redistribute information
Also whenever i will have treated a remark, i will replace it with “Treated”. And if I have to delete the row I’ll write “Delete”. After re importing the data I will make a count of these changes to see if I still have the correct amount of cells and changes that have been done
a.Context: BoxMalfunction, BoxOpenedBefore, NoExperiment, Agonistic, Guat;Ap;Xian, CornLeak, BetweenGroupEncounter, ContactCalling,
b.SpecialBehaviour Oerw;Vo;Exp, Sey;Ap;AfterOpen, Oerw;Vo;Exp,Nge;Vo;Exp, Sirk;ApAfter30, Sirk;Av;Oerw, Oerw;Lo,Sey;Sf;Oort,Oort;At;Kom, Kom;Ap;AfterOpen, Sey;Ch,Sirk, Xin;Hesitation. Xia;Sf;Piep, Pom;Sf;Xian,Kom;Sf;Oort, Sey;Sf;Sirk, Xia;Sf;Piep,Piep;Sf,XIa, Oort;At;Kom, Sey;Rt;Sho;Ap, Sho;Rt;Ginq;Ap, Buk;Sf;Ndaw, Sho;Rt;Ndaw;Ap, Oort;Sf;Kom, Ginq;Sho;Ap;After30, Ndaw;Sc,Buk;Sf, Ndaw;Ap;After30, Kom;Ap;After30, Xia;Piep;Ap;After30, Pom;Bi;Xian, Sho;Ndaw;Av;Buk, Kom;Sf;Oort, Kom;St;Oort,Oort;St;Kom, Sey;Hi;Sirk, Obse;Ap;Piep;Av,Piep;Sf;Xia, Sirk;ApWhenPartnerLeft, Sey;Hh;Sirk, Xia;Sf;Piep;Sc, Xia;Piep;ShareFood, Piep;Ap;After30,Xia;Mu;Piep, Oort;St;Sirk;Ja,Sey;Sf, Pom;Sf;Xian, Ndaw;ApWhenPartnerLeft, Xian;At;Pom,Gaya;Su, Xian;Sf;Pom, Xian;Hesitation, Xia;ApWhenPartnerLeft,Sirk;Hesitation, Ginq;Hesitation, Sey;Ap;Kom;Av, Oort;Sc;Kom, Xian; Pom, Pom;Ap;Xian, Pom;Ap;Xian,Xian;Rt, Sey;Ap;Sirk;Rt, Sey;St;Sirk;Ig, Xia;Asf;Piep, Piep;ApWhenPartnerLeft, Sho;Ap;After30, Ginq;ApWhenPartnerLeft, Pom;Sf;Xian;Sf;Pom, Xian;ApWhenPartnerLeft, Piep;Ch;Sirk, Sey;St;Sirk, Ndaw;Ap;After30, Xian;Ap;After30, Xian;St;Prai, Pom;Sf;Xian;Vc, Kom;Ap;After30, Kom;ApproachWithPartner, Oort;ApWhenPartnerLeft, Sho;Ap;After30,Ginq;Ap;After30, Ginq;ApproachWithPartner, Ndaw;Hesitation, Oerw;Hesitation, Oerw;ApWhenPartnerLeft. Piep;Ap;After30, Sirk;Ap;After30, Xia;Ap;After30, Ouli;Gr;BBOuli, Oerw;Ap;After30, Sirk;Hesitation, Sey;Ap;Sirk;Av, Ouli;Ap;Xia;Av, Xin;Ap;After30, Sho;Sf;Ginq;Sc, Xia;ApWhenPartnerLeft, Sey;Ap;Sirk;Ja, Nge;Oerw;ShareFood, Nge;Ap;Oerw;Oerw;At,Obse;At;Nge,
c.GotCorn: No;Nge, No;Piep, No;Xian, No;Oort, No;Sirk, No;Kom, No;Ndaw, No;Kom, No;Oort, No;Xia, No;Buk, No;Sho, No;Sey,No;Piep, No;Ginq
IntruderID: Sey, Oerw, Guat, Kom, Gris, Sho, Oerw; Ouli, Guat; Gri, Xop, Obse, Oort, Obse; Sey, Ginq; Ghid, Xia, Grif, Sey, Gree; Gran, Godu; Gub, Gran, Oerw; Nak, Ghid, Buk, Oup
DyadDistance: 6, 7, 8, 9 , 1
Audience: UnidentifiedAudience, Ouli; Riss, Gris, Sey, Sey; Piep; Sirk, Oup Ome
IDIndividual1: Piep, Oort; Kom, Ndaw; Buk, Sho; Ginq, Ndaw, Buk, Xian, Pom, Oort; Kom, Buk; Ndaw, Sirk; Sey, Xin; Ouli, Oerw; Nge
DyadResponse: Tolerance, Not approaching; Losing interest, Losing interest; Intrusion
## Rows: 2,742
## Columns: 19
## $ Time <chr> "09:47:50", "09:50:07", "09:53:11", "09:54:28", "09:5…
## $ Date <dttm> 2022-09-27, 2022-09-27, 2022-09-27, 2022-09-27, 2022…
## $ Group <chr> "Baie Dankie", "Baie Dankie", "Baie Dankie", "Baie Da…
## $ MaleID <chr> "Nge", "Nge", "Nge", "Nge", "Nge", "Nge", "Nge", "Nge…
## $ FemaleID <chr> "Oerw", "Oerw", "Oerw", "Oerw", "Oerw", "Oerw", "Oerw…
## $ FemaleCorn <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 7, 7, 7,…
## $ DyadDistance <dbl> 2, 2, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 2, 2, 1, 1, 0,…
## $ DyadResponse <chr> "Tolerance", "Tolerance", "Tolerance", "Tolerance", "…
## $ OtherResponse <chr> "No Response", "No Response", "No Response", "No Resp…
## $ Audience <chr> "Obse; Oup; Sirk", "Obse; Oup; Sirk", "Oup; Sirk", "S…
## $ IDIndividual1 <chr> "No individual", "No individual", "No individual", "N…
## $ IntruderID <chr> "No Intrusion", "No Intrusion", "No Intrusion", "No I…
## $ Remarks <chr> "No Remarks", "No Remarks", "Treated", "Treated", "No…
## $ MaleCorn <dbl> 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,…
## $ Intrusion <dbl> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0,…
## $ AmountAudience <dbl> 3, 3, 2, 1, 2, 2, 2, 1, 1, 2, 6, 6, 3, 2, 2, 2, 2, 2,…
## $ Context <chr> "NoContext", "NoContext", "BoxMalfunction", "BoxOpene…
## $ SpecialBehaviour <chr> "NoSpecialBehaviour", "NoSpecialBehaviour", "Oerw;Vo;…
## $ GotCorn <chr> "Yes", "Yes", "No;Nge", "Yes", "Yes", "Yes", "Yes", "…
## Number of NA entries in Context: 0
## Number of NA entries in SpecialBehaviour: 0
## Number of NA entries in GotCorn: 0
## Number of NA entries in BexClean: 0
Time : I considered looking at the time sections in which we did the expermiment. I will thus look at the time ranges (max and min in the day / latest and earliest time) before separating the day in different sections to have an idea in which part of the day most of the experiments occured. This will not be used in my analysis, but if I wanted to, I could interesting to compare the amount of experimentations made per day and have a line indicating the time of sunrise.
The Minimum Time in the dataset is 06:03:26* while the Maximum Time is at 16:36:59
In my box experiment I have this variable called time that tells me when the experiment was done. I don’t think I need this information per se. I was wondering if it could be easy and interesting to see from when to when the time occurs and then separate this time in a few sections like early, monring, morning, miday, afternoon, end of the day
a.6 to 8 : Early morning b.8 to 10: Morning c.10 to 12: Noon d.12 to 14: Afternoon e.14 to 17: End of the day
Last, I want to create a variable called Hour that will take the value in Time and round it to the hour in which it is ex: from 06:00 to 06:59 -> 6, from 07:00 to 07:59 -> 7 etc…
This will allow me to see when most of the trials occured with more detail and I will be to see in which hour most of the trial happened. Nevertheless Period will be better for an improved readability
## Unique Female IDs: Sirk Ginq Piep Oerw Xin Ndaw Xia Sey Ouli Nge Oort Xian Guat Kom
## Unique Male IDs: Sey Sho Xia Nge Ouli Buk Piep Sirk Xin Oerw Kom Pom Oort
| Buk | Kom | Nge | Oerw | Oort | Ouli | Piep | Pom | Sey | Sho | Sirk | Xia | Xin | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Ginq | 6 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 277 | 0 | 0 | 0 |
| Guat | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 5 | 0 | 0 | 0 | 0 | 0 |
| Kom | 0 | 0 | 0 | 0 | 29 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| Ndaw | 259 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| Nge | 0 | 0 | 0 | 19 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| Oerw | 0 | 0 | 164 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| Oort | 0 | 337 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 4 |
| Ouli | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 159 |
| Piep | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 575 | 0 |
| Sey | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 17 | 0 | 0 |
| Sirk | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 570 | 0 | 0 | 0 | 0 |
| Xia | 0 | 0 | 0 | 0 | 0 | 0 | 35 | 0 | 0 | 0 | 0 | 0 | 0 |
| Xian | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 259 | 0 | 0 | 0 | 0 | 0 |
| Xin | 0 | 0 | 0 | 0 | 0 | 27 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
## Unique Dyads: Sey Sirk Sho Ginq Xia Piep Nge Oerw Xin Ouli Buk Ndaw Buk Ginq Kom Oort Pom Xian Pom Guat Xin Oort
| Var1 | Freq |
|---|---|
| Buk Ginq | 6 |
| Buk Ndaw | 259 |
| Kom Oort | 366 |
| Nge Oerw | 183 |
| Pom Guat | 5 |
| Pom Xian | 259 |
| Sey Sirk | 587 |
| Sho Ginq | 277 |
| Xia Piep | 610 |
| Xin Oort | 4 |
| Xin Ouli | 186 |
## Unique Male-Female Combinations:
## # A tibble: 11 × 2
## Male Female
## <chr> <chr>
## 1 Sey Sirk
## 2 Sho Ginq
## 3 Xia Piep
## 4 Nge Oerw
## 5 Xin Ouli
## 6 Buk Ndaw
## 7 Buk Ginq
## 8 Kom Oort
## 9 Pom Xian
## 10 Pom Guat
## 11 Xin Oort
## [1] "Wrong Rows:"
## [1] 613 614 615 616 617 931 2710 2711 2712 2713
## [1] "Wrong Dyads:"
## [1] "Buk Ginq" "Buk Ginq" "Buk Ginq" "Buk Ginq" "Buk Ginq" "Buk Ginq"
## [7] "Xin Oort" "Xin Oort" "Xin Oort" "Xin Oort"
They are a 10 wrong dyads that I will have to identify in the dataset and manually correct, those wrong dyads to change and identify are: -Buk Ginq - 6 occurences -Xin Oort - 4 occurences
I will change the occurences of Buk Ginq to Sho Ginq for row 613 to 617 and row 931. I know these trials are with Sho Ginq because the comments mentioned Sho in them while Male(ID) gave Buk which was a mistake
For the rows from 2710 to 2713 since, Ouli is in the audience it is unlikely that we had trials with the dyad Xin Ouli. Also I think they are little chances that the names of both individuals were entered wrong. I will replace these occurences where we had Xin Oort by Kom Oort
I thus want Buk to be replaced in male ID in rows 613 to 617 and row 913 with Sho and, Xin to be replaced by Kom in rows 2710 to 2713 in Male ID before updating Dyad
If Rows 613 to 617 and 931 are coded with Buk for MaleId and Ginq for FemaleId replace Male by Sho
## Rows to correct Sho Ginq:
## [1] 613 614 615 616 617 931
## Rows to correct Kom Oort:
## [1] 2710 2711 2712 2713
## Wrong Rows After Correction:
## integer(0)
## Wrong Dyads After Correction:
## character(0)
## All dyads are now correct.
## Unique Dyads after correction: Sey Sirk Sho Ginq Xia Piep Nge Oerw Xin Ouli Buk Ndaw Kom Oort Pom Xian Pom Guat
| Var1 | Freq |
|---|---|
| Buk Ndaw | 259 |
| Kom Oort | 370 |
| Nge Oerw | 183 |
| Pom Guat | 5 |
| Pom Xian | 259 |
| Sey Sirk | 587 |
| Sho Ginq | 283 |
| Xia Piep | 610 |
| Xin Ouli | 186 |
## Number of rows changed to Sho Ginq: 6
## Number of rows changed to Kom Oort: 4
**Experiment Day (across all dyads): This variable should count the number of unique experiment days across all dyads. For instance, if multiple dyads have trials on the same date, that date should be considered a single experiment day.
Dyad Day (within each dyad): This variable should count the number of unique experiment days for each dyad separately. The first day of trials for a dyad should be Day 1, the second distinct date of trials should be Day 2, and so on.
DaysSinceStart: Tracks the total number of days since the first experiment, counting every calendar day, including gaps between experiments.
ExperimentDay: Counts unique experiment dates, with each distinct date assigned a consecutive day number.
Trial: Numbers the trials sequentially within each dyad, starting from 1 for each dyad.
## # A tibble: 2,742 × 5
## Date Dyad Trial DyadDay TrialDay
## <date> <chr> <int> <int> <int>
## 1 2022-09-29 Buk Ndaw 1 1 1
## 2 2022-09-29 Buk Ndaw 2 1 2
## 3 2022-09-29 Buk Ndaw 3 1 3
## 4 2022-10-04 Buk Ndaw 4 2 1
## 5 2022-10-13 Buk Ndaw 5 3 1
## 6 2022-10-13 Buk Ndaw 6 3 2
## 7 2022-10-13 Buk Ndaw 7 3 3
## 8 2022-10-13 Buk Ndaw 8 3 4
## 9 2022-10-13 Buk Ndaw 9 3 5
## 10 2022-10-13 Buk Ndaw 10 3 6
## # ℹ 2,732 more rows
## Total Number of Unique Experiment Days: 92
##
## Combined Summary:
| Dyad | Amount of Trials | Number of Days |
|---|---|---|
| Xia Piep | 610 | 49 |
| Sey Sirk | 587 | 53 |
| Kom Oort | 370 | 31 |
| Sho Ginq | 283 | 35 |
| Buk Ndaw | 259 | 39 |
| Pom Xian | 259 | 19 |
| Xin Ouli | 186 | 27 |
| Nge Oerw | 183 | 22 |
| Pom Guat | 5 | 1 |
## Change in Rows: -5
## Placement columns were created successfully.
## Maximum Distance: 10
## Minimum Distance: 0
Reminder: The different behaviors that are coded in DyadResponse are: Distracted, Female aggress male, Male aggress female, Intrusion, Loosing interest, Not approaching, Tolerance and Other
I will create some tables to have a better understanding of the state of the column dyadresponse and the different existing combinations at this point
Also I will create the hierarchy before implementing in the dataset
## Total number of rows in DyadResponse: 2737
## Number of rows with a single entry in DyadResponse: 2466
## Number of rows with multiple entries in DyadResponse: 271
## Number of rows with exactly 2 entries in DyadResponse: 258
## Number of rows with more than 2 entries in DyadResponse: 13
##
##
## Table: Rows with Single Responses in DyadResponse
##
## |DyadResponse | Frequency|
## |:-------------------|---------:|
## |Tolerance | 1809|
## |Not approaching | 465|
## |Male aggress female | 73|
## |Intrusion | 51|
## |Losing interest | 37|
## |Female aggress male | 18|
## |Other | 9|
## |Distracted | 4|
##
##
## Table: Rows with Multiple Responses in DyadResponse
##
## |Combination | Frequency|
## |:-------------------------------------------------|---------:|
## |Intrusion;Tolerance | 49|
## |Losing interest;Not approaching | 47|
## |Intrusion;Not approaching | 30|
## |Male aggress female;Tolerance | 27|
## |Looks at partner;Tolerance | 23|
## |Losing interest;Tolerance | 22|
## |Female aggress male;Tolerance | 19|
## |Looks at partner;Not approaching | 9|
## |Other;Tolerance | 7|
## |Distracted;Not approaching | 6|
## |Not approaching;Tolerance | 4|
## |Distracted;Losing interest | 3|
## |Female aggress male;Male aggress female;Tolerance | 3|
## |Female aggress male;Not approaching | 3|
## |Distracted;Tolerance | 2|
## |Female aggress male;Intrusion;Tolerance | 2|
## |Male aggress female;Not approaching | 2|
## |Distracted;Intrusion;Not approaching | 1|
## |Distracted;Intrusion;Tolerance | 1|
## |Female aggress male;Intrusion | 1|
## |Intrusion;Losing interest | 1|
## |Intrusion;Losing interest;Not approaching | 1|
## |Intrusion;Male aggress female | 1|
## |Intrusion;Not approaching;Tolerance | 1|
## |Looks at partner;Losing interest;Not approaching | 1|
## |Looks at partner;Male aggress female | 1|
## |Looks at partner;Male aggress female;Tolerance | 1|
## |Looks at partner;Other;Tolerance | 1|
## |Losing interest;Not approaching;Tolerance | 1|
## |Not approaching;Other | 1|
## Unique Combinations and Counts for More than 2 Entries:
## Female aggress male & Male aggress female & Tolerance 3
## Female aggress male & Intrusion & Tolerance 2
## Distracted & Intrusion & Not approaching 1
## Distracted & Intrusion & Tolerance 1
## Intrusion & Losing interest & Not approaching 1
## Intrusion & Not approaching & Tolerance 1
## Looks at partner & Losing interest & Not approaching 1
## Looks at partner & Male aggress female & Tolerance 1
## Looks at partner & Other & Tolerance 1
## Losing interest & Not approaching & Tolerance 1
## Rows with multiple entries in DyadResponse: 3 4 13 53 58 64 68 71 72 84 88 89 97 98 113 133 137 192 194 244 275 277 278 284 289 298 302 304 318 322 326 327 328 331 341 346 347 367 368 387 388 391 397 409 454 468 469 482 497 510 511 522 529 536 597 600 626 629 687 696 706 711 713 726 740 748 749 757 760 761 763 768 769 780 786 818 829 841 843 845 861 866 877 888 892 898 912 919 934 937 938 954 955 962 973 997 1004 1015 1027 1042 1048 1060 1104 1141 1143 1164 1165 1190 1209 1213 1217 1225 1229 1236 1240 1243 1245 1247 1254 1267 1271 1284 1288 1293 1308 1310 1311 1326 1327 1330 1336 1342 1344 1353 1378 1398 1399 1407 1416 1461 1464 1473 1478 1484 1488 1489 1495 1510 1520 1521 1528 1558 1595 1598 1602 1607 1616 1629 1630 1633 1636 1639 1642 1645 1654 1658 1679 1708 1714 1718 1724 1741 1750 1752 1784 1785 1788 1795 1800 1801 1802 1804 1805 1806 1807 1808 1815 1816 1819 1820 1821 1822 1824 1825 1828 1829 1830 1831 1832 1837 1843 1856 1857 1858 1889 1897 1968 1984 2023 2066 2088 2090 2091 2100 2101 2102 2103 2104 2105 2109 2112 2113 2147 2148 2149 2152 2153 2166 2175 2184 2185 2192 2193 2194 2199 2201 2202 2205 2218 2226 2237 2248 2257 2290 2294 2343 2353 2354 2359 2360 2362 2363 2370 2388 2395 2396 2400 2465 2466 2467 2493 2534 2601 2609 2626 2629 2630 2643 2651 2657 2730
## Unique Responses and Counts (Sorted Alphabetically):
## Distracted 4
## Distracted Losing interest 3
## Female aggress male 18
## Female aggress male Intrusion 1
## Female aggress male Not approaching 3
## Intrusion 51
## Losing interest 37
## Losing interest Intrusion 1
## Male aggress female 73
## Male aggress female Intrusion 1
## Male aggress female Looks at partner 1
## Male aggress female Not approaching 2
## Not approaching 465
## Not approaching Distracted 6
## Not approaching Distracted Intrusion 1
## Not approaching Intrusion 30
## Not approaching Looks at partner 9
## Not approaching Losing interest 47
## Not approaching Losing interest Intrusion 1
## Not approaching Losing interest Looks at partner 1
## Not approaching Other 1
## Other 9
## Tolerance 1809
## Tolerance Distracted 2
## Tolerance Distracted Intrusion 1
## Tolerance Female aggress male 19
## Tolerance Female aggress male Intrusion 2
## Tolerance Intrusion 49
## Tolerance Looks at partner 23
## Tolerance Looks at partner Other 1
## Tolerance Losing interest 22
## Tolerance Male aggress female 27
## Tolerance Male aggress female Female aggress male 3
## Tolerance Male aggress female Looks at partner 1
## Tolerance Not approaching 4
## Tolerance Not approaching Intrusion 1
## Tolerance Not approaching Losing interest 1
## Tolerance Other 7
## Number of occurrences of 'Looks at partner' in DyadResponse: 36
## Number of rows that have been changed: 36
## Number of 'Other' occurrences in DyadResponse: 18
| Line | MaleID | FemaleID | OtherResponse |
|---|---|---|---|
| 39 | Buk | Ndaw | No Response |
| 286 | Kom | Oort | Ooet scream while at the box and Kom get the corm of both |
| 300 | Kom | Oort | Both at boxes; Kom touching the box of Oort |
| 304 | Kom | Oort | Kom touching Oort’s box. Oort came after 30 sec to her own box. |
| 700 | Nge | Oerw | No Response |
| 815 | Pom | Xian | No Response |
| 826 | Pom | Xian | No Response |
| 1484 | Sey | Sirk | wait for her. moved back and waited for sirk to approach and came back |
| 1659 | Sho | Ginq | No Response |
| 1741 | Sho | Ginq | Sho took one of her corn |
| 2101 | Xia | Piep | aggression |
| 2102 | Xia | Piep | Xia ate corn from both boxes |
| 2103 | Xia | Piep | Xia stolen some pieces from piep box |
| 2104 | Xia | Piep | Xia immediately go for piep box and eat her corn then eat his. |
| 2105 | Xia | Piep | piep got 2 pieces out of 3. Xia took one corn |
| 2170 | Piep | Xia | fem ag male |
| 2184 | Xia | Piep | Neither ID approach |
| 2729 | Xin | Ouli | opened for oerw |
There are 18 cases where was Other as a Response in the experiment. Here are the modifications I will make for them:
-First lets remember that all of these DyadResponse contain Other which will make their identification easier
First I will remove the lines with the combinations Other & No Response
## Number of combinations with 'Other' in DyadResponse and 'No Response' in OtherResponse: 5
## Row numbers: 39, 700, 815, 826, 1659
## Number of rows that have been changed: 5
## Number of occurrences of 'Other' in DyadResponse left: 13
## Number of rows changed: 0
I will create a Hierarchy for DyadResponse in order to treat cases where multiple behaviours were produced within a trial in order to reduce the amount of responses and clear any discrepancies that could be found
Aggression & Tolerance: I will keep Agression as we defined tolerance as the absence of any sign of aggression between individuals of a Dyad. Each Time we wrote tolerance it means that the monkeys touched the boxes at the same time. I will create a variable called SimulatenousTouch to record every time the box did the action at the same time. Then I will in every case of aggression and tolerance within a trial keep Aggression: Aggression > Tolerance. (Note that aggression is not yet a variable in itself but instead we have Male aggress female & Female aggress male as occurences of aggression)
## Number of cases with Tolerance & Male aggress female: 31
## Number of cases with Tolerance & Female aggress male: 24
## Number of rows changed: 55
## Number of lines remaining with both Aggression and Tolerance: 0
## Number of cases with Tolerance & Losing interest: 23
## Number of rows changed: 23
## Number of lines remaining with both Tolerance and Losing interest: 0
## Number of cases with Tolerance & Intrusion: 53
## Number of rows changed: 53
## Number of lines remaining with both Tolerance and Intrusion: 0
## Initial duplicates of 'Not Approaching': 571
## Final duplicates of 'Not Approaching': 634
## Final number of 'Not approaching' occurrences after replacing 'Losing interest': 634
## Number of rows with 'Losing interest' remaining: 0
## Initial duplicates of 'Not Approaching': 634
## Final duplicates of 'Not Approaching': 634
## Final number of 'Not approaching' occurrences after replacing 'Losing interest': 634
## Number of rows with 'Losing interest' remaining: 0
## Initial duplicates of 'Not Approaching; Not Approaching': 0
## Final duplicates of 'Not Approaching; Not Approaching': 0
## Number of cases with Intrusion & Not Approaching: 34
## Number of lines remaining with both Intrusion and Not Approaching: 0
## Initial duplicates of 'Not Approaching; Not Approaching': 0
## Final duplicates of 'Not Approaching; Not Approaching': 0
## Number of changes made to replace 'Distracted' with 'Not approaching': 17
## Number of changes in 'Not approaching' entries: 8
## Number of cases with Male Aggress Female & Not Approaching: 2
## Number of cases with Female Aggress Male & Not Approaching: 3
## Total number of cases with Aggression & Not Approaching: 5
## Number of changes made to remove 'Not Approaching' with 'Male aggress female': 2
## Number of changes made to remove 'Not Approaching' with 'Female aggress male': 3
## Number of remaining cases with Male Aggress Female & Not Approaching: 0
## Number of remaining cases with Female Aggress Male & Not Approaching: 0
## Number of cases with Male Aggress Female & Tolerance: 31
## Number of cases with Female Aggress Male & Tolerance: 24
## Number of changes made to remove 'Male aggress female' with 'Tolerance': 31
## Number of changes made to remove 'Female aggress male' with 'Tolerance': 24
## Number of remaining cases with Male Aggress Female & Tolerance: 0
## Number of remaining cases with Female Aggress Male & Tolerance: 0
## Number of cases with Male aggress female: 77
## Number of cases with Female aggress male: 22
## Total number of aggression cases: 99
## Number of cases with Tolerance & Not Approaching: 30
## Number of rows changed: 30
## Number of lines remaining with both Tolerance and Not Approaching: 0
## Count of 'Not approaching Not approaching': 58
## Count of 'Not approaching Not approaching': 0
## Count of 'Not approaching Not approaching': 0
## Remaining duplicates of 'Not approaching Not approaching': 0
| Keyword | Frequency | |
|---|---|---|
| 10 | Tolerance | 1882 |
| 8 | Not approaching | 572 |
| 4 | Intrusion | 83 |
| 5 | Male aggress female | 74 |
| 12 | Tolerance Intrusion | 49 |
| 13 | Tolerance Not approaching | 29 |
| 1 | Female aggress male | 18 |
| 3 | Female aggress male Not approaching | 3 |
| 11 | Tolerance Intrusion | 3 |
| 7 | Male aggress female Not approaching | 2 |
| 2 | Female aggress male Intrusion | 1 |
| 6 | Male aggress female Intrusion | 1 |
| 9 | Not approaching Intrusion | 1 |
| 14 | Tolerance Not approaching Intrusion | 1 |
## Warning in RColorBrewer::brewer.pal(n, pal): n too large, allowed maximum for palette Set3 is 12
## Returning the palette you asked for with that many colors
| DyadResponse | n |
|---|---|
| Tolerance | 1882 |
| Not approaching | 572 |
| Intrusion | 83 |
| Male aggress female | 74 |
| Tolerance Intrusion | 49 |
| Tolerance Not approaching | 29 |
| Female aggress male | 18 |
Now that all the variables have been treated and cleaned I want to modify thew variables in Bex to keep only the ones I need and to put them in the right order
First lets print the name of all the variables before reordering them and keeping only the variables of interest
## tibble [2,719 × 37] (S3: tbl_df/tbl/data.frame)
## $ Time : chr [1:2719] "08:18:23" "08:37:12" "08:37:58" "08:16:30" ...
## $ Date : chr [1:2719] "2022-09-29" "2022-09-29" "2022-09-29" "2022-10-04" ...
## $ Group : chr [1:2719] "Ankhase" "Ankhase" "Ankhase" "Ankhase" ...
## $ MaleID : chr [1:2719] "Buk" "Buk" "Buk" "Buk" ...
## $ FemaleID : chr [1:2719] "Ndaw" "Ndaw" "Ndaw" "Ndaw" ...
## $ FemaleCorn : chr [1:2719] "0" "0" "0" "3" ...
## $ DyadDistance : chr [1:2719] "5" "5" "5" "5" ...
## $ DyadResponse : chr [1:2719] "Tolerance" "Tolerance" "Intrusion" "Intrusion" ...
## $ OtherResponse : chr [1:2719] "No Response" "No Response" "No Response" "No Response" ...
## $ Audience : chr [1:2719] "Ginq; Gubh" "Ginq; Gubh" "Ginq; Gubh" "Ghid; Gil; Ginq; Gom" ...
## $ IDIndividual1 : chr [1:2719] "No individual" "No individual" "Buk; Ndaw" "Buk" ...
## $ IntruderID : chr [1:2719] "No Intrusion" "No Intrusion" "Ginq; Gubh" "Ginq" ...
## $ Remarks : chr [1:2719] "No Remarks" "No Remarks" "Treated" "Treated" ...
## $ MaleCorn : chr [1:2719] "0" "0" "0" "7" ...
## $ Intrusion : chr [1:2719] "0" "0" "1" "0" ...
## $ AmountAudience : chr [1:2719] "2" "2" "2" "4" ...
## $ Context : chr [1:2719] "NoContext" "NoContext" "NoContext" "NoContext" ...
## $ SpecialBehaviour : chr [1:2719] "NoSpecialBehaviour" "NoSpecialBehaviour" "NoSpecialBehaviour" "NoSpecialBehaviour" ...
## $ GotCorn : chr [1:2719] "Yes" "Yes" "Yes" "Yes" ...
## $ Period : chr [1:2719] "6 to 8" "6 to 8" "6 to 8" "6 to 8" ...
## $ Hour : chr [1:2719] "08:00:00" "08:00:00" "08:00:00" "08:00:00" ...
## $ Day : chr [1:2719] "16" "16" "16" "21" ...
## $ Month : chr [1:2719] "2022-09" "2022-09" "2022-09" "2022-10" ...
## $ Male : chr [1:2719] "Buk" "Buk" "Buk" "Buk" ...
## $ Female : chr [1:2719] "Ndaw" "Ndaw" "Ndaw" "Ndaw" ...
## $ Dyad : chr [1:2719] "Buk Ndaw" "Buk Ndaw" "Buk Ndaw" "Buk Ndaw" ...
## $ DaysSinceStart : chr [1:2719] "16" "16" "16" "21" ...
## $ ExperimentDay : chr [1:2719] "6" "6" "6" "7" ...
## $ ExperimentDay_Verified: chr [1:2719] "6" "6" "6" "7" ...
## $ Trial : chr [1:2719] "1" "2" "3" "4" ...
## $ DyadDay : chr [1:2719] "1" "1" "1" "2" ...
## $ TrialDay : chr [1:2719] "1" "2" "3" "1" ...
## $ PlacementMale : chr [1:2719] "0" "0" "0" "7" ...
## $ PlacementFemale : chr [1:2719] "0" "0" "0" "3" ...
## $ Proximity : chr [1:2719] "4-5" "4-5" "4-5" "4-5" ...
## $ DyadResponse_sorted : chr [1:2719] "Tolerance" "Tolerance" "Intrusion;Not approaching" "Intrusion;Losing interest" ...
## $ MultipleResponses : chr [1:2719] "Single Response" "Single Response" ">1 Response" ">1 Response" ...
## Time Date Group
## 0 0 0
## MaleID FemaleID FemaleCorn
## 0 0 0
## DyadDistance DyadResponse OtherResponse
## 0 0 0
## Audience IDIndividual1 IntruderID
## 0 0 0
## Remarks MaleCorn Intrusion
## 0 0 0
## AmountAudience Context SpecialBehaviour
## 0 0 0
## GotCorn Period Hour
## 0 0 0
## Day Month Male
## 0 0 0
## Female Dyad DaysSinceStart
## 0 0 0
## ExperimentDay ExperimentDay_Verified Trial
## 0 0 0
## DyadDay TrialDay PlacementMale
## 0 0 0
## PlacementFemale Proximity DyadResponse_sorted
## 0 0 0
## MultipleResponses
## 0
Before reordering these variables I may have to re update them to make sure all the changes that occurred between the creation of this variable and now, did not create any errors or mistakes
Finally I will keep the following variables in the next order: Date, Month, Day, Time, Period, Trial, Male, Female, Dyad, DyadDistance, Proximity, DyadResponse, Intrusion, IntruderID SpecialBehaviour, Audience, AmountAudience, Context, Special Behaviour
## [1] "DaysSinceStart" "ExperimentDay" "Date" "Month"
## [5] "DyadDay" "TrialDay" "Trial" "Time"
## [9] "Hour" "Period" "Male" "Female"
## [13] "Dyad" "DyadDistance" "Proximity" "DyadResponse"
## [17] "IDIndividual1" "Intrusion" "IntruderID" "SpecialBehaviour"
## [21] "Audience" "AmountAudience" "Context"
####Dyad, Distance & Date
Trials of grpahs, I will have to check all of them
My goal here is too see if each dyad have an general evolution of their dyad distance trough time and how many varaition do they have
## `geom_smooth()` using formula = 'y ~ x'
## Warning: Removed 65 rows containing non-finite outside the scale range
## (`stat_smooth()`).
## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## `geom_smooth()` using formula = 'y ~ x'
## Warning: There was 1 warning in `mutate()`.
## ℹ In argument: `KinPresent = ifelse(...)`.
## Caused by warning in `grepl()`:
## ! l'argument pattern a une longueur > 1 et seul le premier élément est utilisé
Because I will look at the audience and individual factors on the individuals of the dyads i must make a lsit of all the individudals for which I will need to extract information.
First lets extract all of the names of individuals in audience
Step 1: Clean and Process the Audience Data: We first clean the Audience column, ensuring that any “NA” values are replaced with “No audience” and split the combined names. Whitespace around names is also removed.
Step 2: Recalculate AmountAudience: This step ensures the AmountAudience variable is correctly updated after splitting the Audience column.
Step 3: Summarize Changes: We count how many times the value “NA” was replaced with “No audience” and display this count using cat().
Step 4: Categorize Audience Members: We categorize audience members into “Male”, “Female”, and other relevant categories based on the length of their IDs. Specific identifiers like “unkam”, “unkaf”, etc., are given appropriate labels.
Step 5: Display Tables by Category: We split the audience count table into sections for “Male”, “Female”, and other categories, and then display these tables separately. The use of knitr::kable() provides a professional table format that is more suitable for RMarkdown.
##
## **Table: Top 5 Male Audience Members**
| Audience Member | Count |
|---|---|
| Pix | 177 |
| Oup | 165 |
| Nko | 123 |
| Gom | 103 |
| Ome | 103 |
##
## **Table: Top 5 Female Audience Members**
| Audience Member | Count |
|---|---|
| Sirk | 241 |
| Oort | 199 |
| Oerw | 158 |
| Piep | 154 |
| Ginq | 145 |
##
## **List: All Male Audience Members**
## Pix, Oup, Nko, Gom, Ome, Sey, Nda, Ndl, Gha, Xia, Sig, Xin, Rid, Gib, Gub, Mui, Kom, Dix, Nak, Buk, Aan, Non, Aal, Sho, Guz, Pie, Hee, Pom, Syl, Nuk, App, Gri, Ott, Ree, Gil, Kno, Xar, Goe, Ask, UnkAM, Nuu, Ros, Vla, Xop, Bet, Roc, Hem, Nge, Eis, Umb, Ren, Rim, War, Atj, Gua, Vul, Bob, Nca, Tot, Her, Ram, Dal, Tch, Ris, Vry, Dak, Uls, Flu, Gab, Win
##
## **List: All Female Audience Members**
## Sirk, Oort, Oerw, Piep, Ginq, Obse, Ghid, Reen, Gubh, Ndaw, Sitr, Godu, Skem, Ndum, Naal, Ncok, Miel, Gobe, Nkos, Ouli, Ndon, Ndin, Puol, Giji, Lewe, Enge, Papp, Eina, Hond, Rimp, Heer, Sari, Gree, Olyf, Bela, Aapi, Guba, Popp, Oase, Xeni, Nurk, Haai, Rivi, Gran, Misk, Nooi, Prat, Gris, Pikk, Pann, Prai, Riss, Riva, Ekse, Rede, Griv, Regi, UnkJ, Xati, Rooi, Gaya, Prim, Rafa, Udup, Xala, Palm, Xinp, Guat, Reno, Beir, Gese, Grim, UnkA, Gale, Pret, Prag, Prui, Raba, Rioj, UnkAF, Utic
##
## ### Top 10 Audience Members for Dyad: Buk Ndaw
##
##
## |Audience Member | Count|
## |:---------------|-----:|
## |Ginq | 90|
## |Ghid | 62|
## |No audience | 59|
## |Nda | 47|
## |Ndl | 42|
## |Sho | 40|
## |Gha | 37|
## |Gom | 36|
## |Gubh | 35|
## |Godu | 29|
##
## ### Top 10 Audience Members for Dyad: Kom Oort
##
##
## |Audience Member | Count|
## |:---------------|-----:|
## |No audience | 95|
## |Sirk | 51|
## |Xia | 50|
## |Piep | 47|
## |Reen | 35|
## |Pix | 29|
## |Sey | 29|
## |Ome | 25|
## |Oup | 22|
## |Xin | 21|
##
## ### Top 10 Audience Members for Dyad: Nge Oerw
##
##
## |Audience Member | Count|
## |:---------------|-----:|
## |No audience | 55|
## |Sirk | 54|
## |Oup | 40|
## |Pix | 28|
## |Sey | 19|
## |Sig | 18|
## |Obse | 16|
## |Sitr | 13|
## |Ouli | 12|
## |Aan | 11|
##
## ### Top 10 Audience Members for Dyad: Pom Xian
##
##
## |Audience Member | Count|
## |:---------------|-----:|
## |No audience | 115|
## |Gri | 31|
## |Gree | 26|
## |Xar | 23|
## |Xeni | 21|
## |Gran | 18|
## |Xop | 17|
## |Prat | 16|
## |Gris | 15|
## |Miel | 15|
## |Roc | 15|
##
## ### Top 10 Audience Members for Dyad: Sey Sirk
##
##
## |Audience Member | Count|
## |:---------------|-----:|
## |No audience | 250|
## |Piep | 63|
## |Oort | 51|
## |Oerw | 48|
## |Pix | 48|
## |Oup | 44|
## |Ome | 39|
## |Reen | 37|
## |Obse | 34|
## |Sitr | 29|
##
## ### Top 10 Audience Members for Dyad: Sho Ginq
##
##
## |Audience Member | Count|
## |:---------------|-----:|
## |Ndaw | 90|
## |Gom | 67|
## |Gubh | 61|
## |No audience | 59|
## |Gha | 54|
## |Ghid | 49|
## |Buk | 48|
## |Godu | 46|
## |Gib | 42|
## |Ndum | 37|
##
## ### Top 10 Audience Members for Dyad: Xia Piep
##
##
## |Audience Member | Count|
## |:---------------|-----:|
## |No audience | 178|
## |Sirk | 124|
## |Oort | 90|
## |Pix | 56|
## |Oerw | 51|
## |Xin | 37|
## |Naal | 32|
## |Sitr | 32|
## |Obse | 31|
## |Dix | 29|
##
## ### Top 10 Audience Members for Dyad: Xin Ouli
##
##
## |Audience Member | Count|
## |:---------------|-----:|
## |Oerw | 45|
## |Oort | 39|
## |Sey | 38|
## |No audience | 34|
## |Oup | 30|
## |Piep | 29|
## |Aal | 17|
## |Xia | 17|
## |Obse | 16|
## |Ott | 16|
## |Pix | 16|
## [1] "Ordered Unique Intruder Names:"
## [1] "Buk" "Ghid" "Ginq" "Godu" "Gran" "Gree" "Gri" "Grif" "Gris" "Guat"
## [11] "Gub" "Gubh" "Guz" "Hee" "Kno" "Kom" "Nak" "Nda" "Nge" "Non"
## [21] "Obse" "Oerw" "Oort" "Ouli" "Oup" "Piep" "Pix" "Sey" "Sho" "Sirk"
## [31] "Xia" "Xin" "Xop"
## [1] "Top 10 Most Frequent Intruder Names:"
## Name Count Gender
## 1 Oerw 23 Female
## 2 Sey 19 Male
## 3 Oort 15 Female
## 4 Obse 10 Female
## 5 Ghid 9 Female
## 6 Buk 5 Male
## 7 Ginq 5 Female
## 8 Gubh 5 Female
## 9 Gri 4 Male
## 10 Piep 4 Female
## [1] "Male Names (3 letters):"
## [1] "Guz" "Nda" "Sho" "Sey" "Pix" "Xia" "Oup" "Nak" "Gri" "Xop" "Nge" "Kno"
## [13] "Kom" "Buk" "Gub" "Xin" "Hee" "Non"
## [1] "Female Names (4 letters):"
## [1] "Ginq" "Gubh" "Ghid" "Oerw" "Obse" "Piep" "Ouli" "Gree" "Gran" "Gris"
## [11] "Grif" "Guat" "Sirk" "Oort" "Godu"
## `summarise()` has grouped output by 'Dyad'. You can override using the
## `.groups` argument.
Behaviours tests
## `summarise()` has grouped output by 'Dyad'. You can override using the
## `.groups` argument.
## Code DOB FirstRecorded AdjustedDOB Age
## 1 Ouli <NA> 2010-11-09 2006-11-09 18.029118
## 2 Xian 2012-11-05 2012-11-05 2012-11-05 12.038577
## 3 Piep 2012-01-01 2013-01-18 2012-01-01 12.884590
## 4 Ginq 2014-10-18 2014-10-18 2014-10-18 10.089187
## 5 Oort 2015-11-20 2015-11-20 2015-11-20 8.999500
## 6 Ndaw 2016-02-08 2016-02-08 2016-02-08 8.780468
## 7 Sirk 2017-10-21 2017-10-21 2017-10-21 7.080228
## 8 Oerw 2018-11-09 2018-11-09 2018-11-09 6.028871
Since most males dispsersed from unkonwn groups we dont have their date of birth. But male usually dispers around 5 years old.
Male DOB: Use FirstRecorded - 5 years, except for a few individuals who have DOB recorded (like Xia and Xin)
Use DOB directly if it exists.
Subtract 5 years from FirstRecorded if DOB is missing.
Calculate age based on AdjustedDOB.
Output a table with the relevant columns and create a graph * similar to the one you had for the females.
NOTE FOR UPDATE, the individual with the code Buk, short for BukuBuku died on the 12 of september 2024
## Code DOB FirstRecorded AdjustedDOB Age
## 1 Sey 2014-01-01 2014-12-31 2014-01-01 10.883180
## 2 Xia 2016-11-14 2016-11-14 2016-11-14 8.013854
## 3 Nge 2016-11-18 2016-11-18 2016-11-18 8.002902
## 4 Kom <NA> 2017-09-04 2012-09-04 12.208327
## 5 Pom 2017-10-19 2017-10-19 2017-10-19 7.085703
## 6 Xin 2017-01-01 2017-10-20 2017-01-01 7.882434
## 7 Sho <NA> 2020-10-13 2015-10-13 9.103541
## 8 Buk <NA> 2021-05-11 2016-05-11 8.525842
## Code DOB FirstRecorded AdjustedDOB Age Gender
## 1 Ouli <NA> 2010-11-09 2006-11-09 18.029118 Female
## 2 Xian 2012-11-05 2012-11-05 2012-11-05 12.038577 Female
## 3 Piep 2012-01-01 2013-01-18 2012-01-01 12.884590 Female
## 4 Ginq 2014-10-18 2014-10-18 2014-10-18 10.089187 Female
## 5 Oort 2015-11-20 2015-11-20 2015-11-20 8.999500 Female
## 6 Ndaw 2016-02-08 2016-02-08 2016-02-08 8.780468 Female
## 7 Sirk 2017-10-21 2017-10-21 2017-10-21 7.080228 Female
## 8 Oerw 2018-11-09 2018-11-09 2018-11-09 6.028871 Female
## 9 Sey 2014-01-01 2014-12-31 2014-01-01 10.883180 Male
## 10 Xia 2016-11-14 2016-11-14 2016-11-14 8.013854 Male
## 11 Nge 2016-11-18 2016-11-18 2016-11-18 8.002902 Male
## 12 Kom <NA> 2017-09-04 2012-09-04 12.208327 Male
## 13 Pom 2017-10-19 2017-10-19 2017-10-19 7.085703 Male
## 14 Xin 2017-01-01 2017-10-20 2017-01-01 7.882434 Male
## 15 Sho <NA> 2020-10-13 2015-10-13 9.103541 Male
## 16 Buk <NA> 2021-05-11 2016-05-11 8.525842 Male
We are going to calculate the raw differnce of agin between each male and female in the dyad displaying the absolute age difference ordered
Before that i will create a variable called DyadData with information the each dyad
## # A tibble: 8 × 6
## Dyad Male Female MaleAge FemaleAge AgeDifference
## <chr> <chr> <chr> <dbl> <dbl> <dbl>
## 1 Sey Sirk Sey Sirk 10.9 7.08 3.80
## 2 Xia Piep Xia Piep 8.01 12.9 4.87
## 3 Nge Oerw Nge Oerw 8.00 6.03 1.97
## 4 Sho Ginq Sho Ginq 9.10 10.1 0.986
## 5 Xin Ouli Xin Ouli 7.88 18.0 10.1
## 6 Buk Ndaw Buk Ndaw 8.53 8.78 0.255
## 7 Kom Oort Kom Oort 12.2 9.00 3.21
## 8 Pom Xian Pom Xian 7.09 12.0 4.95
## # A tibble: 8 × 4
## Male Female AgeDifference AgeDirection
## <chr> <chr> <dbl> <chr>
## 1 Sey Sirk 3.80 Male Older
## 2 Xia Piep 4.87 Female Older
## 3 Nge Oerw 1.97 Male Older
## 4 Sho Ginq 0.986 Female Older
## 5 Xin Ouli 10.1 Female Older
## 6 Buk Ndaw 0.255 Female Older
## 7 Kom Oort 3.21 Male Older
## 8 Pom Xian 4.95 Female Older
## # A tibble: 8 × 6
## Dyad TotalTrials ToleranceCount NoToleranceCount ToleranceProportion
## <chr> <int> <dbl> <dbl> <dbl>
## 1 Buk Ndaw 258 145 113 0.562
## 2 Kom Oort 367 305 62 0.831
## 3 Nge Oerw 182 121 61 0.665
## 4 Pom Xian 257 171 86 0.665
## 5 Sey Sirk 586 401 185 0.684
## 6 Sho Ginq 281 158 123 0.562
## 7 Xia Piep 603 473 130 0.784
## 8 Xin Ouli 185 108 77 0.584
## # ℹ 1 more variable: NoToleranceProportion <dbl>
## # A tibble: 8 × 10
## Dyad Male Female AgeDifference AgeDirection TotalTrials ToleranceCount
## <chr> <chr> <chr> <dbl> <chr> <int> <dbl>
## 1 Sey Sirk Sey Sirk 3.80 Male Older 586 401
## 2 Xia Piep Xia Piep 4.87 Female Older 603 473
## 3 Nge Oerw Nge Oerw 1.97 Male Older 182 121
## 4 Sho Ginq Sho Ginq 0.986 Female Older 281 158
## 5 Xin Ouli Xin Ouli 10.1 Female Older 185 108
## 6 Buk Ndaw Buk Ndaw 0.255 Female Older 258 145
## 7 Kom Oort Kom Oort 3.21 Male Older 367 305
## 8 Pom Xian Pom Xian 4.95 Female Older 257 171
## # ℹ 3 more variables: NoToleranceCount <dbl>, ToleranceProportion <dbl>,
## # NoToleranceProportion <dbl>
## `geom_smooth()` using formula = 'y ~ x'
##
## Pearson's product-moment correlation
##
## data: dyad_summary$AgeDifference and dyad_summary$ToleranceProportion
## t = 0.19274, df = 6, p-value = 0.8535
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## -0.6628703 0.7420961
## sample estimates:
## cor
## 0.07844459
##
## Call:
## lm(formula = ToleranceProportion ~ AgeDifference, data = dyad_summary)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.099630 -0.096714 -0.001366 0.041322 0.165240
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.657686 0.062557 10.513 4.35e-05 ***
## AgeDifference 0.002536 0.013155 0.193 0.854
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.1076 on 6 degrees of freedom
## Multiple R-squared: 0.006154, Adjusted R-squared: -0.1595
## F-statistic: 0.03715 on 1 and 6 DF, p-value: 0.8535
We are going to conduct a Spearman’s rank correlation because we don’t have a normal distribution and this methods considers the order of the data
Comparing with Previous Results Earlier Pearson’s Correlation: Correlation Coefficient (r): Approximately 0.0784 P-value: 0.8535 Interpretation: No significant linear relationship. Spearman’s Correlation: Correlation Coefficient (rho): 0.4286 P-value: 0.2992 Interpretation: Suggests a moderate monotonic relationship but not statistically significant. Why the Difference? Pearson’s Correlation: Measures linear relationships. Sensitive to outliers and requires normally distributed variables. Spearman’s Correlation: Measures monotonic relationships. Less sensitive to outliers and does not require normal distribution.
##
## Spearman's rank correlation rho
##
## data: dyad_summary$AgeDifference and dyad_summary$ToleranceProportion
## S = 48, p-value = 0.2992
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
## rho
## 0.4285714
## `geom_smooth()` using formula = 'y ~ x'
Using a Spearman correlation we found a spearmans correlation coefficient (rho) which is equal to 0,4285714 indicating a moderate positive association between ranks of age difference and tolerance proportion. It seems that as the age difference increases, the tolerance proportion increases
But looking at the p-value of 0,2992 we can see that the correlation is not statistically significant and there is a 29,92% probability that the observed correlations occured by chance
X-axis: AgeDifference
Y-axis:
The correlation coefficient (0.0784) is very weak
Also the p.value, 0.8535 is not significant
It seems that they are not significant correlation between age difference and tolerance proportion
## `geom_smooth()` using formula = 'y ~ x'
##
## Welch Two Sample t-test
##
## data: ToleranceProportion by AgeDirection
## t = -1.4071, df = 4.5304, p-value = 0.2242
## alternative hypothesis: true difference in means between group Female Older and group Male Older is not equal to 0
## 95 percent confidence interval:
## -0.27457658 0.08425425
## sample estimates:
## mean in group Female Older mean in group Male Older
## 0.6315716 0.7267327
*p-value of 0.2242 suggests that there is no significant difference in tolerance proportions between dyads where the male is older and those where the female is older at the conventional alpha level of 0.05.
Model 3: Absolute age + Age direction ~ tolerance
First lets do a regression model
##
## Call:
## lm(formula = ToleranceProportion ~ AgeDifference + AgeDirection,
## data = dyad_summary)
##
## Residuals:
## 1 2 3 4 5 6 7 8
## -0.04738 0.14899 -0.05564 -0.04934 -0.08396 -0.04513 0.10302 0.02944
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.605582 0.069737 8.684 0.000335 ***
## AgeDifference 0.006127 0.012566 0.488 0.646498
## AgeDirectionMale Older 0.102800 0.075076 1.369 0.229212
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.1005 on 5 degrees of freedom
## Multiple R-squared: 0.2772, Adjusted R-squared: -0.01193
## F-statistic: 0.9588 on 2 and 5 DF, p-value: 0.4442
## `geom_smooth()` using formula = 'y ~ x'
##
## Call:
## lm(formula = ToleranceProportion ~ AgeDifference * AgeDirection,
## data = dyad_summary)
##
## Residuals:
## 1 2 3 4 5 6 7 8
## -0.06620 0.14940 -0.03185 -0.05146 -0.08012 -0.04772 0.09805 0.02991
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.608345 0.077844 7.815 0.00145 **
## AgeDifference 0.005475 0.014107 0.388 0.71769
## AgeDirectionMale Older 0.030261 0.272127 0.111 0.91681
## AgeDifference:AgeDirectionMale Older 0.023947 0.085541 0.280 0.79340
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.1113 on 4 degrees of freedom
## Multiple R-squared: 0.2911, Adjusted R-squared: -0.2406
## F-statistic: 0.5475 on 3 and 4 DF, p-value: 0.676
## `geom_smooth()` using formula = 'y ~ x'
In order to calculate the Elo Rating, I will have to create different files
As males and females hierarchy are distinct I will have to calculate their hierarchies separately
Speciffically for these codes I may have seqcheck() and elo.seq() functions errors as: -First interaction occurred before presence (approx) meaning that there’s already an interaction recorded before the male was added into that group in the life history
I will now create FinalAgonistic.csv using Agonistic data from 2016 to 2023 and maybe the latest agonistic file until May 2023. Also I Focal data from June 2022 until the latest focals (date to check)
FinalAgonistic.csv will combine the above input files after filtering one-on-one interactions (excluding support interactions). It will involve cleaning and merging the raw data from agonistic interactions and focals to prepare it for the winner-loser calculations.
I will use the information of each first experiment day to know when to stop the elo calculations. For each dyad I will take the beggining of the month of the experiment.
BD1 > SEPT 2022
AK > SEPT 2022
BD2 > DEC 2022
NH > MAR 2023
Also, because the lastest dyad starts it’s first trial after the 10.03.2023 I will remove all the elo data after this date
Note that male and female hierarchy are not comparable on the same scale in vervet monkeys so I will have to conduct 8 Elo Ration calcuations
First we are going to investigate the relationship between distance and tolerance proportion
Summary and Bar chart for Dyad distance and tolerance binomial
## `summarise()` has grouped output by 'DyadDistance'. You can override using the
## `.groups` argument.
##
## Shapiro-Wilk normality test
##
## data: BexClean$DyadDistance
## W = 0.85818, p-value < 2.2e-16
##
## Shapiro-Wilk normality test
##
## data: BexClean$ToleranceBinomial
## W = 0.58055, p-value < 2.2e-16
##
## Wilcoxon rank sum test with continuity correction
##
## data: DyadDistance by ToleranceBinomial
## W = 869192, p-value = 8.901e-06
## alternative hypothesis: true location shift is not equal to 0
##
## Welch Two Sample t-test
##
## data: DyadDistance by ToleranceBinomial
## t = 4.5617, df = 1474.5, p-value = 5.494e-06
## alternative hypothesis: true difference in means between group 0 and group 1 is not equal to 0
## 95 percent confidence interval:
## 0.1868450 0.4687605
## sample estimates:
## mean in group 0 mean in group 1
## 1.860215 1.532412
tep 7: Reporting the Results In your write-up, you would report something like this:
Normality assumption: Both DyadDistance and ToleranceBinomial were found to be non-normally distributed based on the Shapiro-Wilk test (p-values < 0.05). Therefore, we proceeded with non-parametric statistical methods.
Wilcoxon Rank Sum Test: The Wilcoxon test revealed a statistically significant difference in the distance between dyads that exhibited tolerance (p-value = 8.901e-06). This suggests that dyads that were tolerant were, on average, closer than those that were not.
t-test for comparison: A t-test supported the Wilcoxon test’s findings, further demonstrating that mean distances were significantly different (p-value = 5.494e-06).
## `summarise()` has grouped output by 'Dyad', 'DyadDistance'. You can override
## using the `.groups` argument.
## Le chargement a nécessité le package : carData
##
## Attachement du package : 'car'
## L'objet suivant est masqué depuis 'package:dplyr':
##
## recode
##
## Shapiro-Wilk normality test
##
## data: BexClean$DyadDistance
## W = 0.85818, p-value < 2.2e-16
##
## Call:
## glm(formula = ToleranceBinomial ~ DyadDistance * Dyad, family = "binomial",
## data = BexClean)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -2.0219 -1.2745 0.6983 0.9295 1.3036
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 0.499924 0.221339 2.259 0.023906 *
## DyadDistance -0.105085 0.075864 -1.385 0.165997
## DyadKom Oort 1.384136 0.281964 4.909 9.16e-07 ***
## DyadNge Oerw -0.382906 0.321745 -1.190 0.234010
## DyadPom Xian 0.228396 0.302935 0.754 0.450882
## DyadSey Sirk -0.222303 0.263713 -0.843 0.399245
## DyadSho Ginq 0.179562 0.334319 0.537 0.591201
## DyadXia Piep 0.787217 0.250780 3.139 0.001695 **
## DyadXin Ouli -0.005553 0.335118 -0.017 0.986779
## DyadDistance:DyadKom Oort -0.224558 0.124341 -1.806 0.070921 .
## DyadDistance:DyadNge Oerw 0.576587 0.173023 3.332 0.000861 ***
## DyadDistance:DyadPom Xian 0.088974 0.098107 0.907 0.364459
## DyadDistance:DyadSey Sirk 0.430648 0.109089 3.948 7.89e-05 ***
## DyadDistance:DyadSho Ginq -0.056796 0.111836 -0.508 0.611554
## DyadDistance:DyadXia Piep 0.111038 0.115011 0.965 0.334316
## DyadDistance:DyadXin Ouli 0.051274 0.102832 0.499 0.618046
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 3357.2 on 2718 degrees of freedom
## Residual deviance: 3196.1 on 2703 degrees of freedom
## AIC: 3228.1
##
## Number of Fisher Scoring iterations: 4
## `geom_smooth()` using formula = 'y ~ x'
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : pseudoinverse used at -0.025
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : neighborhood radius 1.025
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : reciprocal condition number 8.5375e-30
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : There are other near singularities as well. 1
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : pseudoinverse used at -0.025
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : neighborhood radius 2.025
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : reciprocal condition number 5.5579e-16
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : There are other near singularities as well. 1
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : pseudoinverse used at -0.03
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : neighborhood radius 2.03
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : reciprocal condition number 3.501e-15
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : There are other near singularities as well. 1
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : pseudoinverse used at -0.025
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : neighborhood radius 1.025
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : reciprocal condition number 2.4783e-30
## Warning in simpleLoess(y, x, w, span, degree = degree, parametric =
## parametric, : There are other near singularities as well. 1
## `geom_smooth()` using formula = 'y ~ x'
#6974 Dyad and Tolerance
## Df Sum Sq Mean Sq
## Dyad 7 0.06992 0.009989
##
## Kruskal-Wallis rank sum test
##
## data: ToleranceProportion by Dyad
## Kruskal-Wallis chi-squared = 7, df = 7, p-value = 0.4289
Step-by-Step Approach 1. Visualizing Behavioral Trends Over Time We’ll start by analyzing how key behaviors (like tolerance, male aggression, female aggression, and not approaching) evolve over time.
Goal: Understand if there’s a consistent trend over time, indicating factors like seasonality, repeated interactions, or changes in dyad dynamics. Approach: We’ll create time-series plots to visualize how these behaviors change across trials and experiment days. 2. Analyzing Behavioral Trends by Dyad Each dyad may exhibit different patterns. By separating the trends by dyad, we can:
Goal: Identify whether some dyads are consistently tolerant, aggressive, or show variations depending on conditions like audience or time. Approach: Create faceted plots, where each facet represents a different dyad, showing the evolution of behaviors like distance, tolerance, and aggression. 3. Investigating the Impact of Audience Understanding how audience composition influences behaviors like tolerance or aggression is crucial.
Goal: Determine if audience types (e.g., kin vs. non-kin) have a statistically significant effect on behaviors. Approach: Use boxplots and bar charts to show the relationship between audience composition (kin vs. non-kin) and behaviors. 4. Postulate and Refine Hypotheses After the initial EDA, we’ll refine the hypotheses and decide on the most suitable statistical models.
Examples of Hypotheses: Tolerance is associated with shorter dyad distances over time. Audience presence (kin vs. non-kin) influences the likelihood of aggression. 5. Preparing for Mixed-Effects Models Once the EDA provides insights, we can set up mixed-effects models that account for:
Repeated measures within each dyad. Random effects to capture variability across dyads or trials. 6. Outlier Detection and Analysis Outliers might distort the analysis, so we’ll need to:
Goal: Identify and handle influential outliers that could bias our models. Approach: Use boxplots and statistical diagnostics (e.g., Cook’s distance) to detect outliers.
## Le chargement a nécessité le package : zoo
##
## Attachement du package : 'zoo'
## Les objets suivants sont masqués depuis 'package:base':
##
## as.Date, as.Date.numeric
## Le chargement a nécessité le package : sna
## Le chargement a nécessité le package : statnet.common
##
## Attachement du package : 'statnet.common'
## Les objets suivants sont masqués depuis 'package:base':
##
## attr, order
## Le chargement a nécessité le package : network
##
## 'network' 1.18.2 (2023-12-04), part of the Statnet Project
## * 'news(package="network")' for changes since last version
## * 'citation("network")' for citation information
## * 'https://statnet.org' for help, support, and other information
## sna: Tools for Social Network Analysis
## Version 2.7-2 created on 2023-12-05.
## copyright (c) 2005, Carter T. Butts, University of California-Irvine
## For citation information, type citation("sna").
## Type help(package="sna") to get started.
##
## Attachement du package : 'data.table'
## Les objets suivants sont masqués depuis 'package:lubridate':
##
## hour, isoweek, mday, minute, month, quarter, second, wday, week,
## yday, year
## Les objets suivants sont masqués depuis 'package:dplyr':
##
## between, first, last
## [1] "2020-01-02" "2024-04-15"
## [1] "2020-01-01" "2024-04-20"
## [1] "2020-01-01" "2024-04-20"
## [1] "2020-01-01" "2024-04-20"
## [1] "2022-03-01" "2023-10-01"
## [1] "2021-09-01" "2023-10-01"
## [1] "2021-09-01" "2023-10-01"
## [1] "2021-09-01" "2023-10-01"
## Presence data supplied, see below for details
## Everything seems to be fine with the interaction sequence...OK
##
## #####################################
##
## Presence starts earlier than data...WARNING
## Presence continues beyond data...WARNING
## IDs in datasequence and presence do not match!
## The following IDs occur in the presence data but NOT in the data sequence:...WARNING
## Ati, BBGubh23, BBNdaw22, BBNkos23, BGug19, BNda21, BNya20, Buk, Gha, Ghi, Gib, Giji, Gil, Gobe, Gom, Gon, Gub, Gugu, Gus, Guz, Hlu, Kek, Mat, Nak, Nca, Nci, Nda, Ndi, Ndik, Ndin, Ndl, Ndo, Ndum, Nge, Nko, Nkun, Nyal, Nyan, Sho, Tch, Twe, Vla, Yan
##
## #####################################
## Presence data supplied, see below for details
## Everything seems to be fine with the interaction sequence...OK
##
## #####################################
##
## Presence starts earlier than data...WARNING
## Presence continues beyond data...WARNING
## IDs in datasequence and presence do not match!
## The following IDs occur in the presence data but NOT in the data sequence:...WARNING
## Aal, Aan, Add, Apa, App, Ard, Ask, Atj, Bas, BBDian22, BBNaal23, BBSari23, BBSkem23, BHee19, Bob, BOer21, Boo, BOul21, Bra, BSie16, Dal, Dix, Dok, Dri, Eie, Eis, Ekse, Ela, Fen, Flu, Glo, Goe, Haai, Han, Hee, Heli, Hem, Her, Hia, Hibi, Hipp, Hot, Kno, Kom, Mat, Mimi, Misk, Mui, Naa, Naga, Nak, Nami, Neu, Nge, Non, Noo, Nucl, Nuk, Nul, Nuu, Oase, Olyf, Ome, Ott, Ouma, Oup, Padk, Pal, Papp, Pepe, Pie, Pikk, Pix, PlainJane, Pom, Popp, Potj, Pro, Prui, Pur, Rat, Ree, Reno, Rhe, Rid, Rimp, Rivi, Rooi, Ros, Samp, Sey, Siel, Sig, Sil, Sitr, Sla, Span, Spe, Syl, Ted, Tot, Tow, Ubu, Umb, War, Win, Xia, Xin, Xiu
##
## #####################################
## Presence data supplied, see below for details
## Everything seems to be fine with the interaction sequence...OK
##
## #####################################
##
## Presence starts earlier than data...WARNING
## Presence continues beyond data...WARNING
## IDs in datasequence and presence do not match!
## The following IDs occur in the presence data but NOT in the data sequence:...WARNING
## Bang, BBGran22, BBGran23, BBGree23, BBOort22, BBRegi23, BBXala23, BBXian22, BBXimp23, Bet, BGua20, BGua21, BPre21, BRen20, BRos20, Cai, Can, Cus, Dak, Fle, Gab, Gale, Gan, Gree, Gri, Grim, Griv, Gua, Hav, Kek, Kny, Lif, Lima, Lip, Lug, NewMale3, Oua, Palm, Pom, Prag, Pri, Prim, Pru, Pye, Rafa, Ram, Ren, Renn, Reva, Rey, Rim, Ris, Riv, Riva, Roc, Roes, Roma, Rosl, Sio, Tam, Tch, Udi, Udup, Ula, Uls, Umt, Utic, Utr, Vla, Vry, Vul, War, Was, Xal, Xar, Xati, Xeni, Xia, Xih, Xin, Xop, Yan
##
## #####################################
## [1] "Ndaw" "Ginq" "Ghid" "Gubh" "Godu" "Ncok" "Ndon" "Guba" "Nkos"
## Ginq Ghid Godu Ndaw Ndon Ncok Guba Gubh Nkos
## 0.911 0.644 0.557 0.483 0.310 0.301 0.294 0.286 0.108
## Ginq Ghid Godu Ndaw Ndon Ncok Guba Gubh Nkos
## 0.911 0.644 0.557 0.483 0.310 0.301 0.294 0.286 0.108
## Individual Elo_Score Quartile
## Ginq Ginq 0.911 1st Quartile
## Ghid Ghid 0.644 1st Quartile
## Godu Godu 0.557 2nd Quartile
## Ndaw Ndaw 0.483 2nd Quartile
## Ndon Ndon 0.310 3rd Quartile
## Ncok Ncok 0.301 3rd Quartile
## Guba Guba 0.294 4th Quartile
## Gubh Gubh 0.286 4th Quartile
## Nkos Nkos 0.108 4th Quartile
## Obse Oort Ouli Puol Aapi Sirk Miel Asis Piep Skem Heer Reen Oerw
## 0.917 0.684 0.679 0.576 0.557 0.463 0.463 0.457 0.421 0.417 0.409 0.405 0.404
## Lewe Naal Rede Hond Numb Nooi Gese Sari Riss Enge Pann Nurk Eina
## 0.399 0.386 0.385 0.366 0.355 0.314 0.272 0.240 0.233 0.231 0.207 0.154 0.133
## Obse Ouli Oerw Oort Puol Aapi Piep Lewe Heer Asis Reen Hond Skem
## 0.870 0.720 0.655 0.614 0.584 0.508 0.430 0.428 0.418 0.411 0.403 0.382 0.377
## Sirk Rede Numb Sari Naal Miel Pann Enge Nurk Gese Eina Nooi Riss
## 0.374 0.338 0.306 0.303 0.282 0.246 0.227 0.206 0.183 0.176 0.169 0.115 0.102
## Individual Elo_Score Quartile
## Obse Obse 0.870 1st Quartile
## Ouli Ouli 0.720 1st Quartile
## Oerw Oerw 0.655 1st Quartile
## Oort Oort 0.614 1st Quartile
## Puol Puol 0.584 1st Quartile
## Aapi Aapi 0.508 1st Quartile
## Piep Piep 0.430 1st Quartile
## Lewe Lewe 0.428 2nd Quartile
## Heer Heer 0.418 2nd Quartile
## Asis Asis 0.411 2nd Quartile
## Reen Reen 0.403 2nd Quartile
## Hond Hond 0.382 2nd Quartile
## Skem Skem 0.377 2nd Quartile
## Sirk Sirk 0.374 3rd Quartile
## Rede Rede 0.338 3rd Quartile
## Numb Numb 0.306 3rd Quartile
## Sari Sari 0.303 3rd Quartile
## Naal Naal 0.282 3rd Quartile
## Miel Miel 0.246 3rd Quartile
## Pann Pann 0.227 4th Quartile
## Enge Enge 0.206 4th Quartile
## Nurk Nurk 0.183 4th Quartile
## Gese Gese 0.176 4th Quartile
## Eina Eina 0.169 4th Quartile
## Nooi Nooi 0.115 4th Quartile
## Riss Riss 0.102 4th Quartile
## Gran Guat Prai Upps Gaya Xala Pret Xinp Gris Beir Prat Regi Xian
## 0.893 0.659 0.526 0.474 0.461 0.454 0.449 0.447 0.447 0.447 0.442 0.433 0.416
## Bela Raba Rioj
## 0.285 0.239 0.076
## Gran Guat Xala Gris Gaya Upps Prai Xian Prat Pret Beir Xinp Bela
## 0.929 0.725 0.635 0.622 0.602 0.595 0.570 0.566 0.465 0.463 0.431 0.402 0.348
## Regi Raba Rioj
## 0.321 0.301 0.115
## Individual Elo_Score Quartile
## Gran Gran 0.929 1st Quartile
## Guat Guat 0.725 1st Quartile
## Xala Xala 0.635 1st Quartile
## Gris Gris 0.622 1st Quartile
## Gaya Gaya 0.602 2nd Quartile
## Upps Upps 0.595 2nd Quartile
## Prai Prai 0.570 2nd Quartile
## Xian Xian 0.566 2nd Quartile
## Prat Prat 0.465 3rd Quartile
## Pret Pret 0.463 3rd Quartile
## Beir Beir 0.431 3rd Quartile
## Xinp Xinp 0.402 3rd Quartile
## Bela Bela 0.348 4th Quartile
## Regi Regi 0.321 4th Quartile
## Raba Raba 0.301 4th Quartile
## Rioj Rioj 0.115 4th Quartile
## [1] "Sho" "Vla" "Buk"
## Buk Vla Sho
## 0.934 0.484 0.072
## [1] "2020-01-02" "2024-04-15"
## [1] "2020-01-01" "2024-04-20"
## [1] "2020-01-01" "2024-04-20"
## [1] "2020-01-01" "2024-04-20"
## [1] "2022-01-20" "2023-09-13"
## [1] "2022-01-20" "2023-09-13"
## [1] "2022-01-20" "2023-09-13"
## [1] "2022-01-20" "2023-09-13"
## Presence data supplied, see below for details
## Everything seems to be fine with the interaction sequence...OK
##
## #####################################
##
## Presence starts earlier than data...WARNING
## Presence continues beyond data...WARNING
## IDs in datasequence and presence do not match!
## The following IDs occur in the presence data but NOT in the data sequence:...WARNING
## Ati, BBGubh23, BBNdaw22, BBNkos23, BGug19, BNda21, BNya20, Buk, Gha, Ghi, Gib, Giji, Gil, Gobe, Gom, Gon, Gub, Gugu, Gus, Guz, Hlu, Kek, Mat, Nak, Nca, Nci, Nda, Ndi, Ndik, Ndin, Ndl, Ndo, Ndum, Nge, Nko, Nkun, Nyal, Nyan, Sho, Tch, Twe, Vla, Yan
##
## #####################################
## Presence data supplied, see below for details
## Everything seems to be fine with the interaction sequence...OK
##
## #####################################
##
## IDs in datasequence and presence do not match!
## The following IDs occur in the presence data but NOT in the data sequence:...WARNING
## Aal, Aan, Add, Apa, App, Ard, Ask, Atj, Bas, BBDian22, BBNaal23, BBSari23, BBSkem23, BHee19, Bob, BOer21, Boo, BOul21, Bra, BSie16, Dal, Dix, Dok, Dri, Eie, Eis, Ekse, Ela, Fen, Flu, Glo, Goe, Haai, Han, Hee, Heli, Hem, Her, Hia, Hibi, Hipp, Hot, Kno, Kom, Mat, Mimi, Misk, Mui, Naa, Naga, Nak, Nami, Neu, Nge, Non, Noo, Nucl, Nuk, Nul, Nuu, Oase, Olyf, Ome, Ott, Ouma, Oup, Padk, Pal, Papp, Pepe, Pie, Pikk, Pix, PlainJane, Pom, Popp, Potj, Pro, Prui, Pur, Rat, Ree, Reno, Rhe, Rid, Rimp, Rivi, Rooi, Ros, Samp, Sey, Siel, Sig, Sil, Sitr, Sla, Span, Spe, Syl, Ted, Tot, Tow, Ubu, Umb, War, Win, Xia, Xin, Xiu
## During the following interactions, IDs were absent according to the presence data:...ERROR
## 427, 430, 432, 479, 524, 526
##
## #####################################
## Presence data supplied, see below for details
## Everything seems to be fine with the interaction sequence...OK
##
## #####################################
##
## IDs in datasequence and presence do not match!
## The following IDs occur in the presence data but NOT in the data sequence:...WARNING
## Bang, BBGran22, BBGran23, BBGree23, BBOort22, BBRegi23, BBXala23, BBXian22, BBXimp23, Bet, BGua20, BGua21, BPre21, BRen20, BRos20, Cai, Can, Cus, Dak, Fle, Gab, Gale, Gan, Gree, Gri, Grim, Griv, Gua, Hav, Kek, Kny, Lif, Lima, Lip, Lug, NewMale3, Oua, Palm, Pom, Prag, Pri, Prim, Pru, Pye, Rafa, Ram, Ren, Renn, Reva, Rey, Rim, Ris, Riv, Riva, Roc, Roes, Roma, Rosl, Sio, Tam, Tch, Udi, Udup, Ula, Uls, Umt, Utic, Utr, Vla, Vry, Vul, War, Was, Xal, Xar, Xati, Xeni, Xia, Xih, Xin, Xop, Yan
##
## #####################################
## [1] "Nkos" "Ginq" "Ndaw" "Ghid" "Gubh" "Godu" "Ncok" "Ndon" "Guba"
## Ginq Ghid Godu Ndaw Ndon Ncok Guba Gubh Nkos
## 0.911 0.645 0.553 0.482 0.312 0.299 0.292 0.284 0.109
## Ginq Godu Ghid Gubh Ncok Guba Ndon Ndaw Nkos
## 1.000 0.602 0.545 0.522 0.468 0.468 0.335 0.272 0.000
## Individual Elo_Score Quartile
## Ginq Ginq 1.000 1st Quartile
## Godu Godu 0.602 1st Quartile
## Ghid Ghid 0.545 2nd Quartile
## Gubh Gubh 0.522 2nd Quartile
## Ncok Ncok 0.468 3rd Quartile
## Guba Guba 0.468 3rd Quartile
## Ndon Ndon 0.335 4th Quartile
## Ndaw Ndaw 0.272 4th Quartile
## Nkos Nkos 0.000 4th Quartile
## Obse Oort Ouli Puol Aapi Asis Sirk Piep Miel Skem Heer Reen Oerw
## 0.917 0.707 0.675 0.576 0.537 0.480 0.468 0.434 0.427 0.406 0.402 0.395 0.394
## Lewe Naal Hond Nooi Rede Numb Gese Sari Enge Riss Nurk Pann Eina
## 0.390 0.377 0.359 0.323 0.323 0.311 0.264 0.253 0.232 0.196 0.156 0.136 0.129
## Obse Ouli Oort Puol Aapi Asis Miel Nooi Lewe Skem Naal Oerw Reen
## 1.000 0.782 0.706 0.656 0.537 0.453 0.448 0.435 0.351 0.351 0.351 0.351 0.351
## Heer Sirk Rede Piep Hond Numb Gese Enge Sari Nurk Riss Eina Pann
## 0.341 0.312 0.269 0.261 0.247 0.232 0.209 0.161 0.151 0.083 0.055 0.024 0.000
## Individual Elo_Score Quartile
## Obse Obse 1.000 1st Quartile
## Ouli Ouli 0.782 1st Quartile
## Oort Oort 0.706 1st Quartile
## Puol Puol 0.656 1st Quartile
## Aapi Aapi 0.537 1st Quartile
## Asis Asis 0.453 1st Quartile
## Miel Miel 0.448 1st Quartile
## Nooi Nooi 0.435 2nd Quartile
## Lewe Lewe 0.351 2nd Quartile
## Skem Skem 0.351 2nd Quartile
## Naal Naal 0.351 2nd Quartile
## Oerw Oerw 0.351 2nd Quartile
## Reen Reen 0.351 2nd Quartile
## Heer Heer 0.341 3rd Quartile
## Sirk Sirk 0.312 3rd Quartile
## Rede Rede 0.269 3rd Quartile
## Piep Piep 0.261 3rd Quartile
## Hond Hond 0.247 3rd Quartile
## Numb Numb 0.232 3rd Quartile
## Gese Gese 0.209 4th Quartile
## Enge Enge 0.161 4th Quartile
## Sari Sari 0.151 4th Quartile
## Nurk Nurk 0.083 4th Quartile
## Riss Riss 0.055 4th Quartile
## Eina Eina 0.024 4th Quartile
## Pann Pann 0.000 4th Quartile
## Gran Guat Xala Gris Gaya Upps Xian Prai Pret Prat Beir Xinp Bela
## 0.932 0.724 0.629 0.617 0.600 0.594 0.562 0.557 0.463 0.462 0.428 0.400 0.340
## Regi Raba Rioj
## 0.320 0.297 0.102
## Gran Guat Xala Gris Gaya Upps Prai Xian Prat Pret Beir Xinp Bela
## 0.929 0.725 0.635 0.622 0.602 0.595 0.570 0.566 0.465 0.463 0.431 0.402 0.348
## Regi Raba Rioj
## 0.321 0.301 0.115
## [1] 0
## [1] 0
## 25% 50% 75%
## 0.44450 0.46900 0.50775
## Individual Elo_Score Quartile
## Gran Gran 0.929 1st Quartile
## Guat Guat 0.725 1st Quartile
## Xala Xala 0.635 1st Quartile
## Gris Gris 0.622 1st Quartile
## Gaya Gaya 0.602 2nd Quartile
## Upps Upps 0.595 2nd Quartile
## Prai Prai 0.570 2nd Quartile
## Xian Xian 0.566 2nd Quartile
## Prat Prat 0.465 3rd Quartile
## Pret Pret 0.463 3rd Quartile
## Beir Beir 0.431 3rd Quartile
## Xinp Xinp 0.402 3rd Quartile
## Bela Bela 0.348 4th Quartile
## Regi Regi 0.321 4th Quartile
## Raba Raba 0.301 4th Quartile
## Rioj Rioj 0.115 4th Quartile
## [1] "2022-02-01" "2023-08-12"
## [1] "2022-01-22" "2023-09-11"
## [1] "2022-01-20" "2023-08-16"
## [1] "2022-01-20" "2023-09-13"
## [1] "2022-01-20" "2023-09-13"
## [1] "2022-01-20" "2023-09-13"
## [1] "Date" "Aal" "Aan" "Aapi" "Add" "Apa"
## [7] "App" "Ard" "Asis" "Ask" "Atj" "Bas"
## [13] "BBDian22" "BBNaal23" "BBSari23" "BBSkem23" "BHee19" "Bob"
## [19] "BOer21" "Boo" "BOul21" "Bra" "BSie16" "Dal"
## [25] "Dix" "Dok" "Dri" "Eie" "Eina" "Eis"
## [31] "Ekse" "Ela" "Enge" "Fen" "Flu" "Gese"
## [37] "Glo" "Goe" "Haai" "Han" "Hee" "Heer"
## [43] "Heli" "Hem" "Her" "Hia" "Hibi" "Hipp"
## [49] "Hond" "Hot" "Kno" "Kom" "Lewe" "Mat"
## [55] "Miel" "Mimi" "Misk" "Mui" "Naa" "Naal"
## [61] "Naga" "Nak" "Nami" "Neu" "Nge" "Non"
## [67] "Noo" "Nooi" "Nucl" "Nuk" "Nul" "Numb"
## [73] "Nurk" "Nuu" "Oase" "Obse" "Oerw" "Olyf"
## [79] "Ome" "Oort" "Ott" "Ouli" "Ouma" "Oup"
## [85] "Padk" "Pal" "Pann" "Papp" "Pepe" "Pie"
## [91] "Piep" "Pikk" "Pix" "PlainJane" "Pom" "Popp"
## [97] "Potj" "Pro" "Prui" "Puol" "Pur" "Rat"
## [103] "Rede" "Ree" "Reen" "Reno" "Rhe" "Rid"
## [109] "Rimp" "Riss" "Rivi" "Rooi" "Ros" "Samp"
## [115] "Sari" "Sey" "Siel" "Sig" "Sil" "Sirk"
## [121] "Sitr" "Skem" "Sla" "Span" "Spe" "Syl"
## [127] "Ted" "Tot" "Tow" "Ubu" "Umb" "War"
## [133] "Win" "Xia" "Xin" "Xiu"
## [1] "Date" "Bang" "BBGran22" "BBGran23" "BBGree23" "BBOort22"
## [7] "BBRegi23" "BBXala23" "BBXian22" "BBXimp23" "Beir" "Bela"
## [13] "Bet" "BGua20" "BGua21" "BPre21" "BRen20" "BRos20"
## [19] "Cai" "Can" "Cus" "Dak" "Fle" "Gab"
## [25] "Gale" "Gan" "Gaya" "Gran" "Gree" "Gri"
## [31] "Grim" "Gris" "Griv" "Gua" "Guat" "Hav"
## [37] "Kek" "Kny" "Lif" "Lima" "Lip" "Lug"
## [43] "NewMale3" "Oua" "Palm" "Pom" "Prag" "Prai"
## [49] "Prat" "Pret" "Pri" "Prim" "Pru" "Pye"
## [55] "Raba" "Rafa" "Ram" "Regi" "Ren" "Renn"
## [61] "Reva" "Rey" "Rim" "Rioj" "Ris" "Riv"
## [67] "Riva" "Roc" "Roes" "Roma" "Rosl" "Sio"
## [73] "Tam" "Tch" "Udi" "Udup" "Ula" "Uls"
## [79] "Umt" "Upps" "Utic" "Utr" "Vla" "Vry"
## [85] "Vul" "War" "Was" "Xal" "Xala" "Xar"
## [91] "Xati" "Xeni" "Xia" "Xian" "Xih" "Xin"
## [97] "Xinp" "Xop" "Yan"
## [1] "Vla" "Buk" "Sho"
## [1] "Sho" "Buk" "Vla"
## [1] "Sey" "Pal" "Umb" "Dok" "Pro" "Kom" "Flu" "Xia" "Nge" "Nak" "Pom" "War"
## [13] "Hee" "Nuu" "Xin" "Kno" "Bra" "Win" "Dal" "Boo"
## [1] "Ted" "Kom" "Nge" "Xia" "Flu" "Umb" "Nak" "Dok" "War" "Sey" "Xin" "Apa"
## [13] "Pom" "Aan" "Kno" "Bob" "Dix" "Bra" "Tot" "Dal" "Win" "Dri"
## [1] "Lif" "Fle" "Tam" "Vul" "Tch" "Utr" "Vry" "Dak" "Pom" "Hav" "Kny" "Cai"
## [1] "Vul" "Fle" "Utr" "Lif" "Tch" "Tam" "Dak" "Vry" "Pom" "Kny" "Hav" "Cai"
## Presence data supplied, see below for details
## Everything seems to be fine with the interaction sequence...OK
##
## #####################################
##
## Presence starts earlier than data...WARNING
## Presence continues beyond data...WARNING
## IDs in datasequence and presence do not match!
## The following IDs occur in the presence data but NOT in the data sequence:...WARNING
## Ati, BBGubh23, BBNdaw22, BBNkos23, BGug19, BNda21, BNya20, Gha, Ghi, Ghid, Gib, Giji, Gil, Ginq, Gobe, Godu, Gom, Gon, Gub, Guba, Gubh, Gugu, Gus, Guz, Hlu, Kek, Mat, Nak, Nca, Nci, Ncok, Nda, Ndaw, Ndi, Ndik, Ndin, Ndl, Ndo, Ndon, Ndum, Nge, Nko, Nkos, Nkun, Nyal, Nyan, Tch, Twe, Yan
##
## #####################################
## Presence data supplied, see below for details
## The following individuals were observed only on one day: Nuu ...WARNING
##
## #####################################
##
## Presence starts earlier than data...WARNING
## Presence continues beyond data...WARNING
## IDs in datasequence and presence do not match!
## The following IDs occur in the presence data but NOT in the data sequence:...WARNING
## Aal, Aapi, Add, App, Ard, Asis, Ask, Atj, Bas, BBDian22, BBNaal23, BBSari23, BBSkem23, BHee19, BOer21, BOul21, Bra, BSie16, Dok, Eie, Eina, Eis, Ekse, Ela, Enge, Fen, Gese, Glo, Goe, Haai, Han, Heer, Heli, Hem, Her, Hia, Hibi, Hipp, Hond, Hot, Lewe, Mat, Miel, Mimi, Misk, Mui, Naa, Naal, Naga, Nami, Neu, Non, Noo, Nooi, Nucl, Nuk, Nul, Numb, Nurk, Oase, Obse, Oerw, Olyf, Ome, Oort, Ott, Ouli, Ouma, Oup, Padk, Pal, Pann, Papp, Pepe, Pie, Piep, Pikk, Pix, PlainJane, Popp, Potj, Pro, Prui, Puol, Pur, Rat, Rede, Ree, Reen, Reno, Rhe, Rid, Rimp, Riss, Rivi, Rooi, Ros, Samp, Sari, Siel, Sig, Sil, Sirk, Sitr, Skem, Sla, Span, Spe, Syl, Ted, Tot, Tow, Ubu, Win, Xiu
##
## #####################################
## Presence data supplied, see below for details
## Everything seems to be fine with the interaction sequence...OK
##
## #####################################
##
## Presence continues beyond data...WARNING
## IDs in datasequence and presence do not match!
## The following IDs occur in the presence data but NOT in the data sequence:...WARNING
## Bang, BBGran22, BBGran23, BBGree23, BBOort22, BBRegi23, BBXala23, BBXian22, BBXimp23, Beir, Bela, Bet, BGua20, BGua21, BPre21, BRen20, BRos20, Can, Cus, Gab, Gale, Gan, Gaya, Gran, Gree, Gri, Grim, Gris, Griv, Gua, Guat, Kek, Lima, Lip, Lug, NewMale3, Oua, Palm, Prag, Prai, Prat, Pret, Pri, Prim, Pru, Pye, Raba, Rafa, Ram, Regi, Ren, Renn, Reva, Rey, Rim, Rioj, Ris, Riv, Riva, Roc, Roes, Roma, Rosl, Sio, Udi, Udup, Ula, Uls, Umt, Upps, Utic, Vla, War, Was, Xal, Xala, Xar, Xati, Xeni, Xia, Xian, Xih, Xin, Xinp, Xop, Yan
##
## #####################################
## Vla Buk Sho
## 1174 1036 790
## Buk Vla Sho
## 1237 995 768
## Buk Vla Sho
## 1242 995 763
## Buk Vla Sho
## 1.000 0.484 0.000
## Buk Vla Sho
## 1.000 0.484 0.000
## Individual Elo_Score Quartile
## Buk Buk 1.000 1st Quartile
## Vla Vla 0.484 3rd Quartile
## Sho Sho 0.000 4th Quartile
## Sey Flu Umb Xia Kom Pom Hee Nuu Aan Apa Xin Nge Nak War
## 1227 1186 1166 1076 1073 1000 1000 1000 1000 919 899 880 832 742
## Sey Flu Umb Kom Xia Pom Hee Nuu Aan Xin Nge Nak War
## 1227 1221 1166 1073 1064 1000 1000 1000 1000 899 880 809 742
## Sey Flu Umb Kom Xia Pom Hee Nuu Aan Xin Nge Nak War
## 1300 1176 1166 1073 1036 1000 1000 1000 1000 899 880 809 742
## Sey Flu Umb Kom Xia Pom Hee Nuu Aan Xin Nge Nak War
## 1.000 0.988 0.874 0.682 0.664 0.532 0.532 0.532 0.532 0.324 0.285 0.138 0.000
## Sey Flu Umb Kom Xia Pom Hee Nuu Aan Xin Nge Nak War
## 1.000 0.988 0.874 0.682 0.664 0.532 0.532 0.532 0.532 0.324 0.285 0.138 0.000
## Individual Elo_Score Quartile
## Sey Sey 1.000 1st Quartile
## Flu Flu 0.988 1st Quartile
## Umb Umb 0.874 1st Quartile
## Kom Kom 0.682 2nd Quartile
## Xia Xia 0.664 2nd Quartile
## Pom Pom 0.532 3rd Quartile
## Hee Hee 0.532 3rd Quartile
## Nuu Nuu 0.532 3rd Quartile
## Aan Aan 0.532 3rd Quartile
## Xin Xin 0.324 4th Quartile
## Nge Nge 0.285 4th Quartile
## Nak Nak 0.138 4th Quartile
## War War 0.000 4th Quartile
## Lif Vul Tch Utr Tam Fle Vry
## 1340 1041 982 936 856 845 NA
## Vul Tch Utr Tam Fle Vry Dak
## 1041 982 936 856 845 NA NA
## Vul Vry Utr Tam Tch Fle Dak
## 1041 1029 936 923 915 816 NA
## Vul Tch Utr Tam Fle Vry Dak
## 1.000 0.699 0.464 0.056 0.000 NA NA
## Vul Tch Utr Tam Fle Vry Dak
## 1.000 0.699 0.464 0.056 0.000 NA NA
## 25% 50% 75%
## 0.056 0.464 0.699
## Individual Elo_Score Quartile
## Vul Vul 1.000 1st Quartile
## Tch Tch 0.699 2nd Quartile
## Utr Utr 0.464 3rd Quartile
## Tam Tam 0.056 4th Quartile
## Fle Fle 0.000 4th Quartile
## Individual Elo_Score Quartile
## Vul Vul 1.000 1st Quartile
## Tch Tch 0.699 2nd Quartile
## Utr Utr 0.464 3rd Quartile
## Tam Tam 0.056 4th Quartile
## Fle Fle 0.000 4th Quartile
## Buk Sho Vla
## 0.616 0.428 0.392
## Flu Kom Sey Umb Hee Boo Nak Pom Nuu Xin Dri Dal Apa
## 0.808 0.803 0.764 0.630 0.596 0.555 0.524 0.518 0.511 0.465 0.452 0.380 0.375
## Bob Nge Aan Xia War Dix Kno
## 0.340 0.314 0.283 0.258 0.225 0.224 0.183
## Lif Cai Vul Hav Utr Vry Tch Dak Kny Tam Fle Pom
## 0.911 0.855 0.587 0.527 0.439 0.409 0.362 0.327 0.320 0.318 0.225 0.219
## [1] "AF"
## [2] "after_rows"
## [3] "AK_quartiles"
## [4] "AKELOF"
## [5] "AKELOF2"
## [6] "AKELOM"
## [7] "AKELOM2"
## [8] "AKF"
## [9] "AKFQ"
## [10] "AKFQ2"
## [11] "AKM"
## [12] "AKMQ"
## [13] "AKMQ2"
## [14] "AKpres"
## [15] "all_combinations_more_than_2"
## [16] "all_correct"
## [17] "AM"
## [18] "anova_model"
## [19] "audience_aggression"
## [20] "audience_count"
## [21] "backup_Bex"
## [22] "BackupbeforeDistanceNA"
## [23] "BD_quartiles"
## [24] "BDELOF"
## [25] "BDELOF2"
## [26] "BDELOM"
## [27] "BDELOM2"
## [28] "BDF"
## [29] "BDFQ"
## [30] "BDFQ2"
## [31] "BDM"
## [32] "BDMQ"
## [33] "BDMQ2"
## [34] "BDpres"
## [35] "before_rows"
## [36] "beforesummary"
## [37] "beforesummaryrowswithmultipleentries"
## [38] "Bex"
## [39] "Bex_beforeNA"
## [40] "bex_dyads"
## [41] "BexClean"
## [42] "BexClean_cleaned"
## [43] "BexClean_expanded"
## [44] "bexclean_na_count"
## [45] "BexClean_updated"
## [46] "BexClean_updated_aggression"
## [47] "BexClean_updated_intrusion"
## [48] "BexClean_updated_losing_interest"
## [49] "BexClean_updated_tolerance_not_approaching"
## [50] "both_missing_after"
## [51] "both_missing_count"
## [52] "both_missing_rows"
## [53] "Boxex"
## [54] "categorize_individual"
## [55] "changed_rows_count"
## [56] "changes_made"
## [57] "changes_made_audience"
## [58] "changes_made_female"
## [59] "changes_made_male"
## [60] "changes_made_other_response"
## [61] "checkbehaviour"
## [62] "checkpresence"
## [63] "combinations_summary"
## [64] "combined_context"
## [65] "combined_data"
## [66] "combined_summary"
## [67] "common_na_count"
## [68] "context_na_count"
## [69] "cor_test"
## [70] "correct_dyads"
## [71] "count"
## [72] "count_actual_remarks"
## [73] "count_aggression_not_approaching"
## [74] "count_aggression_tolerance"
## [75] "count_combinations"
## [76] "count_female_aggress_male"
## [77] "count_female_aggress_male_not_approaching"
## [78] "count_female_aggress_male_tolerance"
## [79] "count_male_aggress_female"
## [80] "count_male_aggress_female_not_approaching"
## [81] "count_male_aggress_female_tolerance"
## [82] "count_no_remarks"
## [83] "count_not_approaching_duplicates"
## [84] "count_tolerance_female_aggression"
## [85] "count_tolerance_intrusion"
## [86] "count_tolerance_losing_interest"
## [87] "count_tolerance_male_aggression"
## [88] "count_tolerance_not_approaching"
## [89] "counts_aggression_not_approaching"
## [90] "counts_aggression_tolerance"
## [91] "create_aggression_variable"
## [92] "createnullot"
## [93] "createnullpresence"
## [94] "custom_colors"
## [95] "d"
## [96] "daily_plot"
## [97] "df_unique_responses"
## [98] "distance_freq"
## [99] "distance_freq_table"
## [100] "dot"
## [101] "DSI"
## [102] "dyad"
## [103] "dyad_counts"
## [104] "dyad_data_clean"
## [105] "dyad_distance_missing_data"
## [106] "dyad_distance_na_row"
## [107] "dyad_kin_analysis"
## [108] "dyad_plot"
## [109] "dyad_response_plot"
## [110] "dyad_summary"
## [111] "dyad_table"
## [112] "dyad_with_age"
## [113] "dyad_with_male_age"
## [114] "dyads"
## [115] "entry"
## [116] "extract_unique_keywords"
## [117] "faceted_monthly_plot_with_lines"
## [118] "female_list"
## [119] "female_missing_after"
## [120] "female_missing_rows"
## [121] "female_names"
## [122] "females"
## [123] "filtered_data"
## [124] "filtered_dyad_response_table"
## [125] "filtered_dyads"
## [126] "final_check"
## [127] "find_and_remove_other_no_response"
## [128] "gender_plot"
## [129] "gender_presence_plot"
## [130] "gender_summary"
## [131] "generate_dyad_tables"
## [132] "got_corn_na_count"
## [133] "grand_means"
## [134] "handle_aggression_tolerance"
## [135] "handle_intrusion_not_approaching"
## [136] "handle_tolerance_intrusion"
## [137] "handle_tolerance_losing_interest"
## [138] "handle_tolerance_not_approaching"
## [139] "initial_count"
## [140] "initial_na_audience"
## [141] "initial_na_count"
## [142] "initial_na_other_response"
## [143] "intruder_data"
## [144] "intruder_name_freq"
## [145] "intruder_names"
## [146] "intrusion_counts"
## [147] "intrusion_plot"
## [148] "intrusion_summary"
## [149] "keyword_counts"
## [150] "keyword_plot"
## [151] "keywords"
## [152] "kin_aggression"
## [153] "kruskal_test"
## [154] "kw"
## [155] "LHdata"
## [156] "LHdata_females"
## [157] "LHdata_females_table"
## [158] "LHdata_males"
## [159] "LHdata_males_table"
## [160] "logistic_model"
## [161] "looks_at_partner_count"
## [162] "makedyads"
## [163] "male_list"
## [164] "male_missing_after"
## [165] "male_missing_rows"
## [166] "male_names"
## [167] "males"
## [168] "max_distance"
## [169] "mean_age"
## [170] "mean_agg"
## [171] "mean_days"
## [172] "mean_female_corn"
## [173] "mean_male_corn"
## [174] "mean_trials"
## [175] "mergeactrec"
## [176] "min_distance"
## [177] "missingdata"
## [178] "model"
## [179] "model3"
## [180] "model3_interaction"
## [181] "monthly_aggression"
## [182] "monthly_stats"
## [183] "multiple_responses"
## [184] "na_count_female"
## [185] "na_count_male"
## [186] "na_counts"
## [187] "NA_DyadDistance"
## [188] "na_dyadresponse_rows"
## [189] "na_row_index"
## [190] "na_summary"
## [191] "NH_quartiles"
## [192] "NHELOF"
## [193] "NHELOF2"
## [194] "NHELOM"
## [195] "NHELOM2"
## [196] "NHF"
## [197] "NHFQ"
## [198] "NHFQ2"
## [199] "NHM"
## [200] "NHMQ"
## [201] "NHMQ2"
## [202] "NHpres"
## [203] "no_intrusion_count"
## [204] "no_remarks_count"
## [205] "not_approaching_counts"
## [206] "not_approaching_data"
## [207] "num_multiple_responses"
## [208] "num_single_responses"
## [209] "only_female_missing_count"
## [210] "only_female_missing_rows"
## [211] "other_occurrences"
## [212] "other_occurrences_df"
## [213] "overall_aggression"
## [214] "overall_aggression_plot"
## [215] "plot"
## [216] "plot_days"
## [217] "plot_distance_freq"
## [218] "plot_female_corn"
## [219] "plot_male_corn"
## [220] "plot_proximity_freq"
## [221] "plot_trials"
## [222] "proximity_colors"
## [223] "proximity_freq"
## [224] "remaining_aggression_tolerance"
## [225] "remaining_count"
## [226] "remaining_female_aggress_male_not_approaching"
## [227] "remaining_female_aggress_male_tolerance"
## [228] "remaining_male_aggress_female_not_approaching"
## [229] "remaining_male_aggress_female_tolerance"
## [230] "remaining_na"
## [231] "remaining_na_audience"
## [232] "remaining_na_count"
## [233] "remaining_na_other_response"
## [234] "remaining_other_count"
## [235] "remaining_tolerance_intrusion"
## [236] "remaining_tolerance_losing_interest"
## [237] "remaining_tolerance_not_approaching"
## [238] "remove_aggression_with_tolerance"
## [239] "remove_duplicate_not_approaching"
## [240] "remove_not_approaching_with_aggression"
## [241] "replace_distracted_with_not_approaching"
## [242] "replace_losing_interest"
## [243] "replacement_values"
## [244] "result"
## [245] "result_aggression"
## [246] "result_remove_aggression"
## [247] "result_remove_not_approaching"
## [248] "result_tolerance_intrusion"
## [249] "result_tolerance_losing_interest"
## [250] "result_tolerance_not_approaching"
## [251] "results"
## [252] "row_num"
## [253] "row_numbers"
## [254] "rows_changed"
## [255] "rows_changed_intrusion"
## [256] "rows_changed_losing_interest"
## [257] "rows_changed_tolerance_not_approaching"
## [258] "rows_to_correct_kom"
## [259] "rows_to_correct_sho"
## [260] "rows_to_delete"
## [261] "rows_with_more_than_two_entries"
## [262] "rows_with_two_entries"
## [263] "simplified_summary_table"
## [264] "single_responses"
## [265] "sort_responses"
## [266] "spearman_test"
## [267] "special_behaviour_na_count"
## [268] "summary_Day"
## [269] "summary_female_corn"
## [270] "summary_male_corn"
## [271] "summary_plot"
## [272] "summary_stats"
## [273] "summary_trial"
## [274] "t_test"
## [275] "t_test_result"
## [276] "table_remarks"
## [277] "tolerance_boxplot"
## [278] "tolerance_counts_long"
## [279] "tolerance_proportion_summary"
## [280] "tolerance_proportions_long"
## [281] "tolerance_summary"
## [282] "tolerance_summary_per_dyad"
## [283] "tolerance_summary_table"
## [284] "top_10_audience"
## [285] "top_10_plot"
## [286] "top_female_count"
## [287] "top_male_count"
## [288] "total_aggression_count"
## [289] "total_aggression_not_approaching"
## [290] "total_aggression_tolerance"
## [291] "total_keyword_occurrences"
## [292] "total_rows"
## [293] "unique_combinations"
## [294] "unique_dyad_response_table"
## [295] "unique_dyads"
## [296] "unique_female_ids"
## [297] "unique_keywords_df"
## [298] "unique_male_ids"
## [299] "unique_responses_summary"
## [300] "unwanted_ids_akm"
## [301] "unwanted_ids_bdm"
## [302] "unwanted_ids_nhm"
## [303] "update_rows"
## [304] "valid_dyads"
## [305] "varname"
## [306] "wilcox_test"
## [307] "wrong_rows"
## [308] "wrong_rows_after_correction"
## [309] "zero_count"
## [1] "2022-01-03" "2024-05-27"
## [1] "2022-06-02" "2024-05-27"
## [1] "2022-01-03" "2023-10-19"
## [1] "2022-01-03" "2023-10-18"
## [1] "2022-01-03" "2023-10-17"
## [1] "2022-06-04" "2023-11-18"
## [1] "2022-01-03" "2023-11-18"
## [1] "date" "focal" "actor" "receiver" "beh" "dur"
## [1] "focal" "date" "OT"
## Date role ID INDEX
## [1,] "2022-08-23" "focal" "Godu" "191"
## [2,] "2022-08-23" "actor" "Godu" "191"
## [3,] "2022-08-23" "receiver" "Gobe" "191"
## [4,] "2022-08-23" "focal" "Godu" "192"
## [5,] "2022-08-23" "actor" "Gobe" "192"
## [6,] "2022-08-23" "receiver" "Godu" "192"
## [7,] "2022-08-23" "focal" "Godu" "193"
## [8,] "2022-08-23" "actor" "Godu" "193"
## [9,] "2022-08-23" "receiver" "Guba" "193"
## [10,] "2022-08-23" "focal" "Godu" "194"
## [11,] "2022-08-23" "actor" "Gobe" "194"
## [12,] "2022-08-23" "receiver" "Godu" "194"
## [13,] "2022-08-23" "focal" "Godu" "195"
## [14,] "2022-08-23" "actor" "Godu" "195"
## [15,] "2022-08-23" "receiver" "Gobe" "195"
## [16,] "2022-08-23" "focal" "Godu" "196"
## [17,] "2022-08-23" "actor" "Godu" "196"
## [18,] "2022-08-23" "receiver" "Gobe" "196"
## [19,] "2023-03-16" "actor" "Nkos" "726"
## [20,] "2023-03-16" "receiver" "Nkos" "727"
## [21,] "2023-03-16" "actor" "Nkos" "728"
## [22,] "2023-03-15" "receiver" "Ndin" "734"
## [23,] "2023-03-15" "receiver" "Ndin" "735"
## [24,] "2023-03-15" "receiver" "Ndin" "736"
## [25,] "2023-03-15" "receiver" "Ndin" "737"
## [26,] "2023-03-15" "receiver" "Ndin" "738"
## [27,] "2023-03-15" "receiver" "Ndin" "739"
## [28,] "2023-03-13" "receiver" "Ndin" "767"
## [29,] "2023-02-18" "receiver" "Nca" "846"
## [30,] "2023-02-18" "receiver" "Nca" "847"
## [31,] "2023-02-18" "receiver" "Nca" "848"
## [32,] "2023-01-25" "receiver" "Nca" "876"
## [33,] "2023-01-25" "receiver" "Nca" "877"
## [34,] "2023-01-25" "receiver" "Nca" "881"
## [35,] "2023-01-25" "receiver" "Nca" "883"
## [36,] "2023-01-25" "receiver" "Nca" "885"
## [37,] "2023-01-21" "receiver" "Nca" "894"
## [38,] "2023-01-21" "receiver" "Nca" "895"
## [39,] "2023-01-19" "receiver" "Ndin" "904"
## [40,] "2023-01-19" "receiver" "Ndin" "905"
## [41,] "2023-01-07" "receiver" "Ndin" "912"
## [42,] "2022-12-20" "receiver" "Ndin" "926"
## [43,] "2022-12-16" "receiver" "Nca" "954"
## [44,] "2022-12-16" "actor" "Nca" "955"
## [45,] "2022-12-16" "receiver" "Nca" "956"
## [46,] "2022-12-10" "receiver" "Nca" "994"
## [47,] "2022-11-19" "receiver" "Nko" "1243"
## [48,] "2022-11-19" "receiver" "Nko" "1244"
## [49,] "2022-11-19" "receiver" "Nko" "1245"
## [50,] "2022-11-19" "actor" "Ncok" "1246"
## [51,] "2022-11-19" "receiver" "Ncok" "1247"
## [52,] "2023-05-03" "actor" "Gom" "1383"
## [53,] "2023-05-03" "receiver" "Gom" "1384"
## [54,] "2023-05-03" "receiver" "Gom" "1385"
## [55,] "2023-05-03" "actor" "Gom" "1386"
## [56,] "2023-06-02" "focal" "Nda" "1485"
## [57,] "2023-06-02" "actor" "Nda" "1485"
## [58,] "2023-06-02" "focal" "Nda" "1486"
## [59,] "2023-06-02" "actor" "Nda" "1486"
## [60,] "2023-06-02" "focal" "Nda" "1487"
## [61,] "2023-06-02" "actor" "Nda" "1487"
## [62,] "2023-06-02" "focal" "Nda" "1488"
## [63,] "2023-06-02" "actor" "Nda" "1488"
## [64,] "2023-06-02" "focal" "Nda" "1489"
## [65,] "2023-06-02" "actor" "Nda" "1489"
## [66,] "2023-06-02" "focal" "Nda" "1490"
## [67,] "2023-06-02" "actor" "Nda" "1490"
## [68,] "2023-06-02" "focal" "Sho" "1491"
## [69,] "2023-06-02" "actor" "Sho" "1491"
## [70,] "2023-06-02" "receiver" "Ndaw" "1491"
## [71,] "2023-06-02" "focal" "Sho" "1492"
## [72,] "2023-06-02" "actor" "Ndaw" "1492"
## [73,] "2023-06-02" "receiver" "Sho" "1492"
## [74,] "2023-06-14" "focal" "Sho" "1524"
## [75,] "2023-06-14" "actor" "Sho" "1524"
## [76,] "2023-07-19" "actor" "Ndum" "1644"
## [77,] "2023-07-19" "receiver" "Ndaw" "1645"
## [78,] "2023-07-19" "actor" "Ndaw" "1646"
## [79,] "2023-07-19" "focal" "Gubh" "1648"
## [80,] "2023-07-19" "actor" "Gubh" "1648"
## [81,] "2023-09-27" "focal" "Buk" "1877"
## [82,] "2023-09-27" "receiver" "Buk" "1877"
## [83,] "2023-09-27" "focal" "Buk" "1878"
## [84,] "2023-09-27" "actor" "Buk" "1878"
## [85,] "2023-09-27" "focal" "Buk" "1879"
## [86,] "2023-09-27" "actor" "Buk" "1879"
## [87,] "2023-08-11" "focal" "Nda" "8093"
## [88,] "2023-08-11" "actor" "Nda" "8093"
## [89,] "2023-07-19" "focal" "Ndaw" "8097"
## [90,] "2023-07-19" "actor" "Ndaw" "8097"
## [91,] "2023-07-19" "actor" "Ndum" "8098"
## [92,] "2022-04-01" "receiver" "Ndl" "8344"
## [93,] "2022-04-01" "receiver" "Ndl" "8345"
## [94,] "2022-08-19" "receiver" "Guba" "8478"
## [95,] "2022-08-23" "focal" "Godu" "8480"
## [96,] "2022-08-23" "actor" "Godu" "8480"
## [97,] "2022-08-23" "receiver" "Gobe" "8480"
## [98,] "2023-03-13" "receiver" "Ndin" "8652"
## [99,] "2023-03-13" "receiver" "Ndin" "8656"
## [100,] "2023-03-14" "receiver" "Ndin" "8658"
## [101,] "2023-03-16" "receiver" "Nca" "8661"
## [102,] "2023-03-16" "receiver" "Nca" "8662"
## [103,] "2023-03-16" "receiver" "Nca" "8664"
## [104,] "2023-03-16" "receiver" "Nca" "8666"
## [105,] "2023-03-16" "receiver" "Nca" "8667"
## [106,] "2023-08-11" "focal" "Buk" "8741"
## [107,] "2023-08-11" "actor" "Buk" "8741"
## [108,] "2023-08-11" "receiver" "Nko" "8742"
## [109,] "2022-01-24" "receiver" "Nkun" "8934"
## [110,] "2022-04-01" "receiver" "Ndl" "9474"
## [111,] "2022-08-19" "receiver" "Gil" "9935"
## [112,] "2022-08-23" "receiver" "Nko" "9939"
## [113,] "2022-08-23" "focal" "Ndon" "9941"
## [114,] "2022-08-23" "actor" "Ndon" "9941"
## [115,] "2022-10-15" "focal" "Buk" "10233"
## [116,] "2022-10-15" "actor" "Buk" "10233"
## [117,] "2022-10-15" "receiver" "Ginq" "10233"
## [118,] "2022-10-15" "focal" "Buk" "10234"
## [119,] "2022-10-15" "actor" "Buk" "10234"
## [120,] "2023-03-13" "receiver" "Ndin" "10447"
## [121,] "2023-03-16" "receiver" "Gobe" "10468"
## [122,] "2023-03-16" "receiver" "Giji" "10471"
## [123,] "2023-03-16" "receiver" "Nca" "10472"
## [124,] "2023-07-19" "receiver" "Ndaw" "10650"
## [125,] "2023-07-19" "focal" "Ndaw" "10651"
## [126,] "2023-07-19" "actor" "Ndum" "10651"
## [127,] "2023-07-19" "receiver" "Ndaw" "10651"
## [128,] "2023-08-11" "receiver" "Nda" "10672"
## [129,] "2023-08-11" "focal" "Ndon" "10673"
## [130,] "2023-08-11" "actor" "Ndon" "10673"
## [131,] "2023-08-11" "receiver" "Ndum" "10673"
## [132,] "2023-08-11" "focal" "Ndon" "10674"
## [133,] "2023-08-11" "actor" "Ndon" "10674"
## [134,] "2023-08-11" "receiver" "Nko" "10674"
## [135,] "2023-08-11" "focal" "Buk" "10675"
## [136,] "2023-08-11" "actor" "Buk" "10675"
## [137,] "2023-08-11" "receiver" "Nko" "10675"
## [138,] "2023-09-27" "receiver" "Ndum" "10764"
## [139,] "2022-07-22" "receiver" "Nko" "12995"
## [140,] "2022-08-19" "receiver" "Nda" "13208"
## [141,] "2022-08-19" "receiver" "Vla" "13212"
## [142,] "2022-08-23" "receiver" "Nko" "13229"
## [143,] "2022-08-23" "focal" "Ndon" "13230"
## [144,] "2022-08-23" "actor" "Ndon" "13230"
## [145,] "2022-08-23" "focal" "Ndon" "13231"
## [146,] "2022-08-23" "actor" "Ndon" "13231"
## [147,] "2022-08-23" "actor" "Ndon" "13232"
## [148,] "2022-08-23" "focal" "Ndon" "13233"
## [149,] "2022-08-23" "actor" "Ndon" "13233"
## [150,] "2022-08-23" "receiver" "Godu" "13234"
## [151,] "2022-08-23" "receiver" "Godu" "13237"
## [152,] "2022-08-23" "focal" "Godu" "13240"
## [153,] "2022-08-23" "actor" "Godu" "13240"
## [154,] "2022-08-23" "focal" "Godu" "13241"
## [155,] "2022-08-23" "actor" "Godu" "13241"
## [156,] "2022-08-23" "receiver" "Guba" "13241"
## [157,] "2022-10-15" "focal" "Buk" "13896"
## [158,] "2022-10-15" "actor" "Buk" "13896"
## [159,] "2022-10-15" "focal" "Buk" "13897"
## [160,] "2022-10-15" "actor" "Buk" "13897"
## [161,] "2022-10-15" "focal" "Buk" "13898"
## [162,] "2022-10-15" "actor" "Buk" "13898"
## [163,] "2022-10-15" "focal" "Buk" "13899"
## [164,] "2022-10-15" "actor" "Buk" "13899"
## [165,] "2022-10-15" "focal" "Buk" "13908"
## [166,] "2022-10-15" "actor" "Buk" "13908"
## [167,] "2022-10-15" "focal" "Ginq" "13909"
## [168,] "2022-10-15" "actor" "Ginq" "13909"
## [169,] "2023-03-04" "receiver" "Ndin" "14351"
## [170,] "2023-03-13" "receiver" "Ndin" "14405"
## [171,] "2023-03-16" "receiver" "Gobe" "14440"
## [172,] "2023-03-16" "receiver" "Nca" "14441"
## [173,] "2023-03-16" "receiver" "Nca" "14442"
## [174,] "2023-03-16" "receiver" "Ndl" "14444"
## [175,] "2023-03-16" "receiver" "Nca" "14446"
## [176,] "2023-03-16" "receiver" "Ndin" "14447"
## [177,] "2023-03-16" "receiver" "Ndl" "14450"
## [178,] "2023-03-16" "receiver" "Nko" "14451"
## [179,] "2023-03-16" "receiver" "Giji" "14452"
## [180,] "2023-05-03" "receiver" "Gom" "14714"
## [181,] "2023-05-03" "receiver" "Gom" "14728"
## [182,] "2023-07-19" "focal" "Ndaw" "14874"
## [183,] "2023-07-19" "actor" "Ndaw" "14874"
## [184,] "2023-07-19" "receiver" "Ndum" "14874"
## [185,] "2023-07-19" "receiver" "Ndum" "14875"
## [186,] "2023-07-19" "focal" "Ndaw" "14876"
## [187,] "2023-07-19" "receiver" "Ndaw" "14876"
## [188,] "2023-08-11" "receiver" "Gom" "14934"
## [189,] "2023-08-11" "receiver" "Nda" "14937"
## [190,] "2023-08-11" "receiver" "Nda" "14939"
## [191,] "2023-08-11" "focal" "Nda" "14940"
## [192,] "2023-08-11" "actor" "Nda" "14940"
## [193,] "2023-08-11" "receiver" "Ndin" "14942"
## [194,] "2023-08-11" "focal" "Ndon" "14943"
## [195,] "2023-08-11" "actor" "Ndon" "14943"
## [196,] "2023-08-11" "receiver" "Ndum" "14943"
## [197,] "2023-08-11" "focal" "Ndon" "14944"
## [198,] "2023-08-11" "actor" "Ndon" "14944"
## [199,] "2023-08-11" "focal" "Buk" "14945"
## [200,] "2023-08-11" "actor" "Ndon" "14945"
## [201,] "2023-08-11" "receiver" "Buk" "14945"
## [202,] "2023-08-11" "actor" "Ndon" "14946"
## [203,] "2023-08-11" "focal" "Ndon" "14947"
## [204,] "2023-08-11" "actor" "Ndon" "14947"
## [205,] "2023-08-11" "receiver" "Ndum" "14947"
## [206,] "2023-08-11" "focal" "Ndaw" "14948"
## [207,] "2023-08-11" "actor" "Ndaw" "14948"
## [208,] "2023-08-11" "receiver" "Ndl" "14950"
## [209,] "2023-08-11" "focal" "Ndl" "14951"
## [210,] "2023-08-11" "actor" "Ndl" "14951"
## [211,] "2023-08-11" "focal" "Nko" "14953"
## [212,] "2023-08-11" "actor" "Nko" "14953"
## [213,] "2023-09-27" "focal" "Ndum" "15168"
## [214,] "2023-09-27" "actor" "Ndum" "15168"
## [215,] "2023-09-27" "actor" "Nkos" "15169"
## [216,] "2023-09-27" "receiver" "Gib" "15172"
## [217,] "2023-10-13" "focal" "Buk" "15245"
## [218,] "2023-10-13" "actor" "Buk" "15245"
## [219,] "2023-10-13" "focal" "Buk" "15246"
## [220,] "2023-10-13" "actor" "Buk" "15246"
## [221,] "2023-10-13" "focal" "Buk" "15247"
## [222,] "2023-10-13" "actor" "Buk" "15247"
## some dyads don't have defined CSI/DSI despite co-residency
## some dyads don't have defined CSI/DSI despite co-residency
## some dyads don't have defined CSI/DSI despite co-residency
## some dyads don't have defined CSI/DSI despite co-residency
## dyad rank zDSI
## 1 Ginq_@_Sho 116 -0.24727
## i1 i2 dyad zDSI
## 1 Ncok Sho Ncok_@_Sho 0.12960
## 2 Ndon Sho Ndon_@_Sho 0.05034
## 3 Gubh Sho Gubh_@_Sho -0.10711
## 4 Ghid Sho Ghid_@_Sho -0.11191
## 5 Ndum Sho Ndum_@_Sho -0.14986
## 6 Ndaw Sho Ndaw_@_Sho -0.15429
## 7 Gub Sho Gub_@_Sho -0.16490
## 8 Gil Sho Gil_@_Sho -0.20710
## 9 Nda Sho Nda_@_Sho -0.22794
## 10 Guz Sho Guz_@_Sho -0.23846
## 11 Nkos Sho Nkos_@_Sho -0.23911
## 12 Buk Sho Buk_@_Sho -0.24610
## 13 Ginq Sho Ginq_@_Sho -0.24727
## 14 Gib Sho Gib_@_Sho -0.25075
## 15 Godu Sho Godu_@_Sho -0.25136
## 16 Gom Sho Gom_@_Sho -0.25383
## 17 Gobe Sho Gobe_@_Sho -0.27323
## 18 Nko Sho Nko_@_Sho -0.27869
## 19 Guba Sho Guba_@_Sho -0.28151
## 20 Gha Sho Gha_@_Sho -0.28211
## 21 Sho Vla Sho_@_Vla -0.28352
## i1 i2 dyad zDSI
## 1 Ginq Gom Ginq_@_Gom 1.30046
## 2 Buk Ginq Buk_@_Ginq 1.29729
## 3 Ghid Ginq Ghid_@_Ginq 0.58689
## 4 Ginq Gubh Ginq_@_Gubh 0.15406
## 5 Ginq Godu Ginq_@_Godu 0.04182
## 6 Ginq Ncok Ginq_@_Ncok 0.03422
## 7 Ginq Ndaw Ginq_@_Ndaw -0.01527
## 8 Ginq Guz Ginq_@_Guz -0.01707
## 9 Ginq Guba Ginq_@_Guba -0.05460
## 10 Gha Ginq Gha_@_Ginq -0.08226
## 11 Ginq Vla Ginq_@_Vla -0.13911
## 12 Ginq Ndon Ginq_@_Ndon -0.16404
## 13 Ginq Gub Ginq_@_Gub -0.16963
## 14 Gil Ginq Gil_@_Ginq -0.18138
## 15 Gib Ginq Gib_@_Ginq -0.21068
## 16 Ginq Gobe Ginq_@_Gobe -0.21976
## 17 Ginq Nda Ginq_@_Nda -0.24456
## 18 Ginq Sho Ginq_@_Sho -0.24727
## 19 Ginq Nko Ginq_@_Nko -0.24823
## 20 Ginq Ndum Ginq_@_Ndum -0.25669
## 21 Ginq Nkos Ginq_@_Nkos -0.26255
## some dyads don't have defined CSI/DSI despite co-residency
## Relative Rank of Specified Dyad:
## [1] dyad rank zDSI
## <0 lignes> (ou 'row.names' de longueur nulle)
##
## DSI Values for Interactions Involving Sho :
## i1 i2 dyad zDSI
## 1 Ncok Sho Ncok_@_Sho 0.12960
## 2 Ndon Sho Ndon_@_Sho 0.05034
## 3 Gubh Sho Gubh_@_Sho -0.10711
## 4 Ghid Sho Ghid_@_Sho -0.11191
## 5 Ndum Sho Ndum_@_Sho -0.14986
## 6 Ndaw Sho Ndaw_@_Sho -0.15429
## 7 Gub Sho Gub_@_Sho -0.16490
## 8 Gil Sho Gil_@_Sho -0.20710
## 9 Nda Sho Nda_@_Sho -0.22794
## 10 Guz Sho Guz_@_Sho -0.23846
## 11 Nkos Sho Nkos_@_Sho -0.23911
## 12 Buk Sho Buk_@_Sho -0.24610
## 13 Ginq Sho Ginq_@_Sho -0.24727
## 14 Gib Sho Gib_@_Sho -0.25075
## 15 Godu Sho Godu_@_Sho -0.25136
## 16 Gom Sho Gom_@_Sho -0.25383
## 17 Gobe Sho Gobe_@_Sho -0.27323
## 18 Nko Sho Nko_@_Sho -0.27869
## 19 Guba Sho Guba_@_Sho -0.28151
## 20 Gha Sho Gha_@_Sho -0.28211
## 21 Sho Vla Sho_@_Vla -0.28352
##
## DSI Values for Interactions Involving Ginq :
## i1 i2 dyad zDSI
## 1 Ginq Gom Ginq_@_Gom 1.30046
## 2 Buk Ginq Buk_@_Ginq 1.29729
## 3 Ghid Ginq Ghid_@_Ginq 0.58689
## 4 Ginq Gubh Ginq_@_Gubh 0.15406
## 5 Ginq Godu Ginq_@_Godu 0.04182
## 6 Ginq Ncok Ginq_@_Ncok 0.03422
## 7 Ginq Ndaw Ginq_@_Ndaw -0.01527
## 8 Ginq Guz Ginq_@_Guz -0.01707
## 9 Ginq Guba Ginq_@_Guba -0.05460
## 10 Gha Ginq Gha_@_Ginq -0.08226
## 11 Ginq Vla Ginq_@_Vla -0.13911
## 12 Ginq Ndon Ginq_@_Ndon -0.16404
## 13 Ginq Gub Ginq_@_Gub -0.16963
## 14 Gil Ginq Gil_@_Ginq -0.18138
## 15 Gib Ginq Gib_@_Ginq -0.21068
## 16 Ginq Gobe Ginq_@_Gobe -0.21976
## 17 Ginq Nda Ginq_@_Nda -0.24456
## 18 Ginq Sho Ginq_@_Sho -0.24727
## 19 Ginq Nko Ginq_@_Nko -0.24823
## 20 Ginq Ndum Ginq_@_Ndum -0.25669
## 21 Ginq Nkos Ginq_@_Nkos -0.26255
## some dyads don't have defined CSI/DSI despite co-residency
## Relative Rank of Specified Dyad:
## dyad rank zDSI
## 1 Buk_@_Ndaw 93 -0.21387
##
## DSI Values for Interactions Involving Buk :
## i1 i2 dyad zDSI
## 1 Buk Ginq Buk_@_Ginq 1.33398
## 2 Buk Ncok Buk_@_Ncok -0.08841
## 3 Buk Gubh Buk_@_Gubh -0.12951
## 4 Buk Ghid Buk_@_Ghid -0.14436
## 5 Buk Ndon Buk_@_Ndon -0.14571
## 6 Buk Gil Buk_@_Gil -0.14716
## 7 Buk Nkos Buk_@_Nkos -0.16039
## 8 Buk Guz Buk_@_Guz -0.20408
## 9 Buk Ndaw Buk_@_Ndaw -0.21387
## 10 Buk Gom Buk_@_Gom -0.22386
## 11 Buk Sho Buk_@_Sho -0.23854
## 12 Buk Vla Buk_@_Vla -0.24987
## 13 Buk Nda Buk_@_Nda -0.25026
## 14 Buk Gub Buk_@_Gub -0.25138
## 15 Buk Gha Buk_@_Gha -0.25768
## 16 Buk Godu Buk_@_Godu -0.26854
## 17 Buk Gib Buk_@_Gib -0.27173
## 18 Buk Nko Buk_@_Nko -0.27217
## 19 Buk Ndum Buk_@_Ndum -0.27741
## 20 Buk Gobe Buk_@_Gobe -0.29369
## 21 Buk Guba Buk_@_Guba -0.30204
##
## DSI Values for Interactions Involving Ndaw :
## i1 i2 dyad zDSI
## 1 Ndaw Ndon Ndaw_@_Ndon 0.22580
## 2 Ghid Ndaw Ghid_@_Ndaw 0.11416
## 3 Ncok Ndaw Ncok_@_Ndaw 0.08163
## 4 Ginq Ndaw Ginq_@_Ndaw 0.00889
## 5 Guz Ndaw Guz_@_Ndaw -0.04503
## 6 Ndaw Ndum Ndaw_@_Ndum -0.06061
## 7 Gha Ndaw Gha_@_Ndaw -0.06664
## 8 Nda Ndaw Nda_@_Ndaw -0.07319
## 9 Gil Ndaw Gil_@_Ndaw -0.12238
## 10 Gub Ndaw Gub_@_Ndaw -0.13637
## 11 Gom Ndaw Gom_@_Ndaw -0.14126
## 12 Ndaw Sho Ndaw_@_Sho -0.15421
## 13 Gobe Ndaw Gobe_@_Ndaw -0.16434
## 14 Gubh Ndaw Gubh_@_Ndaw -0.17220
## 15 Ndaw Vla Ndaw_@_Vla -0.19111
## 16 Ndaw Nko Ndaw_@_Nko -0.19502
## 17 Guba Ndaw Guba_@_Ndaw -0.19675
## 18 Ndaw Nkos Ndaw_@_Nkos -0.20852
## 19 Buk Ndaw Buk_@_Ndaw -0.21387
## 20 Gib Ndaw Gib_@_Ndaw -0.22387
## 21 Godu Ndaw Godu_@_Ndaw -0.23676
## some dyads don't have defined CSI/DSI despite co-residency
## Relative Rank of Specified Dyad:
## dyad rank zDSI
## 1 Sey_@_Sirk 1627 -0.13823
##
## DSI Values for Interactions Involving Sey :
## i1 i2 dyad zDSI
## 1 Nuk Sey Nuk_@_Sey 2.04830
## 2 Obse Sey Obse_@_Sey 0.49546
## 3 Naal Sey Naal_@_Sey 0.32224
## 4 Piep Sey Piep_@_Sey 0.22918
## 5 Enge Sey Enge_@_Sey 0.22537
## 6 Riss Sey Riss_@_Sey 0.22128
## 7 Nuu Sey Nuu_@_Sey 0.16281
## 8 Rid Sey Rid_@_Sey 0.11481
## 9 Sey Skem Sey_@_Skem 0.08656
## 10 Nurk Sey Nurk_@_Sey 0.03152
## 11 Eis Sey Eis_@_Sey 0.00934
## 12 Sey Sig Sey_@_Sig -0.01001
## 13 Ouli Sey Ouli_@_Sey -0.02732
## 14 Oerw Sey Oerw_@_Sey -0.05418
## 15 Hee Sey Hee_@_Sey -0.06746
## 16 Reen Sey Reen_@_Sey -0.07363
## 17 Oort Sey Oort_@_Sey -0.07400
## 18 Popp Sey Popp_@_Sey -0.08431
## 19 Non Sey Non_@_Sey -0.09041
## 20 Rede Sey Rede_@_Sey -0.09524
## 21 Rivi Sey Rivi_@_Sey -0.09546
## 22 Pom Sey Pom_@_Sey -0.09636
## 23 Lewe Sey Lewe_@_Sey -0.10435
## 24 Miel Sey Miel_@_Sey -0.10503
## 25 Ree Sey Ree_@_Sey -0.10718
## 26 Pie Sey Pie_@_Sey -0.10971
## 27 Hond Sey Hond_@_Sey -0.11007
## 28 Papp Sey Papp_@_Sey -0.11558
## 29 Reno Sey Reno_@_Sey -0.11558
## 30 Sey War Sey_@_War -0.12044
## 31 Ome Sey Ome_@_Sey -0.12144
## 32 Nak Sey Nak_@_Sey -0.12167
## 33 Oup Sey Oup_@_Sey -0.12340
## 34 Haai Sey Haai_@_Sey -0.12397
## 35 Heli Sey Heli_@_Sey -0.12397
## 36 Rimp Sey Rimp_@_Sey -0.12397
## 37 Puol Sey Puol_@_Sey -0.12625
## 38 Pann Sey Pann_@_Sey -0.12807
## 39 Eina Sey Eina_@_Sey -0.12892
## 40 Sey Xin Sey_@_Xin -0.12927
## 41 Heer Sey Heer_@_Sey -0.12930
## 42 Aal Sey Aal_@_Sey -0.12983
## 43 Aapi Sey Aapi_@_Sey -0.12983
## 44 App Sey App_@_Sey -0.12983
## 45 Misk Sey Misk_@_Sey -0.12983
## 46 Mui Sey Mui_@_Sey -0.12983
## 47 Nooi Sey Nooi_@_Sey -0.12983
## 48 Sey Spe Sey_@_Spe -0.12983
## 49 Sari Sey Sari_@_Sey -0.13026
## 50 Sey Syl Sey_@_Syl -0.13125
## 51 Sey Xia Sey_@_Xia -0.13257
## 52 Kom Sey Kom_@_Sey -0.13273
## 53 Gese Sey Gese_@_Sey -0.13559
## 54 Aan Sey Aan_@_Sey -0.13633
## 55 Apa Sey Apa_@_Sey -0.13823
## 56 Asis Sey Asis_@_Sey -0.13823
## 57 Dok Sey Dok_@_Sey -0.13823
## 58 Flu Sey Flu_@_Sey -0.13823
## 59 Goe Sey Goe_@_Sey -0.13823
## 60 Nge Sey Nge_@_Sey -0.13823
## 61 Numb Sey Numb_@_Sey -0.13823
## 62 Pix Sey Pix_@_Sey -0.13823
## 63 Sey Sirk Sey_@_Sirk -0.13823
## 64 Sey Sitr Sey_@_Sitr -0.13823
## 65 Sey Umb Sey_@_Umb -0.13823
## 66 Sey Unk Sey_@_Unk -0.13823
##
## DSI Values for Interactions Involving Sirk :
## i1 i2 dyad zDSI
## 1 Sirk Sitr Sirk_@_Sitr 7.05968
## 2 Pix Sirk Pix_@_Sirk 5.57675
## 3 Sirk Syl Sirk_@_Syl 5.55807
## 4 Sirk Skem Sirk_@_Skem 3.83152
## 5 Oerw Sirk Oerw_@_Sirk 3.54644
## 6 Ree Sirk Ree_@_Sirk 3.31635
## 7 Papp Sirk Papp_@_Sirk 3.06645
## 8 Sirk Spe Sirk_@_Spe 2.87040
## 9 Non Sirk Non_@_Sirk 2.81655
## 10 Misk Sirk Misk_@_Sirk 2.37208
## 11 Sig Sirk Sig_@_Sirk 2.37208
## 12 Aapi Sirk Aapi_@_Sirk 2.33655
## 13 Rid Sirk Rid_@_Sirk 2.33655
## 14 Pie Sirk Pie_@_Sirk 1.97418
## 15 Lewe Sirk Lewe_@_Sirk 1.58832
## 16 Heli Sirk Heli_@_Sirk 1.33990
## 17 Naal Sirk Naal_@_Sirk 1.33990
## 18 Flu Sirk Flu_@_Sirk 1.09476
## 19 Miel Sirk Miel_@_Sirk 0.41075
## 20 Eina Sirk Eina_@_Sirk 0.40126
## 21 Piep Sirk Piep_@_Sirk 0.37257
## 22 Aal Sirk Aal_@_Sirk 0.36010
## 23 Nooi Sirk Nooi_@_Sirk 0.36010
## 24 Obse Sirk Obse_@_Sirk 0.36010
## 25 Popp Sirk Popp_@_Sirk 0.36010
## 26 Rede Sirk Rede_@_Sirk 0.14906
## 27 Heer Sirk Heer_@_Sirk 0.04753
## 28 Oort Sirk Oort_@_Sirk 0.03657
## 29 Nurk Sirk Nurk_@_Sirk 0.00463
## 30 Riss Sirk Riss_@_Sirk -0.01246
## 31 Enge Sirk Enge_@_Sirk -0.03894
## 32 Nge Sirk Nge_@_Sirk -0.05263
## 33 Apa Sirk Apa_@_Sirk -0.05935
## 34 Pann Sirk Pann_@_Sirk -0.06232
## 35 Ouli Sirk Ouli_@_Sirk -0.06409
## 36 Dok Sirk Dok_@_Sirk -0.07652
## 37 Nuu Sirk Nuu_@_Sirk -0.08265
## 38 Sari Sirk Sari_@_Sirk -0.09068
## 39 Nak Sirk Nak_@_Sirk -0.09869
## 40 Gese Sirk Gese_@_Sirk -0.10295
## 41 Kom Sirk Kom_@_Sirk -0.10500
## 42 Sirk Xia Sirk_@_Xia -0.11532
## 43 Sirk War Sirk_@_War -0.11830
## 44 Sirk Xin Sirk_@_Xin -0.11943
## 45 Pom Sirk Pom_@_Sirk -0.12681
## 46 Sirk Umb Sirk_@_Umb -0.12840
## 47 Aan Sirk Aan_@_Sirk -0.13823
## 48 App Sirk App_@_Sirk -0.13823
## 49 Asis Sirk Asis_@_Sirk -0.13823
## 50 Eis Sirk Eis_@_Sirk -0.13823
## 51 Goe Sirk Goe_@_Sirk -0.13823
## 52 Haai Sirk Haai_@_Sirk -0.13823
## 53 Hee Sirk Hee_@_Sirk -0.13823
## 54 Hond Sirk Hond_@_Sirk -0.13823
## 55 Mui Sirk Mui_@_Sirk -0.13823
## 56 Nuk Sirk Nuk_@_Sirk -0.13823
## 57 Numb Sirk Numb_@_Sirk -0.13823
## 58 Ome Sirk Ome_@_Sirk -0.13823
## 59 Oup Sirk Oup_@_Sirk -0.13823
## 60 Puol Sirk Puol_@_Sirk -0.13823
## 61 Reen Sirk Reen_@_Sirk -0.13823
## 62 Reno Sirk Reno_@_Sirk -0.13823
## 63 Rimp Sirk Rimp_@_Sirk -0.13823
## 64 Rivi Sirk Rivi_@_Sirk -0.13823
## 65 Sey Sirk Sey_@_Sirk -0.13823
## 66 Sirk Unk Sirk_@_Unk -0.13823
## some dyads don't have defined CSI/DSI despite co-residency
## Relative Rank of Specified Dyad:
## [1] dyad rank zDSI
## <0 lignes> (ou 'row.names' de longueur nulle)
##
## DSI Values for Interactions Involving Xia :
## i1 i2 dyad zDSI
## 1 Reen Xia Reen_@_Xia 0.88590
## 2 Nooi Xia Nooi_@_Xia 0.80836
## 3 Eina Xia Eina_@_Xia 0.11689
## 4 Nge Xia Nge_@_Xia 0.02650
## 5 Naal Xia Naal_@_Xia -0.03937
## 6 Sirk Xia Sirk_@_Xia -0.07954
## 7 Hond Xia Hond_@_Xia -0.08624
## 8 Skem Xia Skem_@_Xia -0.09264
## 9 Aal Xia Aal_@_Xia -0.10576
## 10 Sitr Xia Sitr_@_Xia -0.11920
## 11 Nurk Xia Nurk_@_Xia -0.12608
## 12 Aapi Xia Aapi_@_Xia -0.13231
## 13 Misk Xia Misk_@_Xia -0.13231
## 14 Oerw Xia Oerw_@_Xia -0.13231
## 15 Miel Xia Miel_@_Xia -0.13772
## 16 Oort Xia Oort_@_Xia -0.14092
## 17 Heer Xia Heer_@_Xia -0.14216
## 18 Xia Xin Xia_@_Xin -0.14450
## 19 Ouli Xia Ouli_@_Xia -0.14635
## 20 War Xia War_@_Xia -0.14750
## 21 Piep Xia Piep_@_Xia -0.15279
## 22 Puol Xia Puol_@_Xia -0.15287
## 23 Sey Xia Sey_@_Xia -0.16591
## 24 Obse Xia Obse_@_Xia -0.16959
## 25 Sig Xia Sig_@_Xia -0.16959
## 26 Nak Xia Nak_@_Xia -0.17517
## 27 Aan Xia Aan_@_Xia -0.18120
## 28 Apa Xia Apa_@_Xia -0.18120
## 29 App Xia App_@_Xia -0.18120
## 30 Asis Xia Asis_@_Xia -0.18120
## 31 Dok Xia Dok_@_Xia -0.18120
## 32 Eis Xia Eis_@_Xia -0.18120
## 33 Enge Xia Enge_@_Xia -0.18120
## 34 Flu Xia Flu_@_Xia -0.18120
## 35 Gese Xia Gese_@_Xia -0.18120
## 36 Goe Xia Goe_@_Xia -0.18120
## 37 Haai Xia Haai_@_Xia -0.18120
## 38 Hee Xia Hee_@_Xia -0.18120
## 39 Heli Xia Heli_@_Xia -0.18120
## 40 Kom Xia Kom_@_Xia -0.18120
## 41 Lewe Xia Lewe_@_Xia -0.18120
## 42 Mui Xia Mui_@_Xia -0.18120
## 43 Non Xia Non_@_Xia -0.18120
## 44 Nuk Xia Nuk_@_Xia -0.18120
## 45 Numb Xia Numb_@_Xia -0.18120
## 46 Nuu Xia Nuu_@_Xia -0.18120
## 47 Ome Xia Ome_@_Xia -0.18120
## 48 Oup Xia Oup_@_Xia -0.18120
## 49 Pann Xia Pann_@_Xia -0.18120
## 50 Papp Xia Papp_@_Xia -0.18120
## 51 Pix Xia Pix_@_Xia -0.18120
## 52 Pom Xia Pom_@_Xia -0.18120
## 53 Popp Xia Popp_@_Xia -0.18120
## 54 Rede Xia Rede_@_Xia -0.18120
## 55 Ree Xia Ree_@_Xia -0.18120
## 56 Reno Xia Reno_@_Xia -0.18120
## 57 Rid Xia Rid_@_Xia -0.18120
## 58 Rimp Xia Rimp_@_Xia -0.18120
## 59 Riss Xia Riss_@_Xia -0.18120
## 60 Rivi Xia Rivi_@_Xia -0.18120
## 61 Sari Xia Sari_@_Xia -0.18120
## 62 Spe Xia Spe_@_Xia -0.18120
## 63 Syl Xia Syl_@_Xia -0.18120
## 64 Umb Xia Umb_@_Xia -0.18120
## 65 Unk Xia Unk_@_Xia -0.18120
##
## DSI Values for Interactions Involving Piep :
## i1 i2 dyad zDSI
## 1 Piep Pix Piep_@_Pix 2.30107
## 2 Piep Skem Piep_@_Skem 2.08812
## 3 Piep Puol Piep_@_Puol 1.71465
## 4 Piep Rid Piep_@_Rid 1.49592
## 5 Piep Popp Piep_@_Popp 0.97634
## 6 Piep Reen Piep_@_Reen 0.91151
## 7 Eis Piep Eis_@_Piep 0.71233
## 8 Papp Piep Papp_@_Piep 0.61339
## 9 Obse Piep Obse_@_Piep 0.59385
## 10 Oerw Piep Oerw_@_Piep 0.56375
## 11 Piep Rede Piep_@_Rede 0.38967
## 12 Piep Sirk Piep_@_Sirk 0.36117
## 13 Miel Piep Miel_@_Piep 0.24335
## 14 Gese Piep Gese_@_Piep 0.23258
## 15 Piep Sitr Piep_@_Sitr 0.21942
## 16 Aal Piep Aal_@_Piep 0.21277
## 17 Piep Rimp Piep_@_Rimp 0.21277
## 18 Piep Spe Piep_@_Spe 0.21277
## 19 Piep Riss Piep_@_Riss 0.21162
## 20 Piep Sey Piep_@_Sey 0.19640
## 21 Oort Piep Oort_@_Piep 0.17692
## 22 Nak Piep Nak_@_Piep 0.17555
## 23 Goe Piep Goe_@_Piep 0.16977
## 24 Piep Pom Piep_@_Pom 0.15447
## 25 Hond Piep Hond_@_Piep 0.13647
## 26 Piep Sari Piep_@_Sari 0.06764
## 27 Ouli Piep Ouli_@_Piep 0.04677
## 28 Kom Piep Kom_@_Piep 0.01215
## 29 Oup Piep Oup_@_Piep 0.00601
## 30 Hee Piep Hee_@_Piep -0.02419
## 31 Piep Umb Piep_@_Umb -0.02437
## 32 Aan Piep Aan_@_Piep -0.04413
## 33 Pann Piep Pann_@_Piep -0.04825
## 34 Dok Piep Dok_@_Piep -0.06775
## 35 Aapi Piep Aapi_@_Piep -0.08760
## 36 Lewe Piep Lewe_@_Piep -0.08760
## 37 Piep Syl Piep_@_Syl -0.08760
## 38 Flu Piep Flu_@_Piep -0.10870
## 39 Piep War Piep_@_War -0.13046
## 40 Eina Piep Eina_@_Piep -0.13469
## 41 Nge Piep Nge_@_Piep -0.13495
## 42 Piep Xin Piep_@_Xin -0.14080
## 43 Piep Xia Piep_@_Xia -0.15279
## 44 Apa Piep Apa_@_Piep -0.18120
## 45 App Piep App_@_Piep -0.18120
## 46 Asis Piep Asis_@_Piep -0.18120
## 47 Enge Piep Enge_@_Piep -0.18120
## 48 Haai Piep Haai_@_Piep -0.18120
## 49 Heer Piep Heer_@_Piep -0.18120
## 50 Heli Piep Heli_@_Piep -0.18120
## 51 Misk Piep Misk_@_Piep -0.18120
## 52 Mui Piep Mui_@_Piep -0.18120
## 53 Naal Piep Naal_@_Piep -0.18120
## 54 Non Piep Non_@_Piep -0.18120
## 55 Nooi Piep Nooi_@_Piep -0.18120
## 56 Nuk Piep Nuk_@_Piep -0.18120
## 57 Numb Piep Numb_@_Piep -0.18120
## 58 Nurk Piep Nurk_@_Piep -0.18120
## 59 Nuu Piep Nuu_@_Piep -0.18120
## 60 Ome Piep Ome_@_Piep -0.18120
## 61 Piep Ree Piep_@_Ree -0.18120
## 62 Piep Reno Piep_@_Reno -0.18120
## 63 Piep Rivi Piep_@_Rivi -0.18120
## 64 Piep Sig Piep_@_Sig -0.18120
## 65 Piep Unk Piep_@_Unk -0.18120
## some dyads don't have defined CSI/DSI despite co-residency
## Relative Rank of Specified Dyad:
## dyad rank zDSI
## 1 Nge_@_Oerw 307 0.0483
##
## DSI Values for Interactions Involving Nge :
## i1 i2 dyad zDSI
## 1 Nge Reen Nge_@_Reen 0.18284
## 2 Nge Nurk Nge_@_Nurk 0.04903
## 3 Nge Oerw Nge_@_Oerw 0.04830
## 4 Nge Sari Nge_@_Sari 0.04401
## 5 Nge Sirk Nge_@_Sirk 0.03686
## 6 Eina Nge Eina_@_Nge 0.01741
## 7 Nge Piep Nge_@_Piep 0.00534
## 8 Naal Nge Naal_@_Nge -0.00841
## 9 Nge Xia Nge_@_Xia -0.00998
## 10 Enge Nge Enge_@_Nge -0.01037
## 11 Nge Pom Nge_@_Pom -0.02736
## 12 Miel Nge Miel_@_Nge -0.03379
## 13 Aapi Nge Aapi_@_Nge -0.03980
## 14 Mui Nge Mui_@_Nge -0.04515
## 15 Nge Skem Nge_@_Skem -0.05617
## 16 Nge Puol Nge_@_Puol -0.05735
## 17 Nge Rede Nge_@_Rede -0.06024
## 18 Nge Obse Nge_@_Obse -0.06763
## 19 Nge Rid Nge_@_Rid -0.06930
## 20 Heer Nge Heer_@_Nge -0.07605
## 21 Nge Popp Nge_@_Popp -0.07815
## 22 Hond Nge Hond_@_Nge -0.09124
## 23 Nak Nge Nak_@_Nge -0.09436
## 24 Goe Nge Goe_@_Nge -0.10159
## 25 Nge Sig Nge_@_Sig -0.10510
## 26 Nge Pann Nge_@_Pann -0.11363
## 27 Nge Riss Nge_@_Riss -0.11938
## 28 Nge War Nge_@_War -0.11980
## 29 Lewe Nge Lewe_@_Nge -0.12154
## 30 App Nge App_@_Nge -0.12631
## 31 Nge Spe Nge_@_Spe -0.12631
## 32 Kom Nge Kom_@_Nge -0.14036
## 33 Nge Oort Nge_@_Oort -0.14730
## 34 Nge Numb Nge_@_Numb -0.15606
## 35 Flu Nge Flu_@_Nge -0.16058
## 36 Eis Nge Eis_@_Nge -0.16352
## 37 Heli Nge Heli_@_Nge -0.17258
## 38 Nge Reno Nge_@_Reno -0.17258
## 39 Nge Oup Nge_@_Oup -0.17471
## 40 Nge Pix Nge_@_Pix -0.17471
## 41 Aal Nge Aal_@_Nge -0.18377
## 42 Nge Papp Nge_@_Papp -0.18377
## 43 Nge Sitr Nge_@_Sitr -0.18543
## 44 Nge Syl Nge_@_Syl -0.18543
## 45 Nge Ouli Nge_@_Ouli -0.18788
## 46 Nge Umb Nge_@_Umb -0.19089
## 47 Nge Xin Nge_@_Xin -0.19399
## 48 Aan Nge Aan_@_Nge -0.19497
## 49 Asis Nge Asis_@_Nge -0.20073
## 50 Apa Nge Apa_@_Nge -0.21884
## 51 Dok Nge Dok_@_Nge -0.21884
## 52 Gese Nge Gese_@_Nge -0.21884
## 53 Haai Nge Haai_@_Nge -0.21884
## 54 Hee Nge Hee_@_Nge -0.21884
## 55 Misk Nge Misk_@_Nge -0.21884
## 56 Nge Non Nge_@_Non -0.21884
## 57 Nge Nooi Nge_@_Nooi -0.21884
## 58 Nge Nuk Nge_@_Nuk -0.21884
## 59 Nge Nuu Nge_@_Nuu -0.21884
## 60 Nge Ome Nge_@_Ome -0.21884
## 61 Nge Ree Nge_@_Ree -0.21884
## 62 Nge Rimp Nge_@_Rimp -0.21884
## 63 Nge Rivi Nge_@_Rivi -0.21884
## 64 Nge Sey Nge_@_Sey -0.21884
## 65 Nge Unk Nge_@_Unk -0.21884
##
## DSI Values for Interactions Involving Oerw :
## i1 i2 dyad zDSI
## 1 Oerw Oup Oerw_@_Oup 3.89023
## 2 Obse Oerw Obse_@_Oerw 3.00544
## 3 Oerw Popp Oerw_@_Popp 2.42334
## 4 Oerw Ome Oerw_@_Ome 2.17448
## 5 Oerw Pix Oerw_@_Pix 1.88224
## 6 Hond Oerw Hond_@_Oerw 1.69670
## 7 Aapi Oerw Aapi_@_Oerw 1.53485
## 8 Oerw Sig Oerw_@_Sig 1.40446
## 9 Oerw Ouli Oerw_@_Ouli 1.15897
## 10 Goe Oerw Goe_@_Oerw 1.00649
## 11 Oerw Oort Oerw_@_Oort 0.90601
## 12 Aal Oerw Aal_@_Oerw 0.86336
## 13 Oerw Puol Oerw_@_Puol 0.58512
## 14 Lewe Oerw Lewe_@_Oerw 0.47813
## 15 Nuk Oerw Nuk_@_Oerw 0.47813
## 16 Oerw Rid Oerw_@_Rid 0.45863
## 17 Heli Oerw Heli_@_Oerw 0.44589
## 18 Oerw Papp Oerw_@_Papp 0.42080
## 19 Misk Oerw Misk_@_Oerw 0.30952
## 20 Mui Oerw Mui_@_Oerw 0.30952
## 21 Oerw Ree Oerw_@_Ree 0.30952
## 22 Oerw Syl Oerw_@_Syl 0.30952
## 23 Flu Oerw Flu_@_Oerw 0.26665
## 24 Oerw Sirk Oerw_@_Sirk 0.16575
## 25 Oerw Piep Oerw_@_Piep 0.13804
## 26 Nurk Oerw Nurk_@_Oerw 0.12979
## 27 Oerw Pom Oerw_@_Pom 0.10203
## 28 Nak Oerw Nak_@_Oerw 0.08630
## 29 Nge Oerw Nge_@_Oerw 0.04830
## 30 Oerw Sari Oerw_@_Sari 0.04523
## 31 Eina Oerw Eina_@_Oerw -0.05145
## 32 Gese Oerw Gese_@_Oerw -0.06286
## 33 Kom Oerw Kom_@_Oerw -0.07635
## 34 App Oerw App_@_Oerw -0.08247
## 35 Oerw Reen Oerw_@_Reen -0.08247
## 36 Oerw Rivi Oerw_@_Rivi -0.08247
## 37 Oerw Sitr Oerw_@_Sitr -0.08247
## 38 Oerw Skem Oerw_@_Skem -0.08247
## 39 Oerw Sey Oerw_@_Sey -0.09895
## 40 Nuu Oerw Nuu_@_Oerw -0.11206
## 41 Hee Oerw Hee_@_Oerw -0.12387
## 42 Heer Oerw Heer_@_Oerw -0.12850
## 43 Oerw Xin Oerw_@_Xin -0.14330
## 44 Oerw Umb Oerw_@_Umb -0.15093
## 45 Dok Oerw Dok_@_Oerw -0.15309
## 46 Oerw War Oerw_@_War -0.15862
## 47 Oerw Xia Oerw_@_Xia -0.16602
## 48 Aan Oerw Aan_@_Oerw -0.18504
## 49 Apa Oerw Apa_@_Oerw -0.21884
## 50 Asis Oerw Asis_@_Oerw -0.21884
## 51 Eis Oerw Eis_@_Oerw -0.21884
## 52 Enge Oerw Enge_@_Oerw -0.21884
## 53 Haai Oerw Haai_@_Oerw -0.21884
## 54 Miel Oerw Miel_@_Oerw -0.21884
## 55 Naal Oerw Naal_@_Oerw -0.21884
## 56 Non Oerw Non_@_Oerw -0.21884
## 57 Nooi Oerw Nooi_@_Oerw -0.21884
## 58 Numb Oerw Numb_@_Oerw -0.21884
## 59 Oerw Pann Oerw_@_Pann -0.21884
## 60 Oerw Rede Oerw_@_Rede -0.21884
## 61 Oerw Reno Oerw_@_Reno -0.21884
## 62 Oerw Rimp Oerw_@_Rimp -0.21884
## 63 Oerw Riss Oerw_@_Riss -0.21884
## 64 Oerw Spe Oerw_@_Spe -0.21884
## 65 Oerw Unk Oerw_@_Unk -0.21884
## some dyads don't have defined CSI/DSI despite co-residency
## Relative Rank of Specified Dyad:
## [1] dyad rank zDSI
## <0 lignes> (ou 'row.names' de longueur nulle)
##
## DSI Values for Interactions Involving Xin :
## i1 i2 dyad zDSI
## 1 Eina Xin Eina_@_Xin 0.14408
## 2 Sari Xin Sari_@_Xin -0.05816
## 3 Pix Xin Pix_@_Xin -0.06234
## 4 Naal Xin Naal_@_Xin -0.08072
## 5 Papp Xin Papp_@_Xin -0.08220
## 6 Sig Xin Sig_@_Xin -0.08220
## 7 Umb Xin Umb_@_Xin -0.08495
## 8 Nooi Xin Nooi_@_Xin -0.08592
## 9 Oort Xin Oort_@_Xin -0.11589
## 10 Kom Xin Kom_@_Xin -0.12313
## 11 Aapi Xin Aapi_@_Xin -0.12839
## 12 Heli Xin Heli_@_Xin -0.13275
## 13 Ree Xin Ree_@_Xin -0.13275
## 14 Reen Xin Reen_@_Xin -0.13541
## 15 Oerw Xin Oerw_@_Xin -0.13946
## 16 Pom Xin Pom_@_Xin -0.14323
## 17 Obse Xin Obse_@_Xin -0.15113
## 18 Reno Xin Reno_@_Xin -0.15626
## 19 Rivi Xin Rivi_@_Xin -0.15626
## 20 Aan Xin Aan_@_Xin -0.15710
## 21 Xia Xin Xia_@_Xin -0.16869
## 22 Nurk Xin Nurk_@_Xin -0.16912
## 23 Pann Xin Pann_@_Xin -0.17156
## 24 Enge Xin Enge_@_Xin -0.17358
## 25 Piep Xin Piep_@_Xin -0.17358
## 26 Puol Xin Puol_@_Xin -0.17535
## 27 Skem Xin Skem_@_Xin -0.17535
## 28 Sirk Xin Sirk_@_Xin -0.17550
## 29 Riss Xin Riss_@_Xin -0.17693
## 30 Miel Xin Miel_@_Xin -0.17890
## 31 App Xin App_@_Xin -0.17971
## 32 Goe Xin Goe_@_Xin -0.17971
## 33 Numb Xin Numb_@_Xin -0.17971
## 34 Hee Xin Hee_@_Xin -0.18260
## 35 Lewe Xin Lewe_@_Xin -0.18306
## 36 Hond Xin Hond_@_Xin -0.18308
## 37 Flu Xin Flu_@_Xin -0.18589
## 38 Apa Xin Apa_@_Xin -0.18653
## 39 Rede Xin Rede_@_Xin -0.18853
## 40 Syl Xin Syl_@_Xin -0.18964
## 41 Sey Xin Sey_@_Xin -0.19398
## 42 Nge Xin Nge_@_Xin -0.19407
## 43 Mui Xin Mui_@_Xin -0.21238
## 44 Ouli Xin Ouli_@_Xin -0.21429
## 45 Aal Xin Aal_@_Xin -0.22667
## 46 Asis Xin Asis_@_Xin -0.22667
## 47 Dok Xin Dok_@_Xin -0.22667
## 48 Eis Xin Eis_@_Xin -0.22667
## 49 Gese Xin Gese_@_Xin -0.22667
## 50 Haai Xin Haai_@_Xin -0.22667
## 51 Heer Xin Heer_@_Xin -0.22667
## 52 Misk Xin Misk_@_Xin -0.22667
## 53 Nak Xin Nak_@_Xin -0.22667
## 54 Non Xin Non_@_Xin -0.22667
## 55 Nuk Xin Nuk_@_Xin -0.22667
## 56 Nuu Xin Nuu_@_Xin -0.22667
## 57 Ome Xin Ome_@_Xin -0.22667
## 58 Oup Xin Oup_@_Xin -0.22667
## 59 Popp Xin Popp_@_Xin -0.22667
## 60 Rid Xin Rid_@_Xin -0.22667
## 61 Rimp Xin Rimp_@_Xin -0.22667
## 62 Sitr Xin Sitr_@_Xin -0.22667
## 63 Spe Xin Spe_@_Xin -0.22667
## 64 Unk Xin Unk_@_Xin -0.22667
## 65 War Xin War_@_Xin -0.22667
##
## DSI Values for Interactions Involving Ouli :
## i1 i2 dyad zDSI
## 1 Ouli Oup Ouli_@_Oup 2.75296
## 2 Obse Ouli Obse_@_Ouli 1.53088
## 3 Aapi Ouli Aapi_@_Ouli 1.15652
## 4 Haai Ouli Haai_@_Ouli 1.08349
## 5 Ouli Ree Ouli_@_Ree 1.07889
## 6 Goe Ouli Goe_@_Ouli 0.89777
## 7 Ouli Sig Ouli_@_Sig 0.68182
## 8 Oort Ouli Oort_@_Ouli 0.61402
## 9 Oerw Ouli Oerw_@_Ouli 0.52929
## 10 Aal Ouli Aal_@_Ouli 0.47144
## 11 Ouli Umb Ouli_@_Umb 0.44640
## 12 Flu Ouli Flu_@_Ouli 0.32559
## 13 Ouli Sey Ouli_@_Sey 0.25403
## 14 Nuk Ouli Nuk_@_Ouli 0.23083
## 15 Ouli Pann Ouli_@_Pann 0.23083
## 16 Ouli Popp Ouli_@_Popp 0.23083
## 17 Nurk Ouli Nurk_@_Ouli 0.20889
## 18 Hond Ouli Hond_@_Ouli 0.19843
## 19 Ouli Papp Ouli_@_Papp 0.17134
## 20 Lewe Ouli Lewe_@_Ouli 0.14039
## 21 Ouli Piep Ouli_@_Piep 0.09514
## 22 Eis Ouli Eis_@_Ouli 0.07848
## 23 Mui Ouli Mui_@_Ouli 0.07848
## 24 Ouli Skem Ouli_@_Skem 0.07848
## 25 Ouli Riss Ouli_@_Riss 0.07687
## 26 Ouli Reen Ouli_@_Reen 0.03871
## 27 Ouli Pom Ouli_@_Pom 0.01590
## 28 Miel Ouli Miel_@_Ouli 0.00227
## 29 Apa Ouli Apa_@_Ouli -0.00914
## 30 Nooi Ouli Nooi_@_Ouli -0.02327
## 31 Ouli Sari Ouli_@_Sari -0.02332
## 32 Nak Ouli Nak_@_Ouli -0.03175
## 33 App Ouli App_@_Ouli -0.04095
## 34 Ouli War Ouli_@_War -0.04513
## 35 Aan Ouli Aan_@_Ouli -0.04688
## 36 Dok Ouli Dok_@_Ouli -0.07485
## 37 Heer Ouli Heer_@_Ouli -0.08750
## 38 Ouli Sirk Ouli_@_Sirk -0.09285
## 39 Hee Ouli Hee_@_Ouli -0.09856
## 40 Nuu Ouli Nuu_@_Ouli -0.10157
## 41 Kom Ouli Kom_@_Ouli -0.12121
## 42 Asis Ouli Asis_@_Ouli -0.13381
## 43 Heli Ouli Heli_@_Ouli -0.13381
## 44 Ome Ouli Ome_@_Ouli -0.13381
## 45 Ouli Rid Ouli_@_Rid -0.13381
## 46 Gese Ouli Gese_@_Ouli -0.14203
## 47 Eina Ouli Eina_@_Ouli -0.15249
## 48 Ouli Xia Ouli_@_Xia -0.17113
## 49 Nge Ouli Nge_@_Ouli -0.19092
## 50 Ouli Xin Ouli_@_Xin -0.21429
## 51 Enge Ouli Enge_@_Ouli -0.22667
## 52 Misk Ouli Misk_@_Ouli -0.22667
## 53 Naal Ouli Naal_@_Ouli -0.22667
## 54 Non Ouli Non_@_Ouli -0.22667
## 55 Numb Ouli Numb_@_Ouli -0.22667
## 56 Ouli Pix Ouli_@_Pix -0.22667
## 57 Ouli Puol Ouli_@_Puol -0.22667
## 58 Ouli Rede Ouli_@_Rede -0.22667
## 59 Ouli Reno Ouli_@_Reno -0.22667
## 60 Ouli Rimp Ouli_@_Rimp -0.22667
## 61 Ouli Rivi Ouli_@_Rivi -0.22667
## 62 Ouli Sitr Ouli_@_Sitr -0.22667
## 63 Ouli Spe Ouli_@_Spe -0.22667
## 64 Ouli Syl Ouli_@_Syl -0.22667
## 65 Ouli Unk Ouli_@_Unk -0.22667
## some dyads don't have defined CSI/DSI despite co-residency
## Relative Rank of Specified Dyad:
## dyad rank zDSI
## 1 Kom_@_Oort 1646 -0.24291
##
## DSI Values for Interactions Involving Kom :
## i1 i2 dyad zDSI
## 1 Kom Numb Kom_@_Numb 0.34287
## 2 Kom Sari Kom_@_Sari 0.28621
## 3 Kom Oup Kom_@_Oup 0.28342
## 4 Kom Rivi Kom_@_Rivi 0.26104
## 5 Kom Piep Kom_@_Piep 0.24222
## 6 Kom Syl Kom_@_Syl 0.22297
## 7 Kom Nuk Kom_@_Nuk 0.22020
## 8 Kom Ree Kom_@_Ree 0.22020
## 9 Kom Sirk Kom_@_Sirk 0.20936
## 10 Kom Pix Kom_@_Pix 0.17771
## 11 Kom Sig Kom_@_Sig 0.12511
## 12 Kom Nuu Kom_@_Nuu 0.10568
## 13 Kom Oerw Kom_@_Oerw 0.08950
## 14 Kom Reen Kom_@_Reen 0.07480
## 15 Kom Nurk Kom_@_Nurk 0.06062
## 16 Kom Non Kom_@_Non 0.04864
## 17 Kom Sitr Kom_@_Sitr 0.03045
## 18 Hond Kom Hond_@_Kom 0.02823
## 19 Kom Papp Kom_@_Papp 0.01283
## 20 Kom Miel Kom_@_Miel 0.00625
## 21 App Kom App_@_Kom 0.00530
## 22 Kom Puol Kom_@_Puol -0.00023
## 23 Aal Kom Aal_@_Kom -0.02834
## 24 Kom Spe Kom_@_Spe -0.03277
## 25 Kom Ome Kom_@_Ome -0.03554
## 26 Kom Riss Kom_@_Riss -0.05090
## 27 Kom Pom Kom_@_Pom -0.05512
## 28 Kom Rid Kom_@_Rid -0.05792
## 29 Asis Kom Asis_@_Kom -0.06166
## 30 Kom Pann Kom_@_Pann -0.09775
## 31 Kom Rede Kom_@_Rede -0.09890
## 32 Kom Umb Kom_@_Umb -0.10905
## 33 Kom Lewe Kom_@_Lewe -0.11364
## 34 Kom Mui Kom_@_Mui -0.11648
## 35 Kom Sey Kom_@_Sey -0.12138
## 36 Hee Kom Hee_@_Kom -0.12357
## 37 Kom Obse Kom_@_Obse -0.13021
## 38 Flu Kom Flu_@_Kom -0.13838
## 39 Kom Nak Kom_@_Nak -0.13962
## 40 Kom Nge Kom_@_Nge -0.16462
## 41 Kom Xin Kom_@_Xin -0.16585
## 42 Kom Popp Kom_@_Popp -0.17970
## 43 Kom Naal Kom_@_Naal -0.18622
## 44 Kom Nooi Kom_@_Nooi -0.18981
## 45 Kom Ouli Kom_@_Ouli -0.20694
## 46 Aapi Kom Aapi_@_Kom -0.20986
## 47 Aan Kom Aan_@_Kom -0.24291
## 48 Ask Kom Ask_@_Kom -0.24291
## 49 Eina Kom Eina_@_Kom -0.24291
## 50 Eis Kom Eis_@_Kom -0.24291
## 51 Enge Kom Enge_@_Kom -0.24291
## 52 Gese Kom Gese_@_Kom -0.24291
## 53 Goe Kom Goe_@_Kom -0.24291
## 54 Haai Kom Haai_@_Kom -0.24291
## 55 Heer Kom Heer_@_Kom -0.24291
## 56 Heli Kom Heli_@_Kom -0.24291
## 57 Hem Kom Hem_@_Kom -0.24291
## 58 Kno Kom Kno_@_Kom -0.24291
## 59 Kom Misk Kom_@_Misk -0.24291
## 60 Kom Oort Kom_@_Oort -0.24291
## 61 Kom Pie Kom_@_Pie -0.24291
## 62 Kom Reno Kom_@_Reno -0.24291
## 63 Kom Rimp Kom_@_Rimp -0.24291
## 64 Kom Skem Kom_@_Skem -0.24291
## 65 Kom Unk Kom_@_Unk -0.24291
## 66 Kom War Kom_@_War -0.24291
## 67 Kom Xia Kom_@_Xia -0.24291
##
## DSI Values for Interactions Involving Oort :
## i1 i2 dyad zDSI
## 1 Oerw Oort Oerw_@_Oort 3.59386
## 2 Oort Ouli Oort_@_Ouli 2.39347
## 3 Ome Oort Ome_@_Oort 1.76212
## 4 Obse Oort Obse_@_Oort 1.39581
## 5 Lewe Oort Lewe_@_Oort 0.77579
## 6 Oort Oup Oort_@_Oup 0.64386
## 7 Oort Puol Oort_@_Puol 0.59016
## 8 Oort Sey Oort_@_Sey 0.27689
## 9 Oort Rid Oort_@_Rid 0.14215
## 10 Hond Oort Hond_@_Oort 0.12945
## 11 Oort Rede Oort_@_Rede 0.10131
## 12 Flu Oort Flu_@_Oort 0.06660
## 13 Oort Piep Oort_@_Piep 0.03664
## 14 Miel Oort Miel_@_Oort 0.01119
## 15 Mui Oort Mui_@_Oort 0.00409
## 16 Nak Oort Nak_@_Oort 0.00259
## 17 Oort Reen Oort_@_Reen 0.00005
## 18 Oort Sig Oort_@_Sig -0.01354
## 19 Aapi Oort Aapi_@_Oort -0.02174
## 20 Goe Oort Goe_@_Oort -0.02321
## 21 Nuu Oort Nuu_@_Oort -0.02364
## 22 Oort Pann Oort_@_Pann -0.02500
## 23 Aan Oort Aan_@_Oort -0.03623
## 24 Heer Oort Heer_@_Oort -0.05615
## 25 Oort Pom Oort_@_Pom -0.05621
## 26 Gese Oort Gese_@_Oort -0.07370
## 27 Hee Oort Hee_@_Oort -0.07945
## 28 Enge Oort Enge_@_Oort -0.08225
## 29 Haai Oort Haai_@_Oort -0.09162
## 30 Oort Popp Oort_@_Popp -0.09162
## 31 Oort Syl Oort_@_Syl -0.09162
## 32 Nooi Oort Nooi_@_Oort -0.09226
## 33 Oort Xia Oort_@_Xia -0.11346
## 34 Eis Oort Eis_@_Oort -0.11859
## 35 Eina Oort Eina_@_Oort -0.11997
## 36 Misk Oort Misk_@_Oort -0.12023
## 37 Aal Oort Aal_@_Oort -0.13127
## 38 Oort Sirk Oort_@_Sirk -0.16467
## 39 Oort Umb Oort_@_Umb -0.16600
## 40 Oort Xin Oort_@_Xin -0.17493
## 41 Oort Sari Oort_@_Sari -0.18103
## 42 Nurk Oort Nurk_@_Oort -0.19096
## 43 Oort Papp Oort_@_Papp -0.20552
## 44 Oort Rimp Oort_@_Rimp -0.20552
## 45 Oort Rivi Oort_@_Rivi -0.20552
## 46 Oort Sitr Oort_@_Sitr -0.20552
## 47 Nge Oort Nge_@_Oort -0.22418
## 48 App Oort App_@_Oort -0.24291
## 49 Asis Oort Asis_@_Oort -0.24291
## 50 Ask Oort Ask_@_Oort -0.24291
## 51 Heli Oort Heli_@_Oort -0.24291
## 52 Hem Oort Hem_@_Oort -0.24291
## 53 Kno Oort Kno_@_Oort -0.24291
## 54 Kom Oort Kom_@_Oort -0.24291
## 55 Naal Oort Naal_@_Oort -0.24291
## 56 Non Oort Non_@_Oort -0.24291
## 57 Nuk Oort Nuk_@_Oort -0.24291
## 58 Numb Oort Numb_@_Oort -0.24291
## 59 Oort Pie Oort_@_Pie -0.24291
## 60 Oort Pix Oort_@_Pix -0.24291
## 61 Oort Ree Oort_@_Ree -0.24291
## 62 Oort Reno Oort_@_Reno -0.24291
## 63 Oort Riss Oort_@_Riss -0.24291
## 64 Oort Skem Oort_@_Skem -0.24291
## 65 Oort Spe Oort_@_Spe -0.24291
## 66 Oort Unk Oort_@_Unk -0.24291
## 67 Oort War Oort_@_War -0.24291
## some dyads don't have defined CSI/DSI despite co-residency
## Relative Rank of Specified Dyad:
## dyad rank zDSI
## 1 Pom_@_Xian 813 -0.32524
##
## DSI Values for Interactions Involving Pom :
## i1 i2 dyad zDSI
## 1 Pom Prai Pom_@_Prai 0.52191
## 2 Dak Pom Dak_@_Pom 0.27290
## 3 Pom Prim Pom_@_Prim 0.13741
## 4 Pom Ren Pom_@_Ren 0.13741
## 5 Pom Vul Pom_@_Vul 0.04341
## 6 Pom Ram Pom_@_Ram 0.01625
## 7 Pom Tch Pom_@_Tch -0.14947
## 8 Pom Regi Pom_@_Regi -0.16089
## 9 Pom Roc Pom_@_Roc -0.19347
## 10 Pom Xinp Pom_@_Xinp -0.20627
## 11 Beir Pom Beir_@_Pom -0.32524
## 12 Bela Pom Bela_@_Pom -0.32524
## 13 Bet Pom Bet_@_Pom -0.32524
## 14 Fle Pom Fle_@_Pom -0.32524
## 15 Gale Pom Gale_@_Pom -0.32524
## 16 Gaya Pom Gaya_@_Pom -0.32524
## 17 Gran Pom Gran_@_Pom -0.32524
## 18 Gree Pom Gree_@_Pom -0.32524
## 19 Gri Pom Gri_@_Pom -0.32524
## 20 Grim Pom Grim_@_Pom -0.32524
## 21 Gris Pom Gris_@_Pom -0.32524
## 22 Griv Pom Griv_@_Pom -0.32524
## 23 Gua Pom Gua_@_Pom -0.32524
## 24 Guat Pom Guat_@_Pom -0.32524
## 25 Kny Pom Kny_@_Pom -0.32524
## 26 Palm Pom Palm_@_Pom -0.32524
## 27 Pom Prat Pom_@_Prat -0.32524
## 28 Pom Pret Pom_@_Pret -0.32524
## 29 Pom Raba Pom_@_Raba -0.32524
## 30 Pom Rafa Pom_@_Rafa -0.32524
## 31 Pom Rim Pom_@_Rim -0.32524
## 32 Pom Rioj Pom_@_Rioj -0.32524
## 33 Pom Ris Pom_@_Ris -0.32524
## 34 Pom Riva Pom_@_Riva -0.32524
## 35 Pom Tam Pom_@_Tam -0.32524
## 36 Pom Udup Pom_@_Udup -0.32524
## 37 Pom Uls Pom_@_Uls -0.32524
## 38 Pom Upps Pom_@_Upps -0.32524
## 39 Pom Utic Pom_@_Utic -0.32524
## 40 Pom Utr Pom_@_Utr -0.32524
## 41 Pom Vry Pom_@_Vry -0.32524
## 42 Pom Xala Pom_@_Xala -0.32524
## 43 Pom Xar Pom_@_Xar -0.32524
## 44 Pom Xeni Pom_@_Xeni -0.32524
## 45 Pom Xian Pom_@_Xian -0.32524
## 46 Pom Xih Pom_@_Xih -0.32524
## 47 Pom Xop Pom_@_Xop -0.32524
##
## DSI Values for Interactions Involving Xian :
## i1 i2 dyad zDSI
## 1 Xeni Xian Xeni_@_Xian 1.40139
## 2 Prat Xian Prat_@_Xian 1.25600
## 3 Xian Xih Xian_@_Xih 1.22527
## 4 Xar Xian Xar_@_Xian 0.36315
## 5 Gaya Xian Gaya_@_Xian 0.33304
## 6 Vry Xian Vry_@_Xian 0.28773
## 7 Vul Xian Vul_@_Xian 0.22195
## 8 Xala Xian Xala_@_Xian 0.19646
## 9 Rioj Xian Rioj_@_Xian 0.11760
## 10 Palm Xian Palm_@_Xian 0.10509
## 11 Xian Xinp Xian_@_Xinp 0.09001
## 12 Uls Xian Uls_@_Xian 0.08436
## 13 Griv Xian Griv_@_Xian 0.07623
## 14 Udup Xian Udup_@_Xian 0.02791
## 15 Bet Xian Bet_@_Xian 0.01999
## 16 Utic Xian Utic_@_Xian 0.01999
## 17 Gris Xian Gris_@_Xian 0.00389
## 18 Gua Xian Gua_@_Xian -0.04560
## 19 Grim Xian Grim_@_Xian -0.05759
## 20 Xian Xop Xian_@_Xop -0.05759
## 21 Pret Xian Pret_@_Xian -0.06671
## 22 Ram Xian Ram_@_Xian -0.08146
## 23 Guat Xian Guat_@_Xian -0.09089
## 24 Bela Xian Bela_@_Xian -0.09840
## 25 Gran Xian Gran_@_Xian -0.10143
## 26 Tch Xian Tch_@_Xian -0.10290
## 27 Prai Xian Prai_@_Xian -0.10539
## 28 Upps Xian Upps_@_Xian -0.11582
## 29 Gree Xian Gree_@_Xian -0.11606
## 30 Ren Xian Ren_@_Xian -0.11606
## 31 Utr Xian Utr_@_Xian -0.16273
## 32 Gale Xian Gale_@_Xian -0.19141
## 33 Gri Xian Gri_@_Xian -0.19141
## 34 Rim Xian Rim_@_Xian -0.19141
## 35 Regi Xian Regi_@_Xian -0.19857
## 36 Beir Xian Beir_@_Xian -0.23086
## 37 Fle Xian Fle_@_Xian -0.24256
## 38 Dak Xian Dak_@_Xian -0.27596
## 39 Kny Xian Kny_@_Xian -0.32524
## 40 Pom Xian Pom_@_Xian -0.32524
## 41 Prim Xian Prim_@_Xian -0.32524
## 42 Raba Xian Raba_@_Xian -0.32524
## 43 Rafa Xian Rafa_@_Xian -0.32524
## 44 Ris Xian Ris_@_Xian -0.32524
## 45 Riva Xian Riva_@_Xian -0.32524
## 46 Roc Xian Roc_@_Xian -0.32524
## 47 Tam Xian Tam_@_Xian -0.32524
## [1] "DaysSinceStart" "ExperimentDay" "Date"
## [4] "Month" "DyadDay" "TrialDay"
## [7] "Trial" "Time" "Hour"
## [10] "Period" "Male" "Female"
## [13] "Dyad" "DyadDistance" "Proximity"
## [16] "DyadResponse" "IDIndividual1" "Intrusion"
## [19] "IntruderID" "SpecialBehaviour" "Audience"
## [22] "AmountAudience" "Context" "AudienceComposition"
## [25] "KinPresent" "Gender" "ToleranceBinomial"
## [28] "predicted_tolerance"
## [1] "DaysSinceStart" "ExperimentDay" "Date"
## [4] "Month" "DyadDay" "TrialDay"
## [7] "Trial" "Time" "Hour"
## [10] "Period" "Male" "Female"
## [13] "Dyad" "DyadDistance" "Proximity"
## [16] "DyadResponse" "IDIndividual1" "Intrusion"
## [19] "IntruderID" "SpecialBehaviour" "Audience"
## [22] "AmountAudience" "Context" "AudienceComposition"
## [25] "KinPresent" "Gender" "ToleranceBinomial"
## [28] "predicted_tolerance" "Tol" "Agg"
## [31] "NotApp"
*MAge, FAge,MElo, FElo,
Variations in Age, Rank, IDSI, and IEloDiff and time spent together (amount of trials on a set period/in general) can explain differneces in tolerance (aggression/not approaching?!) rates in a forced proximity box experiment conducted in wild vervet monkeys
The effects may be moderated by /> Seasonality, Male tenure, Female pregnancy
…
Lower variations in age and rank and higher inital social bonds and overall amount of time spent otgether leads to higher rates of tolerance
Check the evolution of the proportion of tolerance per Dyad
Create new dataset with set varialbles
Create Dataset and export as CSV
Pom Elo from a month after experiment maybe since we started experiment as soon as he arrived in Noha
ElO
Paper to read: Factors affecting tolerance persistence after grooming interactions in wild female vervet monkeys, Chlorocebus pygerythrus
therotical background Background: How do individual fogr relationship , how does life history diff impact how these relationships
Goal /Cooperation
Evolution of tolerance, creation of relationship, fwhat foactors change that, food eovolution of social bonds, goal, group living male female realtionhsip, friendhsips
why are picky
Tolerance
* I may consider, in parallel of my hypothesis, to separate the data in *4 seasons* to make a preliminary check of a potential effect of seasonality. Nevertheless the fact that we did not use anywithout tools to mesure the weather and the idea to make a categorization in 4 seasons without considering the actua quite arbitrary. I may do it but with no intention to include this in my scientific report.
l temperature, food quantitiy and other elements related to seasonailty make this categorizationn a categorization where 12 months of data will be separated in 4 categories
#Lines to check unique values in MaleFemaleID to see if they are any problems with it # Unique values in MaleID unique_male_ids <- unique(BexClean$MaleID)
unique_female_ids <- unique(Bex$FemaleID)
What factors influence the rate at which individuals (vervets) learn to tolerate each other in a controlled box experiment?
Ex: The rate at which individuals (vervets) learn to tolerate each other in a box experiment is influenced by social factors (audience, social network, behavior of the partner) and idioyncratic factors (age, rank)
The presence of a higher number of high-ranking individuals in the audience will negatively correlate with the level of tolerance achieved among vervets in the box experiment. This is expected to result in higher frequencies of aggressive behaviors, intrusions, and loss of interest, particularly from lower-ranking individuals.
Vervets tolerance levels in the box experiment will be influenced by their partner’s display of agonistic behaviors. Specifically, partners who exhibit more frequent agonistic behaviors towards their partner will lead to decrease in their motivation to participate in future trials.
During the box experiment, vervet dyads will establish an “optimal” distance for interaction, characterized by a higher frequency of tolerance compared to other distances. This optimal distance is expected to signify that the individuals tolerate each other more effectively at this specific proximity .
The age and rank of individual vervets within the group will influence the success of the trials in the box experiment. Specifically, older and higher-ranking individuals are expected to exhibit lower rates of success compared to dyads consisting of younger and lower-ranked individuals. This decrease in success is anticipated to be associated with a higher frequency of aggressive behaviors displayed by older and higher-ranking individuals towards their partners. (I’m not sure this hypothesis makes sens, I have the feeling age and rank must have an influence but I don’t know how to put it, I will think about it)
Seasonality is expected to impact the motivation of vervet dyads to participate in the box experiment. We hypothesize that dyads will have lower motivation, as indicated by a reduced number of trials, during the summer months compared to the winter months. This difference in motivation is likely influenced by temperature and food availability. To test this hypothesis, we will categorize the data into four seasonal periods, each spanning four months, and analyze whether there is a significant effect of seasonality on the motivation to engage in the trials.
Variables Needed:
DyadResponse (specifically, “aggression” responses) Amountaudience (to measure the number of individuals in the audience) Audience…15 (to identify the names of individuals in the audience for calculating dominance ranks) Elo rating of the individuals based on the ab libitum data collected in IVP (which I have to calculate asap)
Statistical Analysis:Logistic Regression, as it could analyze the influence of high-ranking individuals on the occurrence of aggression in dyad responses. This will help determine whether the presence of high-ranking individuals affects the likelihood of aggression.
Variables Needed:
Statistical Analysis: Logistic Regression as it could be used to assess how the occurrence of aggression in dyad responses is influenced by the partner’s gender-specific agonistic behaviors.
Variables Needed:
Statistical Analysis: generalized Linear Model (GLM) to investigate whether there is an optimal distance that leads to a higher likelihood of tolerance (Tolerance = 1).
Variables Needed:
Statistical Analysis: Logistic Regression Logistic regression can be employed to determine whether the age and rank of individual vervets within dyads have an impact on the likelihood of tolerance (Tolerance = 1).
Variables Needed:
Statistical Analysis:
ANOVA or Kruskal-Wallis Test: Depending on the distribution of your trial data, you can use either ANOVA (if the data are normally distributed) or the Kruskal-Wallis test (for non-normally distributed data) to assess the impact of seasonality on the number of trials. If significant differences are found, you can follow up with post-hoc tests to identify which seasons differ from each other. Please note that the effectiveness of these analyses may depend on the distribution of your data and specific research objectives. You may also consider conducting exploratory data analysis (e.g., visualization) to gain a better understanding of your dataset before performing these analyses. Additionally, if you have specific questions about data preprocessing or variable transformations, feel free to ask for further guidance. –> I took this from ChatGPT, I have to look more into it
REMARKS: So here are a few updates I made in the document. I also planned to send my cleaned data to Radu (the statistician of UNINE) as he was keen to help me find the right test. Of course I will also look again in Bshary’s and Charlotte’s work with the boxes and improve these suggestions that are quite simple for now
Also I still have to clean the last grpahs about male/female aggression as I didn’t finish that yet. I juste wanted to share my hypothesis and ideas for statistics so I can soon go into the “serious” work
Anyway, thank you in advance for your help <3
Michael
But: intro need triangle shape: broad to narrow end wiht research question> tolerance importance > animal reign, actual knowledge/ direction knowledge we need > show how my experiment goes in that way How to adress the gap, answer with research question
Then explain why choosing vervet monkeys, (IVP in methods), sociality, experiments made
Methods
IVP, research area, (goal, house, type people)
Population: groups, dyads, male/female, ranks..
Box material: boxes, remotes, batteries, camera, tripod, corn (no marmelade ;), (water spray, security reason, non agressive way to select individuals and not engage with mokeys when reachrging boxes with corn), pattern, previous distances, tablets, box experiment form
Tablets
(No observers mentionned)
Habituation boxes > individuals trained to recognice boxes, they have differernt levels of habituation
Patterns > appendix, mention similar to habituation, use to recognize box but efficieny depeds of experience)
Selection dyads > assigment from elo rating (different rank), if above average bond no dyad made, if not possible, availibilty of monkey also factor !! Non random can be a problem, think about why and how you selected data We created variations in dyads made by different sex, rank and not above average bonde (calculate bondeness)
Amount corn, do you want to mention it> maybe important Calculate corn during and placement cf paper on corn /food motivation
Corn (daily intake vervet % made from corn, cf site we saw, cf screenshot, comapre paper previousely made an all)
1st dyad trial (BD) > appendix
Videos > details appendix
Finding dyads > appendix
Placement to attract them > meniton if statiscial made on placement corn
Trials (1 session = max 15 trials/in total) (session could be broken in different sub sessions to reach 15 trials max)
If agression > 1m / If 2x tolerance < 1m , also if not approaching > 1m ( if no tolerance increase distance except if intrusion) (borgeaud > expectation fo aggression)
Time of the day > appendix
Territory? > appendix
Amount sessions p day/week, how we chose the moment to follow them >appendix
Problems/ unplanned events: weather, BGE’s, not finding the monkeys (group, dyad or individual), dispersal of males, river crossing, inacessibility (experiments or boxes), low vision (experiments or monkeys),> appendix
(Where do i mention the confounding variables?) > look in litterature, if something that could affect and already reported in papers check, oterhwise exclude “normal life” factors for both monekys and Experimenter
Types of experimental plan
Statistical tests (for each hypothesis)
Analysis
Results
Interpretation
Conclusion
Inter sexual food competition/or tolerance: which type of adaptation
Tolerance
Adaptation to a new context (box experiment setting and repeated encounters with a specific individual)
Evolution of tolerance in competitive context
Social cognition > meaning influence of the group or audience on the choice of an individual of a Dyad
Primate social structure and dynamics
Foundation of tolerance, solving of competition related to food, mechanisms underlying male-female vervet relationships, dyadic interactions
Insights on the evolution of complex decision making in a social context and adaptive mechanisms related to food competition x
In summary, studying dyadic interactions, particularly between males and females in vervet monkeys, not only enhances our understanding of evolutionary origins and adaptive advantages of complex cognition but also provides valuable insights into human social behavior, cooperation, and cognitive foundations shared across species. These studies bridge the gap between animal behavior research and cognitive science, offering interdisciplinary perspectives on the complexities of social interactions and their implications for both animal and human societies.
• Pisor, A. C., & Surbeck, M. (2019). The evolution of intergroup tolerance in nonhuman primates and humans. Evolutionary Anthropology: Issues and ReViews. Advance online publication. https://doi.org/10.1002/evan.21793 (Pisor & Surbeck, 2019)
| Date | Time | Data | Group |
|---|---|---|---|
| 2022-09-27 | 1899-12-31 09:47:50 | Box Experiment | Baie Dankie |
| 2022-09-27 | 1899-12-31 09:50:07 | Box Experiment | Baie Dankie |
| 2022-09-27 | 1899-12-31 09:53:11 | Box Experiment | Baie Dankie |
| 2022-09-27 | 1899-12-31 09:54:28 | Box Experiment | Baie Dankie |
| 2022-09-27 | 1899-12-31 09:55:19 | Box Experiment | Baie Dankie |
| 2022-09-27 | 1899-12-31 09:56:56 | Box Experiment | Baie Dankie |
| GPSS | GPSE | MaleID | FemaleID |
|---|---|---|---|
| -28.010549999999999 | 31.191050000000001 | Nge | Oerw |
| -28.010549999999999 | 31.191050000000001 | Nge | Oerw |
| -28.010549999999999 | 31.191050000000001 | Nge | Oerw |
| -28.010549999999999 | 31.191050000000001 | Nge | Oerw |
| -28.010549999999999 | 31.191050000000001 | Nge | Oerw |
| -28.010549999999999 | 31.191050000000001 | Nge | Oerw |
| Male placement corn | MaleCorn | FemaleCorn | DyadDistance | DyadResponse |
|---|---|---|---|---|
| NA | 3 | NA | 2m | Tolerance |
| NA | 3 | NA | 2m | Tolerance |
| NA | 3 | NA | 1m | Tolerance |
| NA | 3 | NA | 1m | Tolerance |
| NA | 3 | NA | 0m | Tolerance |
| NA | 3 | NA | 0m | Tolerance |
| OtherResponse | Audience | IDIndividual1 | IntruderID |
|---|---|---|---|
| NA | Obse; Oup; Sirk | NA | NA |
| NA | Obse; Oup; Sirk | NA | NA |
| NA | Oup; Sirk | NA | NA |
| NA | Sirk | NA | NA |
| NA | Sey; Sirk | NA | NA |
| NA | Sey; Sirk | NA | NA |
| Remarks |
|---|
| NA |
| NA |
| Nge box did not open because of the battery. Oerw vocalized to MA when he ap to the box to open it. |
| Sey came to the boxes once they were open |
| NA |
| NA |
| Observers | DeviceId |
|---|---|
| Josefien; Michael; Ona; Zonke | {7A4E6639-7387-7648-88EC-7FD27A0F258A} |
| Josefien; Michael; Ona; Zonke | {7A4E6639-7387-7648-88EC-7FD27A0F258A} |
| Josefien; Michael; Ona; Zonke | {7A4E6639-7387-7648-88EC-7FD27A0F258A} |
| Josefien; Michael; Ona; Zonke | {7A4E6639-7387-7648-88EC-7FD27A0F258A} |
| Josefien; Michael; Ona; Zonke | {7A4E6639-7387-7648-88EC-7FD27A0F258A} |
| Josefien; Michael; Ona; Zonke | {7A4E6639-7387-7648-88EC-7FD27A0F258A} |